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SUMMARY 


An essential feature of the interaction phenomena between the 
compression shock and the boundary layer is the coexistence of 
parallel streams of subsonic and supersonic velocities. The com- 
pression disturbance of the shock is propagated ‘‘ahead” of the 
shock through the subsonic portion of the boundary layer, and 
“softening” of the shock takes place. It is the purpose of the 
present paper to demonstrate definitely this essential feature of 
boundary layer-shock interaction by using a simplified model. 
The boundary layer is replaced by a uniform subsonic stream of 
finite width, bounded on one side by the solid wall and on the 
other side by the interface with the uniform supersonic stream of 
infinite extent. The flow fields of two cases are analyzed in de- 
tail: (a) a compression wave in the supersonic stream incident 
upon the subsonic layer, and (b) outgoing compression waves 
generated by a sudden change in the slope of the solid wall. 
Both types of disturbances are assumed te be small, so that lin- 
earization of the differential equations is possible. A general 
result is that the distance of upstream propagation is directly 
proportional to the width of the subsonic layer. Therefore, 
when the boundary layer is turbulent and the subsonic layer is 
extremely thin, there is negligible softening of the shock. 


INTRODUCTION 


= PHENOMENON OF INTERACTION between the 
shock and the boundary layer in supersonic flow 
over a solid surface, as first explicitly demonstrated by 
Ackeret, Feldmann, and Rott,' and then by Liepmann,? 
has aroused great interest among the investigators of 
fluid mechanics. However, because of the great com- 
plexity of the phenomenon that involves both viscous 
heat-conducting effects and the effects of compressi- 
bility, it has not been possible to give a complete ac- 
count based upon the first principles. The above- 
cited authors themselves have tried to elucidate some 
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features of their experimental observations. A more 
conscientious effort on the qualitative relations in the 
interaction of shock and boundary layer is that of 
Lees.* Lagerstrom* introduced the effects of viscosity 
and heat-conduction into the discussion of general 
compressible flow outside of the boundary layer. How- 
ever, in order to make the mathematics involved tract- 
able, he limits himself to the consideration of small dis- 
turbances and is thus able to linearize the differential 
equations. Under this assumption, the velocity at 
every point of the flow field has to be extremely close 
to the free-stream velocity—i.e., the flow is completely 
supersonic. On the other hand, it is the general belief 
of the investigators of the interaction phenomenon that 
one of the essential features is the coexistence of sub- 
sonic velocities in the boundary layer with the super- 
sonic velocities in the main stream. This important 
fact is most clearly demonstrated by Howarth in a re- 
cent paper.’ Howarth considers the problem of effects 
of a disturbance generated in a uniform supersonic 
stream that is bounded on one side by an‘ interface with 
a uniform subsonic stream. Both regions are thus of 
semi-infinite extent. By assuming small disturbances 
and a perfect compressible fluid, he shows that the dis- 
turbances are propagated upstream through the sub- 
sonic region and that there is a region of expansion im- 
mediately behind the point of incidence of a compres- 
sion wave in the interface. 

Howarth’s results are in general agreement with ex- 
perimental observations on the shock and boundary- 
layer interaction. Thus the concept that at least one 
of the essential features of this phenomenon is the co- 
existence of subsonic flow with supersonic flow is 
proved quantitatively. The purpose of the present 
paper is to continue this line of attack and to improve 
Howarth’s model for the real flow. The boundary layer 
is simulated by a uniform subsonic stream of finite 
width bounded on one side by a solid wall and on the 
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other side by the interface with a uniform supersonic 
stream of semi-infinite extent. The fluid is assumed 
to be nonviscous and nonheat-conducting and the dis- 
turbances are assumed to be small. Two cases will be 
considered. First, the case of an incident compression 
wave from the supersonic stream on the interface will 
be studied. Second, the wall is assumed to have a 
break with a small constant slope after the break. For 
both cases, the pressure distribution at the interface 
and at the wall will be calculated. The results can 
then be used to interpret the experiments on the inter- 
action of shock and the boundary layer. 


INCIDENT WAVE 


Let the .-axis be parallel to the undisturbed flows and 
coincide with the undisturbed interface (Fig. 1). For 
y > 0, the flow is supersonic, and quantities there will 
be denoted by the subscript 1. Let the line y = —) 
represent the solid wall. For —b < y < 0, the flow is 
subsonic, and quantities there will be denoted by the 
subscript 2. Furthermore, the pressure, the density, 
and the Mach Number will be denoted by ?, p, and 47, 
respectively. Then 14 > 1 and Mz, < 1. The lin- 
earized differential equation for the disturbance velocity 
potential ¢; in the supersonic stream is 


(My? — 1)(d°%¢,/0x?) — (d%¢,/dy?) = 0 (1) 


The general solution of this equation can be written as 


f(x +N Me _~ ly) + g(x — V My? —ly) (2) 


no = 


fi = 
where f and g are arbitrary functions. If the undis- 
turbed flow is in the positive x-direction, then waves 
the first term on the right of Eq. 
The second term rep- 


represented by 
(2) are the incoming waves. 
resents the outgoing waves originating from the inter- 
face. 

To fix ideas, the form of the function f will now be 


specialized. Let 
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Fic. 1. Compression wave with deflection angle ¢ incident upon 
the interface of parallel subsonic and supersonic streams, 

- / ° - . l . 

fix + VM? —1ly) = Ure 3 M2 — 
e~Ais+VMU—19), for x + VM: — ly>0 (3 

fix + VM, —1y) = 0, forx + 


V M,? —_ Ly <0 


where U is the undisturbed velocity in the supersonic 
stream, and e, 8areconstants. The reason for choosing 
this form for f will be clear presently. By differentiat- 
ing the potential g; with respect to x and y, one obtains 
the x and y components “ and v of the disturbance 
velocity in the supersonic stream. Therefore, the dis- 
turbance velocities of the incoming waves in the x and 
y directions are 


— Uy (1/V Mi? — 1) e- 8 + VMI» 


and 


U,ee ba tv 'y, forx + VMi?-—1ly >0 


For x + V 1/,2 —1y < Othere is no disturbance. By 


making 8 — 0, the incoming waves then degenerate 
into a single compression wave with a deflection angle 
equal to e, and the origin of the x, y plane is the point 
of incidence of this compression wave with the inter- 
face of supersonic and subsonic streams. 


For convenience of calculation, the discontinuous form of the disturbance velocities due to incoming waves call 


be written in the equivalent Fourier-integral form. 


” ae 


Us 


g(x — N M: —ly)= uf )A,(0) sin A(x — *V 
0 4 


where A;(A) and B,(A) are the undetermined Fourier coefficients. 


interest are the disturbance velocities at the interface. 
of Eq. (4). Thus, 


> 0) fs B 
= — s Xz 
<O "sh wae" * 


Then the potential of the outgoing wave should be also written in the form of a Fourier-integral. 


This is possible by noting that 


— sin rzb dx (4 


62+? f 


Thus, 


, 


M2 — 1y) + B,(A) cos Mx — VM? — iy) ban (5 


As will be shown, the quantities of particular 
These can be obtained from Eqs. (3) and (5) with the aid 
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f {XA,(A) cos Ax — AB,(A) sin nx} (6) 


7 ee ee ee 
.o= Oy  f Yar i 08 + a sin ela 





VM; - if {XAi(A) cos Ax — AB,(A) sin EN (7) 
0 


The linearized differential equation for the disturbance velocity potential g: in the subsonic stream is 
(1 — My*)(d? g2/d x2) + (d? w/dv*) = 0 (Ss) 


The appropriate form of solution for the case concerned is 
g = uf cosh XV 1 — AMe?(y + b){A2(A) sin Ax + Ba(d) cos Ax} dr (9) 
0 


where A.(A) and Bs(A) are again the undetermined Fourier coefficients. That this form of g satisfies the differen- 
tial equation can be easily verified. The boundary condition on the wall is 
Voya —p = 0 (10) 
This condition is satisfied also by ¢ specified by Eq. (9). At the interface, y = 0, the disturbance velocities of the 
subsonic side are thus 
Us, q = v. | cosh XV 1 — 1422b{A2(A) cos Ax — AB2(A) sin Ax} dd (11) 


and 


Yoy= -9 = U2V1- ue f sinh XV 1 — 1/22b}AA2(A) sin Ax + AB2(A) cos dx } dd (12) 
: 0 


In the solutions for supersonic stream and subsonic equations from Eq. (13) are 
stream, there are four unknowns A,, Ao, B,, and Bo. Me? poo 
These have to be determined by the conditions at the A,(A) — V2 cosh (AbV 1 — M3*)A2(X) 
interface. Since the disturbances are assumed to be ™ 8 
(15) 


(8? + d?) 


small, the values of quantities at the interface can be ; 
always taken at y = 0, the undisturbed interface. The 
physical conditions to be satisfied at the interface are 
equal static pressure and equal inclination of the flow. 
Since under undisturbed conditions the static pressures ‘ | 
in the supersonic stream must be equal to the static — - : I - = | (16) 
pressure in the subsonic stream, only the change in rV Mi —- e+e 

static pressure due to the disturbance velocities needs Similarly, from Eq. (14) one has 


rV My — 1 


Me? — 
B,() — =f cosh (kb 1 — M2*)B2(d) = 


1 


to be considered. The change in static pressure is y, . 
equal to —pl’u under the assumption of small disturb- 4 ,(,) + ’ .- = dom 1 — M,2)B2(d) = 
ances; thus , VM? - 1 
—p Uj) = +0 = —p2eUo2u2,y = ~0 (13) _ ; : | - e : (17) 
o wae ae | nV My? — 1LA(B* + d’) 
lhe condition for equal flow inclination at the interface 
is, Similarly, V1— MM? . = 
: B\(A) - — === sinh (ABV 1 — M2?)A2(A) = 
ye + 0/ U1 = Vey = ~0/U2 (14) Vay = 1 
€ 
Eqs. (13) and (14) are then the required conditions at = ( = ) (18) 
: rV M2 — 1 \6? +»? 


the interface. 

By substituting Eqs. (6), (7), (11), and (12) into Eggs. (15) to (18) are the four equations for the four 
Eqs. (13) and (14), one can then equate the correspond- unknowns A, Ao, Bi, and Bz. Solution of these equa- 
ing Fourier coefficients on the two sides of these bound- tions then determines the four quantities. For in- 
ary equations. With some reductions, the resulting stance, 
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M.? s ? | ais M,* . @ 9 Ye 
eins :) cosh? (AbV 1 — M,?) + ( - sinh? (AbV 1 — Ao) Aon —— , x 
_ M,?, M,? ae l rN M,? a 1(B? + d?) 
B/ Me? 3 1— Me. pean 
( :) cosh (AbV 1 — AMp*) + V5 sinh (AbV 1 — Mz’) 19 
AME M;? - 1 


Ye 


M.°\? ) 1— M,’ 
| (3 *) cosh? (ABV 1 — Me") + ( ete ) sinh? (0% [= | B,(X) = x 
M a M,? —_ | rn ] —_ M,?(B" + d?) 


1 
M,? 6 jl— M.? . ” 
( .) cosh (ABV 1 — AM") Pr sinh (AbV 1 — M2?) <0) 
M,? A} VAM? - 1 


There are similar expressions for A; and B;. With the Fourier coefficients so determined, the solution of the prob- 


lem is then complete. 

However, the quantities that are of particular interest are the pressure distributions along the interface and 
along the solid wall. If Ap is the change in pressure from the undisturbed value, then Ap = —pUu. At the 
interface, let the change in pressure be denoted by (Ap),;; then Eq. (11) gives 


(Ap); = —pl 2 cosh (AV | Mfo2b) }XAa(X) cos Ax — ABa(A) sin Ax j dd 2} 
Fe 


The integral can be evaluated after the A, and B, from Eqs. (19) and (20) are substituted into the integrand. 
simplify the resultant expression, it is convenient to introduce the following notations: 


0 | (#2) ( =] | (42) 4 (: = ve) 29 
cos My? es. | / Me M2 — 1 oe 
¢ ( a i ey 1 &= 4 
sin - = 2 (92 
m= Ve —1)/ ore M,? — 1 


wee = / ey. i a a 
8 = Nae) / Nae M2 — 1 ai 


with 


and 


Then Eq. (21) can be written as 


(Ap), le 6°*5 nN he 
; — = cos cos Ax + — -sin Ax7pdy\ + — sin? 0 X 
(1/2)p.U;2 xVM —1 2Jo VB? + 82 + 2 f 2 
: (B cos Ax + AX sin Ax)dd Lk. 
/ e —— + =sin@ X 
(8? + d*)[cosh (26V 1 — MW22A) + cos 6] = 2 


[ }A sinh (26V 1 — M,.°A) cos Ax — B sinh (26V 1 — M.A) sin Ax} da ae 
0 (B2 + X2)[cosh (2601 — Me?d) + cos 6] ‘ 


Now the problem can be greatly simplified by making 6 — U—.e., the incoming waves are finally reduced to a 


simple compression wave with a deflection angle «. At this limit, the first integral of Eq. (25) becomes the unit 


step function 1(x)-—-1.e., 
I(x) l, «>0) 
2b 
0, x < Of 


By the substitution \/8 = z, the first part of the second integral is 


f [ B cos \xdX J cos (Bxz)dz 
lim = lim ; 
8079 (8? + )*)[cosh (26V 1 — Mov) + cos 6] a0 7 (2? + 1)[cosh (2h6V 1 — M2282) + cos 6] 


l ft dz T 
1+cos@Jo 22+1 £422(1+cos0) 4cos*(6/2 


Denoting the limiting value of ( Ap); by ( Ap),*, 
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Ap);* te | 9 1g) + Yer, 4 Be of sin (S&/m)dS + 
<i TT ae ; —] cos- (§ sin* sin - ar emanate 
/2)oU2 WM —1 _~™ 2 2 Qr 0 S[cosh S + cos 6] . 
-_ “sinh S cos (SE/m)dS a 
sin 0 ————ene FBT) 
27 0 S [cosh S + cos 6] 
where 
S = 2rV1 — M2? (28) 


t= ex/2bV1 — M2? (29) 


é 
£ is thus the characteristic nondimensional distance parameter. The first integral in Eq. (27) is an odd function 
of &, while the second integral is an even function of —. As shown in the Appendix, for positive values of £, these 


integrals have the following values: 


4 « . ge / . rs) »—- 6 w)t (0/m)§ 
| [ sin(Sé dS - l _— oa x pant Ne J is a e~ \ (30) 
2rJo  S[cosh S + cos 6] S cos?(6/2) Ir sin On =0 lon + 1+ (0/7) on + 1 — (6/r) 
ind 
° : - * \ ‘ x ; O/myé pO /még 
| / sinh S cos (SE m)dS 7 Tenens nes é 4 é GS (31) 
2rJ' S[cosh S + cos A} Irn 9 lon + | + (0/m) In + | (@ r)$ ' 
By substituting the values of these integrals into Eq. where 
27), the pressure distribution at the interface for posi- l2n 4 (0/m) lig 
tive & is, finally, F.(t; 0) = : ——— (35) 
rn=02n + | (0/2) 
Ap),* te ; s — . 
>) U2 = [1 + sin @F,(£;4)], E>0 (32) To determine the pressure distribution along the 
2\o : / Bi . a - , 
sit VM, wall, one obtains first, from Eq. (9), the x-component 
where the function F; (£; 4) is of disturbance velocity mu: at y = —6, the wall. Thus 
‘ 1 \¢c, . c > 
with values of Az and B, from Eqs. (19) and (20), 
, | ie ole ie oo 
F,(E; @) z. (33) _» — — 
Ta axt02n + 1 +. C/o) leyn ng = Us {\Ao(A) cos Ax — ABA) sin Ax} dr 
7 


By observing the proper symmetry of the integrals in The pressure change at the wall (Ap), can then be 
Eq. (27), one can easily obtain the pressure distribution calculated as —p2Ueu, (vy = —b). Now let 8 > 0soas 
to obtain the value (Ap),,.* for a single compression 


for negative values of €: 
wave of deflection «. The procedure of calculation is 


Ap),* ide ae : rae 
an i sin 0F,(E;8), &€ <0 (34) similar to that for the pressure distribution along the 
1 2)aUir VAM —1 interface. The final result is 
S . 2. ££. SB 
F +, cosh — sin ds -, sinh = cos—— ds 
Ap), le | 6/2 2) 1 A 9/2 2 “ (36) 
1/)\p 1,2 JS 5 Toe . + sin = cos ( ) i 
~)pi ly V M,2 -— 112 2\r/ So ~ S[cosh.S + cos 6] 2 2\r//0 S[cosh S + cos 0} 
é <As 


Che.first integral in Eq. (36) is again an odd function of € and the second integral is an even function of £ 
given in the Appendix, for positive values of &, these integrals can be evaluated as follows: 





Ss. Si 
e » cosh > sin dS ~ @/x)k /a)t 
~ [ a T . l A l Cs 1 ) tp (2m + 1) e sci _ _¢é ave 
tJ0 Slcosh S + cos 4} , 6 5 0 ~ a 6 (37) 
2 cos? 5 COS 5 #=0 22n+1l+- 2n+1-- 
~ ~ TT us 
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@ > 
—(0/mjé A rye 


Ss SE 
. sinh - cos — dS | 
9 a » 9 ] \ — (25 + 1)£ é 
2 a (—1)me~ an + 08 fora cecemitnd 
6 ) 4 26 
0 pres 2n+1- 





Sfeosh S + cos 6 , 

” [ rs, =5 
2 

Then, for positive values of £, the pressure distribution can be written as 


* 
(Ap) te oe 
OP) = —— [1 — 2cos : F;(é; 0)], &€ >0 39) 


(1 2)pU2 VM? i 2 


for > 0. 


where 
—[2n +1 + (6/m) iE 


1 @ 
F3(€; 0) = (9p? (40 
~ i 2n+1 + (0/7) 


By utilizing the symmetry properties of the integrals in Eq. (36), one can immediately write down the expression 


for the pressure distribution along the wall for negative values of £: 


(Ap) * be 6 " 
/ = = a 2 cos Fy(é; 0), & « 0 +1) 
1/2)aUi? VM? —-1 - 
where 
| x e [2n +1 — (6/m) ]\é 
Fy(&; 0) = (— 1)” = +2) 
, - 2n + 1 — (0/x) 


Eqs. (32), (34), (39), and (41) then give the formulas for computing the pressure distributions of interest. 
For comparison with experimental observations on the shock-boundary layer interaction, it is of interest to de 


termine the position of the displaced interface between the supersonic and the subsonic flows. 


ment from the undisturbed position be denoted by 7. 


stituting Eqs. (19) and (20) into Eq. (12) so that 


13(-*) inh (2d Vf.) +A So Ms*( a cosh (2\by 1) \} in A\xd\ 


Let the displace- 
The slope dn/dx of the disturbed interface is found by sub- 


3° + A*)icosh (2ADV | A/s*) + cos A} 
yy l— M; s. 5 a 
wh -} sinh (2\bvY | M.?2 By — --+— cosh (2\dV/1 — WV,’ cos Axd\ 
. 8? + 2) {cosh (246V 1 — Me?) + cos 6} 
Introducing &, S and passing to the limit 8 — 0, as above 
—. ~_ , aE ox 
, | ». cosh S sin as " sin as », sinh S cos dS 
dn or) ae T T A T 
ee sin = { ’ : = — COs | 
dx T 2 a i ) a [cosh S++ cos 0 | /0 Ss [cosh S + COS a)\ 2 To S [cosh S + cos 9 | 
It can be shown that, for — > 0, 
" "® cosh S sin(.Sé/m)dS T ; : ; 
I & = 4 = — receoth# [F, &: 4 — F, g; 4H 
wre [cosh S + cos 6] 2(1 + cos A) 
With Eqs. (30) and (31) and /(—€ —I(&), 
dn dx = 2e sin 6 Fi(é; 6), —>O0] 12) 
%e sin 0 F(t: 0), & < OF 


Noting that the displacement vanishes far ahead of the disturbance, integration, 


sults gives 


with respect to £, of these re- 





¢ 


(40) 


ssion 


0 de- 
slace- 
y sub- 


iS hxdd 


13) 


lese re- 
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- l e -(20 + 1 — (0/x) J/g 

: if. ) ie ») a” & s 
[n(&; 6)| = 2e sin 6 >» <<. 2 
20V'1 -- JM,? Ws | 0 





% 44) 
| "heii lhahdipeibaro | 
= 2 sin 0- - ten ta 7 
é\- #\- 6\? 
eer 2n-- 1 —- m+ i+: (2n +1 +~) 
Tv 7 / 7 
Since fi, fe ~ + © as &— O, the interface slope is infinite at the point of incidence of the compression Wave in 
the interface. However, it is easy to show that, where 6 is a positive quantity, 
’ : (en + vat " _ A6 HW , 65) 
dn dy | e —(2n + 1) sinh~— — — cosh ( 
( : = = +e sin @ | 7 Tt T 
dx/: 245 HF oa T 
, n=0 2n + 1)° — (0/r 
Then, the limit as 6 — O is 
dy dy ; A l x 
( ) — ( = —-4e sin a( ) i (45) 
dx ==+0 dx §=—0 ¥ n a 2n me l — (O/ar 
This is a negative quantity. The significance of these calculations will be discussed in a later section, together 


with the numerical results. 


INCLINED WALL 


In the problem of the inclined wall, the solid wall is assumed to have a small inclination eforx > 0. Forx < 0, 
the wall is parallel to the undisturbed streams of supersonic and subsonic velocities (Fig. 2). The x-axis is again 
the undisturbed interface between the semi-infinite uniform supersonic stream of Mach Number M, and the uni- 
form subsonic stream of width 6 and Mach Number J/:. Since now the disturbances all originate from the in- 
clined wall, there are no incoming waves from the supersonic stream. The appropriate disturbance velocity po 
tential y, for the supersonic region is then 


= 

gi = U; / }a,(A) sin A(x —V A? — ly) + d(A) cos Mx — V M2 — ly) }dr (46) 
J 0 

where a, and /; are again the undetermined Fourier coefficients. For the disturbance velocity potential g in the 

subsonic region, one can take 


= - S/ay ») 8 \ 
oe = v| — Z sinh \V 1 Ms*(y + 9)3 ———— SAS —— sin xy dy + 
™V 1 — M2Jo d (eo? + » B*? + 2° 


| cosh V1 — M.7(y + b) { a2(X) sin Ax + bo(X) cos ne} (47) 
0 


Where a, and 6) are the Fourier coefficients to be determined. The y-component 7 of the disturbance velocity at 
the wall, obtained by differentiating the expression in Eq. (47) with respect to y and then setting y = —b, is solely 


due to the first integral. Then, if 8 — 0, one has 


= 0, << 


This means that for x > 0 the wall is inclined at a small angle ¢ from the free stream. Therefore the boundary 
condition at the wall is satisfied by choosing g: as given by Eq. (47). 

To determine the Fourier coefficients a), a2, b;, and b2, the physical conditions at the interface as expressed by 
Eqs. (13) and (14) again have to be used. In other words, the velocity components at the interface from the 
supersonic side and from the subsonic side must be first obtained from the expressions for g; and ge, and then these 
values are substituted into Eqs. (13) and (14). By equating the corresponding Fourier terms, the following equa- 
tions for the unknowns 4), d, };, and b. are obtained: 





Va’ eR i ) 43) 


»y ~ V1 = M6){a,0)] = a. 
ay(A; Me cosh (AV 1 Mb) [a2(d) | Me moe a B+ 
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MOOT nag ee Mz? ¢ sinh (AV 1 — M,% 3 
(A) — 47,2 ©08 — 4Ms’b) [b2(0)] = a7 rV1 — Me \(B? + X?) - 
1— M? - ecosh (AV 1 — Mz) 3 
ay(A) + or sinh (A V1 — Me"b) [bo(A)] = — ; a) 
M;? — 1 TV M2 | (8? + X?) 
1 — M2 - e cosh (AW 1 — Me2) 
by(A) — 2 sinh (A V1 — Web) [a2 (A)] = 7 ( ais : 51) 
M,;? — | TY M;? =a, B2 +4 2 1 


The solution of these equations yields the values for a, and 62 given by the following expressions: 


M,? J RS M,’ ; 
:) cosh?(AbV 1 — Ms?) 4+ Ve , Sinh® (Ab'V 1 — Mo?) jae(A) = 
1” i ae 


M 
é B/Me\  1./Mye — era) 1 — M2") . : = | 
_-— + Sale 1) +— sinh(2\oV 1 — Me? 59 
TY My? — 1(B8? 4 shea ty: — M,?\\An? My — 1s" si = 


ui 


t | M* lL [M? - ACG) + 1— M.") | ; 
rVMP— 182 +29L Me 2Y1— tha) © ee — 1f ONT — ALD] OS 


Similar expressions for a, and 6; can be also computed. 

With the Fourier coefficients so determined, thé desired pressure distributions along the interface and the wall 
can be easily calculated by following exactly the procedure outlined in the problem of an incident wave. No de- 
For the pressure distribution along the interface 


| (#2) OV Me) + sinh? OV vi) 
WV cosh? (AbV 1 — Ae") re , sinh AOV 1 — Ae") [b2(A) = 


tails of calculation will be given here but only the final results. 


(a >’), 
( Ap),* Je 1 cot (6/2) _ l 6 A\ ’ 
ae l : Fy(&; 0) cos 3 cot ) Fale; 8 |, —>0 D4 
(1/2)p.l; VM? 1 2sin (0/2) 2 2 » 
and 
( Ap) Ze ] 6 oo l cot (4/2 ; _ 
—— a cos 3 cot’ Fl ' 6) + . F &: 0) |, E < 0 ) 
(1/2)aiU;2 Wy? — 1L2 2 2 2 sin (0/2 


The pressure distribution along the wall is given by (8 — 0 


(Ap) o* Ze 6. 
= I + 2 cot rie) | E>0 Hf 
(1 2)p,l 5” \ AM, a | Z 
and 
Ap) wo* 2e 6 iad 
Phe cae 2 cot rs F,(§;0), —€< 0 oi 
(1/2)p,U 7 VM? - 1 - 


In Eqs. (54) to (57), the functions F\, Fo, F3, and Fare functions defined by Eqs. (33), (35), (40), and (42). The 
variable @ is connected with the Mach Numbers .J/; and W/2 as specified by Eqs. (22) to (24). These equations 
then supply the desired information for the problem of an inclined wall. The numerical results and the discus- 


sions of the results will be presented in the next section. 


between @ and 1/;, 2 are given by Fig. 3. It is seen 
that, if the representative Mach Numbers for the 
boundary-layer flow and the free stream are taken t 
be 0.8 and 2, then the corresponding value of @ in the 


NUMERICAL RESULTS AND DISCUSSION 


The parameter @ occurs in all expressions for the 
pressure distributions. As shown by Eq. (22), it is 


connected with the Mach Numbers ./; and M2 and can 
be taken to have the range from 0 to 7, When 0 <_ flow model is approximately 37/4. This fact should be 


M. < 1, the value of @ increases to * when the ratio kept in mind when the theoretical results are compared 
M,/M, tends to zero. When M, > 1 but finite, @ tends with the experiments on the shock and boundary-layer 
to 0 when M; tends to 1. The general relationships interactions. 
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The series for the functions /\, 2, F3, and Fy as given 
by Eqs. (33), (35), (40), and (42) are suitable for numer- 
ical computation when & is large, since the convergence 
is then extremely rapid. For small values of £, modi- 
fied series as given in the Appendix is more appropriate. 
F, and Fy have a logarithmic infinity at £ = 0, but Fs; 
and F; are finite at & = 0. Table 1 gives the tabulated 
values for these functions. 


The results of computation for the case of an incom- 
ing compression wave are shown in Figs. 4, 5, 6, and 7. 
In these graphs, the abscissas are the characteristic 
nondimensional length & along the streams defined by 
Eq. (29). The physical distance x for any value of & 
is thus proportional to bV 1 — \,°. The width b of 
the undisturbed subsonic stream is then the measure of 
length in the problems concerned. The ordinates of 


these graphs are « 
[Ap/(1 2)p,U,"| (Ze \ ‘M2 — 1) 


Since, under the general assumption of small disturb- 
ances, a compression wave of ¢€ deflection is associated 
with a pressure rise equal to 


(1 2) p, U1?(2e \ M,? — 1) 


the ordinates of these graphs are then the ratio of pres- 
sure rise to the pressure rise caused by the simple in- 
coming compression wave. It is then immediately 
seen that the pressure rise far downstream of the point 
of incidence of the wave on the interface is twice as 
large as the pressure rise for a simple compression wave. 
When an incoming compression wave is incident upon 
i solid wall, the compression wave is reflected as a com- 
pression wave, and the pressure rise after the reflection 
is twice as large as the pressure rise before the reflection. 
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Fic. 3. The Mach Number parameter @ as function of Mach 


Numbers j/, and ./, in the supersonic and the subsonic streams, 
respectively 


ing compression wave is reflected as if by a solid wall 
without the subsonic stream. 

Another feature of these graphs is that the pressure 
distribution along the interface has a peak at the point 


lherefore, except for the local interaction, the incom- of incidence (£ = 0) of the compression wave. There- 
TABLE 1 
a=) @ = r/4 
t Ff i F; F, F, F. F F, 

0 x x 0.2512 0.2512 ro x 0.1940 0.3496 
0.05 0.5927 0.5927 0.2431 0.2431 0.5090 0.6761 0.1858 0.3416 
0.10 0.4825 0.4825 0.2343 0.2343 0.4010 0.5716 0.1771 0.3328 
0.25 0.3369 0.3369 0.2107 0.2107 0.2616 0.4335 0.1542 0.3079 
0.40 0.2629 0.2629 0.1887 0.1887 0.1942 0.3611 0.1334 0.2835 
0.50 0.2275 0.2275 0.1735 0.1735 9.1627 0.3251 0.1204 0. 2674 
0.75 0.1649 0.1649 0.1391 0.1391 0.1104 0.2580 0.0907 0.2291 
1.00 0.1246 0.1246 0.1121 0.1121 0.0783 0.2094 0.0693 0.1932 
1.75 0.0561 0.0561 0.0548 0.0548 0.0290 0.1157 0.0282 0.1135 
2.50 0.0266 0.0266 0.0260 0.0260 0.0118 0.0656 0.0112 0.0646 
1.00 0.0059 0.0059 0.0018 0.0213 0.0212 

f = wr/2 . — 8 ee f= 3r 4 -— — - 
E F F, F; F; F, F, F; Fy, 

0 x o 0.1569 0.5538 oo « 0.1312 1.1776 
0.05 0.4520 0.9461 0.1486 0.5457 0.4092 1.6056 0.1229 1.1698 
0.10 0.3459 0.8278 0.1401 0.5340 0.3050 1.4895 0.1144 1.1608 
0.25 0.2111 0.6803 0.1175 0.5147 0.1762 1.3252 0.0930 1.1338 
0.40 0.1498 0. 5869 0.0988 0.4836 0.1189 1.2321 0.0756 1.1050 
0.50 0.1217 0.5444 0.0872 0.4654 0.0935 1.1802 0.0651 1.0817 
0.75 0.0771 0.4620 0.0631 0.4202 0.0552 1.0872 0.0445 1.0323 
1.00 0.0513 0.3993 0.0448 0.3766 0.0345 1.0125 0.0298 0.9777 
1.75 0.0156 0.2676 0.0152 0.2640 0.0087 0.8256 0.0084 0.8200 
2.50 0.0054 0.1821 0.0049 0.1815 0.0025 0.6820 0.0023 0.6804 
+00 0.0006 0.0864 0.0864 0.0002 0.4688 0.4688 
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Fic. 4. Pressure increment (Ap) * for the case of incident wave 
plotted against the distance x from the point of incidence of the 
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Fic. 5. Pressure increment (Ap) * for the case of incident wave 
plotted against the distance x from the point of incidence of the 
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Fic. 6. Pressure increment ( Ap) * for the case of incident wave 
plotted against the distance x from the point of incidence of the 
wave upon the interface. @ = 7/2. 


fore before the point of incidence the flow is compressed, 
but after the point of incidence the flow is expanded. 
In the supersonic stream then, there is a series of com- 
pression wavelets sloping downstream, ahead of the in- 
coming compression wave. After the incident com- 
pression wave, there is an expansion region. The pres- 
sure distribution along the wall has no peak and is a 
smooth rising curve. Most of the pressure rise along 
the wall occurs, however, before the point of incidence 


of the compression wave. In fact for all points ahead 
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Fic. 7. Pressure increment ( Ap) * for the case of incident wave 

plotted against the distance x from the point of incidence of the 


wave upon the interface. @ = 37/4. 














Fic. 8. Incident wave case, the displacement 7» of the inter- 
face (or increment of the thickness of the subsonic stream) plotted 
against the distance x from the point of incidence of wave upon the 
interface. 


of the point of incidence, the wall pressure is higher than 
the pressure at the interface. The region where there 
is strong interaction between the supersonic stream 
and the subsonic stream is of the order of tens of widths 
b of the subsonic stream. The length of this region is 
increased as 6 approaches 7. 

Fig. 8 shows the results of computations for the inter- 
face displacement. A noticeable feature of these curves 
is the rapid change from an infinite slope at & = 0 to 
small values at & # 0. This is especially marked for 
~ > 0. Eq. (45) shows that the slope for positive x is 
less than the slope for negative x. Therefore, as far as 
the subsonic flow is concerned, the corner at the point 
of incidence of the compression wave is a compression 
corner. This result agrees with computations for 
(Ap);*, inasmuch as it implies compression upstream 
of & = O, an infinite rise in pressure at & = 0, and ex- 
pansion downstream of § = 0. It is also seen that the 
displacement increases as @ increases. Recalling that 
6— ras M2/M, — 0, for 0 < Ms < 1, this means the 
displacement increases as J/, decreases—i.e., as the 
disturbance propagates upstream with greater strength 
(see Figs. 4-7), as would be expected. 
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One further comparison with the pressure calculations 
may be made. For & < 0, the wall pressure exceeds 
interface pressure except in an extremely small region 
about the singularity. The pressure gradient is such, 
then, as to be consistent with an increasing interface 
slope as § — 0. For £ > O the gradient is from the 
interface to the wall, consistent with a decreasing slope 
as £ increases. 

When a shock occurs*over a laminar boundary layer, 
the case of the so-called \-shock (Fig. 9), experiments 
show that there is a compression region in the free 
supersonic stream ahead of the shock and an expansion 
region after the shock. Furthermore, the pressure 
rise along the wall is known to be continuous. These 
features of the A-shock are then completely demon- 
strated by the incident wave discussed above. 
It is evident that the analysis presented is not a quanti- 
tative explanation of \-shock, since the shock is a rather 
strong disturbance, while small disturbances and linear- 
ization of the differential equation are the basis of the 
present analysis. Furthermore, in the analysis, the 
viscous and the heat-conducting effects are completely 
Perhaps a more important omission is the 
The theoreti- 


case 


neglected. 
velocity gradient in the boundary layer. 
cal calculation assumes uniform supersonic and sub- 
sonic streams, and therefore the vorticity is concen- 
trated at the interface instead of distributed throughout 
the subsonic layer, as is the case for the boundary layer. 
However, the close qualitative agreement between the 
theory and the experimental observations on the A- 
shock indicates the importance of the effects of inter- 
action between the supersonic and the subsonic region 
in the shock-boundary-layer phenomenon. 

One might now ask: Why is the interaction between 
the shock and a turbulent boundary layer different 
from that of the \-shock? In the case of a turbulent 
boundary layer, one does not observe the extended 
compression region before the shock and the expansion 
If the theory presented is of 
any real value, The 
authors believe the solution of this paradox lies in the 
turbulent 


region after the shock. 
it must explain this case also. 
thickness of the subsonic layer in case of 
boundary layer. Because of the much steeper velocity 
gradient resulting from the extremely intensive turbu- 
lent flow exchange, the subsonic part of the turbulent 
boundary layer is generally extremely thin. For usual 
(dimensions the subsonic layer in a turbulent boundary 
layer is only of the order of hundreths of an inch thick, 
while the subsonic layer in a laminar boundary layer is 
of the order of tenths of an inch thick. Therefore, the 
characteristic interaction effect between the subsonic 
and the supersonic flows, if it occurs, will be limited to 
a region of only a few tenths of an inch for the turbulent 
boundary layer. This would be difficult to observe, 
especially if, as is usually the case, the shock is not 
Stationary but oscillates slightly in position. Further- 
nore, the large velocity gradient in the turbulent layer 
will certainly modify the flow considerably from that 
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Flow pattern of a d-shock. Interaction of shock and 
laminar boundary layer. 
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Pressure increment (Ap)* for case of inclined wall 
= 


Fic. 10. 
plotted against the distance x from the break in the wall. 
m/4. 
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Fic. 11. Pressure increment (Ap)* for case of inclined wall 


plotted against the distance x from the break in the wall. @ = 


w/2. 
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Fic. 12. Pressure increment (4p)* for case of inclined wall 


plotted against the distance x from the break in the wall. @ = 


3/4. 
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Fic. 13. 
airfoil. 


calculated theoretically, both because of the strong in 
fluence of vorticity and because of the effects of vis 
cosity and heat-conductance. Therefore, while the 
theory presented cannot give a clear explanation for the 
interaction between shock and the turbulent boundary 
layer, the theory and the experiments are not contradic- 
tory. 

The numerical results for the second problem of the 
inclined wall are presented in Figs. 10, 11, and 12. 
Here the general features of the pressure distributions 
are the same as for the problem of the incident wave. 
Only now the role of the wall and the interface is 


changed. The pressure distribution along the wall 
now has a peak at the starting point of the inclination, 
while the pressure at the interface rises steadily. At 


distances far downstream of the break in the wall, the 
pressure is equal to that due to a wedge of vertex e, the 
angle of inclination of the wall. This is, of course, ex- 
pected. But the important result is the rather large 
compression ahead of the change in flow direction 
caused by the inclination. One can compare these cal- 
culated results with the measured pressure distribution 
over the upper surface near the trailing edge of an air- 
foil in a supersonic stream (Fig. 13). At the trailing 
edge, the flow has to turn through an angle so that com- 
patibility conditions for the flows along the upper sur- 
face and the flow along the lower surface can be satis- 
fied. This change in the direction of flow is the same 
as that caused by the inclination of wall shown in Fig. 2. 
Therefore, from the results of the calculations, one can 
expect a rise of pressure over the surface of airfoil ahead 
of the trailing edge above that predicted by a theory 
without taking into account the effects of the boundary 
layer. This is actually observed.* Of course, when 
the boundary layer near the trailing edge is turbulent 
instead of laminar, the subsonic layer will be so thin as 
to cause negligible forward propagation of pressure. 
Then the pressure rise before the trailing edge will be 
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Flow pattern near the trailing edge of a supersonic 
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greatly reduced, and the simple classical theory for 
pressure distribution over a supersonic airfoil will be 
more accurate. 
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Appendix 


(I) EVALUATION OF INTEGRALS 


The infinite integrals in the expressions for the pres- 
sure distributions can be evaluated by contour integra- 


tion. For instance, 


SE Sé 
: sin dS . i sin — ds 


l / T _ Jf / T 
2r Jo Sl[cosh S + cos 6] 42 /-2 S[cosh S + cos 6 
l f és ds 
$a » S[cosh S + cosé 


Now consider S as a complex variable and take the 
contour as (1) the real axis indented at the origin, and 
(2) semicircle in the upper-half of the complex plane 
with infinite radius. It is easy to show that the con- 
tribution to the contour integral along the semicircle is 


zero. The contour encloses the following simple poles 
S = 4[(22 + 1)r = 6], nu = 1, 2, 3, 
Thus, 


7 sin — ds | 
- T — +. 
es 2L2(1 + cos 8) 
S[cosh S + cos @] 
x Pe 2n +1 + @/w)Jé 
A a once all 
n=01[(2n + 1)r + 6] sinh ¢ [(2n + 1)x + 4] 
ry ela +1— (0/xr)]é ‘c 
n=0 1[(2n + 1)r — 0] sinhz [(2n + 1)x — 4] 





(Formula continued at top of facing page) 
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l 1 “2 


—(2n + 1)& x 
S cos? (0/2) 


: - e 
2r sin On =0 


{ e” wyt e @ wt \ 


12n + 1+ (0/r) 


i oe 
2n + 1 — (0/r)f 


This is the result given in Eq. (30). 

The other integrals can be evaluated in a similar 
way. The results are already given in Eqs. (31), (37), 
and (38). 


II) CONTINUITY OF WALL PRESSURE 


From the computed results it appears that the wall 
pressure is continuous for the incident wave case and 
the interface pressure is continuous for the inclined 
wall case. From Eggs. (39), (41) and (50), (52), it is 
seen that this condition implies that 


2) sec (0/2) 


(A) 


F3(0; 6) — F (0; 4) i (1 
[his may be verified in several ways. The simplest 


seems to be the following: 


IF LS _iyeenite d+ 0k o 
ag rTn=0 
Le li+@ r)é 
7 7 1l+ e~% 
dF, pe - Ck 
dé a ris" 


Now F; (om; @) = F,( 0 5 @) = 0, so 


, | "a elt + O/o)F 4 Ql! — O/a)hs 
£;0) + Fa(é; 0) = — ee dt 
Jt Il+e 


~ 
m 


and 


cosh (6/7)t 


oe & 
/ dt 
TS 0 cosh ¢ 
1 © cosh (6/7)t 
dt 
27 © cosh ft 


Now consider ¢ as a complex quantity and integrate 
cosh (6/7)t/cosh ¢ around a rectangular contour (1) 
along the real axis from (—R, 0) to (R, 0), (2) normal 
to the real axis along (R, 0) to (R, 7) and (—R, 7m) to 
(—R, 0), and (3) parallel to the real axis from (R, 7) 
to(—R, r). The only pole enclosed is at (¢ + is) = 
1(r/2) and the residue is —i cos (6/2). 


By making R—> 
», Eq. (A) can be shown to be true. 


F3(0; 0) + F4(0; 6) = 
Finally, 


F,(0; 0) + F,(0; 0) = 


SERIES FOR COMPUTING F-FUNCTIONS AT SMALL 
VALUES OF ARGUMENT 


ITT) 


The series for Fi, Fo, F3, and Fy as given by Egs. 
(33), (35), (40), and (42) are suitable for numerical 
computation only when |é| is large. When |& is small, 
the convergence of these series is slow and computation 
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tedious. . On the other hand, the physical problem re- 
quires the values of these functions at small values of 
£. The following series for these functions are thus 
more convenient. 

For & > 0 


+ 1+ (0/w) 


inc 
F\(£; 0) = : nae 
02 + 1 + (6/r) 


Thus 
dF, ae oes  (0/e) hi 
=—- 2 ew it! es 
dé 7TH 0 
— {1 (0/m) jé (e 2) 5 
n=0 
Le + Ok 
- xr li-—e™* 
Since F, (©; 0) = 0, 
* ox [1 + (0/w)]t 
] e 
F,(é: 9) = / =a at 
TJS l—e~- 
Now let 
ee" = S 
Then 
F l f dS 
BJ i FA =- SS 
where » = e~* and k = (1/2) — (0/2). For @ such 


» 
thatO <@< 7,0<k < 1/2. Also,O <7 < 1. 


Now 


1 f’{a—S)4+ S]ds 
i ad ae ee 
27 Jo S“(1 -— S$ 


) 
L (S42 dS 
Sete S*'Qaete 1-85 


(1 — k) 


Thus 


— [S°-” log (1 — S)]B + 


"log (1 — S) 
a-» fee 


The remaining integration may be effected by a further 
Let 1 — Then one has to 


@ —p 
2 (-9+( ) x 
nn = () n 


2rF; = 


ds 


’ 
S=¢e 
ws = Ss 


change of variable. 
evaluate 


I log z 
- = dz = 
@—n) (1 — 2) 


2" log 2dz 


—)"tt- iy} 
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—k\. ; . : a form that is convenient for evaluation :< all v: 
wit ( : ) a a - ° m that is convenient for evaluation at small values 
oO é. 
ait a oe ae In similar fashion modified series for F2, F;, and F, 
= may be determined. Then, 
n n! , ; 
a oe : eo I 
Finally, 2rF, = 7" KS — log (1 aden — k) X 
' ll —k f 
, +9 l 
2 F, = 7" — log (1 - nt — P ae ce 
Juels YY (-1)" +1 ond 
s. / n —k jc _ n)"*? * . nn+1—k 
Zu (—)}) oii = log (In) — 
n + : , 
where n, kare asabove. 22>, 27F3 have the same forms 
l la —»*t - 1}! as 2rF,, 2rF,, respectively, except that k =(1/2) + 
Wa =P EP (0/2m). 
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Theory and Tests on the Plastic Stability of 
Plates and Shells’ 


P. P. BISLAARD?t 
Cornell University 


SUMMARY 


Buckling of plates and shells superimposes an additional state 
of stress on the original one. Assuming that plastic deformation 
is governed only by the amount of elastic shearing energy at the 
point in question, it is shown that the assumption of ‘‘plastic de- 
’ leads to smaller buckling stresses than that of “plastic 
“plastic deformation,’’ it is 


formation’ 
flow.”’ Using the assumption of 
found that the theory is in good agreement with tests by Koll- 
brunner and by the N.A.C.A. It is found that the ratio between 
plastic and elastic buckling stress is highest when the proportions 
of the deviator components, as superimposed by buckling, differ 
most from the proportions of the original deviator components. 
The stress-strain relations are derived by writing the infinitely 
small excess strains as total differeatials and computing the par- 
tial derivatives of the strains with respect to the stresses. The 
differential equation for plate buckling is derived, and results 
of its application to several kinds of loading and boundary con- 
ditions are tabulated. The theory for the buckling of shells is 
also given. The author presets a comparison of his theory with 
the theories of Ilyushin, Stowell, and Handelman and Prager. 


INTRODUCTION 


angen THEORY of the plastic stability of plates 
was given by the author in 1938. It appears, 
from several American publications after World War II 
which deal omly with special cases, that this theory is 
not well known in America. Its results are in accord 
with tests by Kollbrunner,* especially with regard to 
the buckle length. They are also in accord with tests 
conducted by the N.A.C.A.4!~* 


FLow LINES IN WIDE PLATES 


In the author’s theory it is assumed that plastic def- 
ormations are governed by the elastic shearing energy. 
Previously it was thought (for example, by Bierett! 
and others) that with inhomogeneous stress distribu- 
tion the yield stress depends not only on the state of 
stress at the point involved but also on the stress dis- 


Presented at the Session on Structures—Static Loads, Seven- 
Meeting, I.A.S., New York, January 24-27, 





teenth Annual 
1949, 

*The author wishes to express his gratitude to Prof. N. J. 
Hoff, Polytechnic Institute of Brooklyn, and S. Levy, National 
Bureau of Standards, for their kind assistance in the preparation 
of preprints and slides for his lecture during his trip to Europe 
and to Mr. Levy for his willingness to be prepared to give the lec- 
ture in case the author had not returned in time. He is also 
greatly indebted to Mr. Levy for linguistic corrections in this 
final paper. 

+ Professor of Advanced Problems in Applied Mechanics, 
Delft, Holland; now Associate Professor, College of Engineering, 
Cornell University. 


tribution in the adjacent material, both in the elastic 
domain. It is shown by the author’s theory of local 
plastic deformations’ that plastic deformation is al- 
ways governed by the amount of elastic shearing energy 
at the proper point and that apparent raising of the 
yield stress at points of stress concentration is caused 
by a change of the state of stress in consequence of minor 
plastic deformations.** In this way the direction of 
flow lines in wide steel plates® and thick bars® was ex- 
plained as was also the difference between upper and 
lower yield stress in thick steel bars. In wide plates 
subjected to pure tension or pure compression, these 
flow lines form angles of 55° with the direction of the 
force. 


The latter may be easily explained as follows. If, 
for example, in a plate that is sufficiently buckling- 
proof a flow line should occur perpendicularly to the 
direction of the compressive force (Fig. 1), the flow in 
the longitudinal direction (x-direction) would cause an 
expansion of the material within the flow lines in the 
lateral direction (y-direction). This would, however, 
be impeded by the adjacent elastic material, so that 
excess stresses a,’ would be superimposed. According 
to Fig. 2, representing all states of stress situated at 
the yield point, this would involve a shifting of the 
representative point of the state of stress from o, = 
os, ty = 0 too, = of + o;’, o, = 2,’, so that the ap- 
pearance of such a flow line would require a higher stress 
than the yield stress o, for pure compression. Only 
when the flow lines embrace an angle of 55° with the di- 
rection of the compressive force does the plastic strain 
in the longitudinal direction of the flow lines happen to 
be zero, so that the flowing in these flow lines is not im- 
peded by the adjacent elastic parts of the plate. Thus, 
in this case, unimpeded flow at the yield stress for pure 
compression is possible, so that these dow lines require 
the smallest force. 

' that local plastic deformations 


which find expression in bands of 
submarine 


It was shown also’ 
of the earth’s crust- 
negative gravity anomalies, 
ridges and troughs, isles series, and chains of mountains 

—have arisen from flow lines in the earth's crust which 
are governed by the same laws as those in a steel plate. 
Later investigations of Vening Meinesz*® and the au- 
thor* confirmed that the entire system of linear pat- 
terns of the earth’s crust can be explained by the au- 
thor’s theory of local plastic deformations. 


geosynclines, 
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Cylindrical buckling of a wide plate 


BUCKLING OF A WIDE PLATE 


In order to make plausible this explanation of flow 
lines in the earth’s crust, it had to be proved that 
buckling of the earth’s crust, as a spherical shell elas- 
tically supported by the substratum, does not occur 
before the yield stress is reached.’ With the then 
existing theories on plastic buckling of steel plates in 
structural engineering (as given by Bleich, Chwalla, 
and Ro&-Eichinger), this was not possible, because ac 
cording to these theories the buckling modulus of a 
plate, under a compressive stress equal to the yield 
stress, is zero. The author can show, however, that 
this must be wrong. With a wide plate, for example, 
which is subjected to compression in the x-direction 
(Fig. 3), the proper buckling will not cause a strain in 
the y-direction. If, however, with buckling, at the con- 
cave side of the plate only a pure compressive stress 
equal to the yield stress a, should occur, the shortening 
in the x-direction would cause there an expansion of the 
y-direction. In the same way, as in the case of flow 
lines perpendicular to the force, however, the condition 
that this expansion should be zero will cause excess 
stresses a,’ and thus also excess stresses o,’ (Fig. 2). 
Moreover, it can be shown that in this case the entire 
plate behaves plastically, by which the excess boundary 
6 which counter- 


/ 


stresses a,’ cause inner moments /’o, 
act the buckling of the plate. Consequently, at the 
yield stress the buckling modulus of a plate is not zero 
as is that of a bar. Apparently this is caused by the 
fact that, with buckling of a plate, another state of 
stress with another ratio of the stress components is 
superimposed on the initial state of stress, and with 
changing ratio of the stress components, or rather of the 
deviator components, the material no longer behaves 
quasi-isotropically, as in the case for which the ratio 


of these components 1s constant. 


BEHAVIOR OF A METAL UNDER VARYING 


STATES OF STRESS 


PLASTIC 


Thus it was necessary to find a more exact theory of 
the plastic stability of plates and shells. A sound base 


for such a theory was the author's theory of local plastic 


deformations,® in which it was shown that from the very 
beginning plastic deformations are governed by the amount 
of elastic shearing energy at the proper point, in contradic- 
tion with former conceptions (Lode!') in which it was 
thought that in the beginning the yielding in layers 
in the flow lines depended on the maximum shearing 
stress. 

Plastic deformations occur by sliding in the crystal- 
lographic planes of the crystallites. Consequently, 
hydrostatic pressure, causing no shearing stresses any- 
where, will not influence the plastic deformations. 


The stress tensor 


\% "xy Tzz( 
) Tyr Sy ra l 
To Tae Os 
may be split into a hydrostatic pressure o = (¢, + 
og, + ¢,)/3 and into the so-called stress deviator 
\ Try Tre 
s i 5) 
Tyz oy ro 2 
(res Ty &G; 
in which, consequently, ¢, = o, — o, ¢, = o, — o, and 


¢, = o, — o. Thus, only the deviator influences the 
plastic deformations. 


strains, being components of the elastic strain deviatot. 


Furthermore, it causes elastic 
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Representative vectors of strain and stress deviators with constant ratios of the deviator components. 
Representative vectors of strain and stress deviators with changing ratios of the deviator components. 
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Fic. 6 (right). 


Comparison of plastic deformation (I'n,), plastic flow (2), and a new possibility (To). 


€xe = Gz, 2G, Eve = &,/2G ( 
(1/2) Yrye = Try/2G, etc. | (4) 


so that all components are proportional to those of the 
stress deviator. According to tests of Lode,'! Ro&- 
Eichinger,'? and the author,” with constant ratio of the 
deviator components the plastic deformation takes 
place practically quasi-isotropically, so that the com- 
ponents of the plastic strain deviator are proportional 
to those of the stress deviator.* The plastic strain 
deviator is, of course, identical to the plastic strain 
tensor, since plastic deformations do not cause a change 
of volume. 

The elastic shearing energy 
plastic deformations, equals the elastic strain energy 
due to the stress deviator, so that 


Ve - (1 2) (Frere + Tyeve + O2€z¢ + Try Vrve i 


TyzVuze + TerVere) = (1/4G) [6,” + a,” + 8, + 
2(T2y? + Ty: + Tz2°) | = $,°/4G 


Ve, Which governs the 


(5) 


At the yield stress, where V, is constant, consequently, 


Si = Va? +4,? + 6.2 + (tr? + ty22 + tx”) (6) 


is constant. 

The stress and strain deviators may be indicated by 
representative vectors Ig and 7», of which the com- 
ponents are the components of these deviators. Hence, 
if considering deformations at the yield stress with S, in 
Eq. (6) constant, the length of the stress vector and that 
of vector 7, of the elastic deformations are also con- 
stant, while with constant ratio of the deviator com- 
ponents the direction of all vectors, including that of 


the plastic strain deviator, will be the same. In that 


* This means that the deformation is the same as with a shear- 
ing strain in all planes which is proportional to the shearing 
Divergence occurs, however, in the direc- 
tion of a sliding along the planes of maximum shearing stress. 
This may be explained by the fact that the crystallites have a 
certain yield stress themselves, so that no sliding takes place 
along planes with a too low shearing stress. 


Stress in these planes. 


case the vectors of the plastic and elastic deformation 
[being, for example, 7), = AB (Fig. 4) and 7), = BC, 
respectively] will become, after further deformation, 
To’ = AB’ and Ty, = B’C’ = BC, respectively, while 
the vector of the stress deviator, being placed at the end 
of the total (i.e., elastic plus plastic) deformation course, 
moves from C to C’, its direction and length remaining 
constant. 


however, the direction of the total deformation 
course changes, as is the case with plastic buckling of 
plates and shells, so that, after having a direction AC 
(Fig. 5) while AB and BC represent the plastic and elas- 
tic deformations, respectively, it proceeds in a direction 
CD, several possibilities exist for the course of the excess 
plastic deformation. The excess deformation CD by 
buckling has to be assumed infinitely small. In the 
author’s theory,'*'* it was assumed for several reasons 
that at the representative point D of the total deforma- 
tion the situation is the same as with a straight defor- 
mation course AD, so that as the length of vector /D of 
the elastic deformations remains constant and equal to 
BC, the excess plastic deformation is represented by 
BE, while the stress deviator is given by ['y, at D. So 
a complete memory was attributed to the material, as 
in a so-called Hencky and in accordance with 
“plastic deformation.” 


body 


what is called 


Hohenemser and Prager'® executed tests at the yield 
stress with steel tubes that were alternately subjected 
to pure tension and torsion. They found that pre- 
ceding free plastic deformations 
during which the ratio of the deviator components does 
not change—are forgotten. In the present case this 
means that the free plastic deformation AB is forgotten 
and that the continuation of the vector of the elastic 
deformation, FD= BC, passes through B, the excess 
plastic deformation being given by BF and the stress 
deviator by min D. As AB is here the total preceding 
plastic deformation, also a Prandtl-Reuss body, having 


i.e., deformations 
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no memory at all and showing so-called ‘‘plastic flow,”’ 
it will behave in the same way. 

The more general assumption of a Hencky body, 
which includes here the Prandtl-Reuss body as a special 
case, was used in the theory. This was done for the 
following reasons. With deformations below the yield 
stress, with which one has to do with buckling, it is 
not sure that, even with free deformation, the ratio 
of the deviator components does not change locally in 
connection with the yielding of the differently orien- 
tated crystallites, one after another. Consequently, 
one is not sure that in the domain of strain hardening 
the preceding plastic deformation AB will not be re- 
membered. Moreover, by taking it into account, one 
is sure not to overestimate the buckling stresses of 
plates, which according to the author's theory were 
nevertheless much higher than according to the exist- 
ing theories of Bleich, Chwalla, and Ro&-Eichinger. 
For with a Hencky body (stress vector I), in D) the 
excess stress deviator, which after an excess total de- 
formation course CD is given by A,Io, is smaller than 
with a Prandtl-Reuss body (stress vector Ig: in D), as 
may be seen directly in Fig. 5. 

As, however, a loss of memory amounts only to a 
neglect of the initial plastic deformation course AB, it 
may eventually be taken into account in the author's 
formulas by equating his value e to zero. This was done 
in his earlier work'* in order to compare his results with 
tests by Kollbrunner and Chase. From Kollbrunner’s 
tests with angles’® he concluded, however, that with 
aluminum all preceding plastic deformations should be 
taken into account. This view was confirmed by the 
good agreement between his theory and the results of 
Kollbrunner’s®* and Stiissi’s*' later tests with plates. 


On the other hand, the author called attention to the 
fact that a complete behavior as a Hencky body will 
generally not be possible."*'’ If, for example, after a 
deformation course AC, Fig. 6, the total deformation 
follows the are CD, the plastic deformation would 
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follow the are BE. The work done by the stress devia- 
tor with a small increment of strain may be written as 


follows: 


dV, = 6,dé, + &,dé, + dé, + Trd¥r + 
Ty2h Vys + Tax 22 (7 


being the inner product of the vectors of stress deviator 
and strain increment. In this case, where the former 
remains perpendicular to the latter, the stress deviator 
apparently does no work, although large plastic de- 
formations as given by the curve BE may be involved, 
The same objection would apply to a deformation ac- 
cording to the tests of Hohenemser and Prager if after 
a course ACE (Fig. 7) the are EH should be followed, 
since here the continuation of the elastic deformation 
vector passes always through B. 

Complete behavior as a Prandtl-Reuss body is also 
unlikely, since in that case, if after the course AC the 
arc CH was followed (Fig. 6), the latter course would 
involve another ratio of the stress deviator components 
(another direction of Io), while no plastic deformations 
would occur. As with constant elastic shearing energy, 
the sum of the squares of the principal shearing stresses 
is also constant, 


oe 


771" const. 


+ 73° = 


> 


a rotation of the stress vector I'y will involve an increase 
of the shearing stresses in one part of the planes anda 
decrease in the other part. It is plausible to assume 
that, in the planes where the shearing stress increases, 
about the same plastic shearing strains will occur as 
when the final state of stress is applied directly, because 
much smaller shearing stresses are sufficient for plastic 
shearing strains in the domain of strain hardening be- 
low the yield stress. This leads to the assumption that, 
for example, with the excess deformation course CD the 
plastic deformation is about BF, so that the direction of 
the total plastic deformation vector A F is between that 
of the original stress vector I) and that of the resulting 
stress vector I'g93. This would mean that only part 
SB of the preceding plastic deformation AB is remem- 
If, at the same time with its rotation, 

i.e., the shearing energy 
of strain hardening, the 


bered (Fig. 6). 
the length of the stress vector 
increases, as in the case 
shearing stress will increase in most planes, so that a 
more complete behavior as a Hencky body could be 
expected. More generally, this will be the case if, 
even at the yield stress, the plastic shearing strain, be- 
longing to the new total deformation, if reached directly 
with constant ratio of the deviator components, is 
larger in most planes than with the original deforma- 
tion—e.g., if the deformation course deviates only 
slightly from the direction of the stress deviator vector. 
This will be possible indeed if one observes that, in 
the same way that Shanley’ showed for bars, the load, 
or here more generally the shortening, increases at the 

same time that the plate buckles. 
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A confirmation of this conception is also given by 
recent tests on tubular specimens of medium steel by 
Fraenkel.'* Representing these tests as shown in 
Fig. 8, where the deformation course is projected on 
the é, — é, plane, it follows from his Fig. 11 that with 
plastic deformation or plastic flow, respectively, the 
stress deviator in the points 1-6 of the total deforma- 
tion course would be represented by the vectors To: 
and Tm, respectively. From Fraenkel’s tests,* how- 
ever, the present author computed the actual vectors 
T), From the vectors in points 3 and 4, of which the 
continuation passes through S; and S,, respectively, it 
follows that the material remembers practically free 
plastic deformations S;B or S,B of 3 to 6 times the 
elastic deformation B-2 

That the results of buckling tests agree with the use 
of preceding plastic deformations (value e) may be 
due to some extent to the fact that the sliding in the 
crystallites, being caused by a migration of Taylor’s 
dislocations, is also going through from one state of 
equilibrium to another, so that with extremely small 
deformations some elastic decrease of the plastic sliding 
in planes with decreasing shearing stress could be pos- 
sible. 

From Fig. 5, where it is consequently assumed that at 
D vector T;, of which the continuation passes through 
A, applies, it follows furthermore that the more the 
direction of the excess deformation course CD deviates 
from that of the initial course AC, the more the direc- 
tion of I’), will deviate from that of Tp) and the larger 
the excess stress deviator vector A,Ip will be. Hence, 
the excess stress and thus the resistance to buckling 
will be the higher the more the proportions of the ex- 
cess deviator components, by buckling, differ from the 
proportions of the original deviator components, before 
buckling.*''** This will be demonstrated by compar- 
ing plastic buckling stresses for several cases of loading 


and boundary conditions (Table 1). 


DERIVATION OF STRESS-STRAIN RELATIONS 


In order to find the buckling stresses in the plastic 
domain, it is assumed that with pure compression an 
arbitrary relation exists between the stress o and the 
plastic strain e, (Fig. 9), so that 


(8) 


With a complete memory of the material the compo- 
nents of the total plastic strain are proportional to the 
components of the resulting stress deviator (in Fig. 5 
the straight line AE is in the continuation of Iq), so 
that by analogy with the elastic domain 


o/eg = E, and da/de, = tan ¢ 


xp = 8,/2Gy, yy = By/2Gy, Exp = 3,/2Gy\ 
(1/2)Yeyp = Tzy/2Gz, etc. 


(9) 


With these relations Eq. (6) yields 


* Mr. Fraenkel kindly informed the author which points of 
his stress and strain curves belong to each other. 
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Fic. 8. Determination of representative vectors of the stress 
deviators for a strain deviator path 0-1-2-3-4-5-6 according to 
tests by Fraenkel (plastic deformation To, plastic flow Te, 
tests Tor). 
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TABLE 1 


Cases of Buckling with Decreasing Difference Between the Ratios 
of the Deviator Components of the Superposed (by Buckling) 
and the Initial State of Strain (Before Buckling). 
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S; = 2G,F, (10) that have been performed since confirm this. With 
. ; pure compression (Fig. 9), it follows from Eqs. (6) and 
in which p= 
(11) that S, = oV2/3 and E, = ¢,V3/2, so that 
ky = 
: - FS S, = (2/3)E = (2/3)E,E 
V exp? + Ge* + ie + (1/2) (vYeyp? + Yue + Yexp’) a (2/ a) i (2/3) iti 


(11) which is in accordance with Eq. (10), as with plastic 
According to tests of Lode,'! Ro&-Eichinger,’* and deformation Poisson’s ratio vy = 1/2 and G, = £,/3. 
others, with constant ratio of the deviator components Furthermore, it follows from Figs. 9 and 10 that tan 
a definite relation S,; = (£,) exists between S, and g, = (2G,/E,) tan g = 2/3tan gy. In order to be able 
FE, (Fig. 10), practically independent of the state of to investigate, besides plates and shells, also three- 
stress. As with a complete memory of the material dimensional cases of plastic stability, the general rela- 
the plastic strains are the same as if the deviator com-__ tions were established between excess strains and excess 
ponents had increased in their resulting ratios, the stresses. Writing the excess strains that occur with 
same relation refers to the present problem. Tests buckling as total differentials gives 


da = = do, + = do, + a * do + = dtm + = dry, + = des | 
di = = da, + = doy + a do, + = dtm + ag dry, + = bres | 
d, = = do, + oe * doy + = do, + rae dtr + = dty, + = Ores | 
12) 
dy = si doz + we de oy + 2 do a+ a ore + dry + = dre | 
dim = ng de, + ig de, + = de. + = dtu + = ite. + — dtu | 
dix = a doz + = do, + ie do, + ie ty + ot tw + = Ores | 


in which the strains may be split in an elastic and plastic part, 


€; = €xe + Erp | (13 


O¢e,/Oo; = (O€/Ooz) + (Vez,/Oe,) 4 
It is evident that 


Oe: 1 O€ye O€z, v OY zye = 1 OV yee _ OY zye ki (14) 
(1 


‘ae ee ae ee a ae == = = () 


do, E od, dy, E Ors G Orn, OTy: 


v being Poisson's ratio in the elastic domain. In order to compute the partial derivatives of the plastic strains 
€.g., O€z/00, = (0/002)(4,/2G,), it is remarked that both ¢, and G, varv. Thus, it is found that 


OG, 2. 06, 1 oS; Gr 15 
Oc, 3 : 


while, according to Fig. 10, 


oO 1 re) 2E; — d(F,/S;) OS; oo Si(dEj/dS;) a Ey dz a 
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With Hence, it follows that 
and Dey _ 2 be _ 1 ds, a, 0 (4) a. (me i. : 
Oo, O0,2G, 2G,00, 2 00, \G, E, 25S,°G, \tan ¢, | 
Oe 0 6 1 06 e © fi 1 5,6, 2G, ery | 
= ae a ae Bp oF ae 7 me a he —— — 1 =— } (17) 
ial 0c, Oc, 2G, 2G,00, 2 00, \G, F,  25S,°G, \tan ¢, de, 
lastic ; 
E./3 OYzyp O Try CsTxy ( 2G, | 
adie = a. aa vas | 
t tan Oo; Oo, Gy SG, \tan 9; ) 
¢ able In the same way it is found that 
three- ' ) 
| rela- OL Bt (= " 1) | 
EXCESS Orr,G, SG, \tan ¢, 
with O€;y = OxTry (= = 1) es OY zp | 
Ory SG, \tan ¢, 0c; | 
> (18) 
Oven _ 1 , Qray? ( 2G, :) | 
Ors G, SiG, \tan ¢, 
‘ 
OV vep _ 2TeyTye (= AS 1) pe: OY xp 
= — = 
OTsy Si°Gp tan Qs Oy: ) 

For convenience the axes of the excess state of stress are placed in the directions of the principal axes of the orig- 
inal one. By this the original shearing stresses r,,, ty, and 72, disappear in Eqs. (17) and (18), with the result 
that all partial derivatives of shearing strains y or with respect to shearing stresses r become zero, except for the 
case in which 7 and y have equal subscripts. 

12 In this way it follows from Eqs. (12), after using Eqs. (13), (14), (17), and (18), with E, = E/e and denoting 
differentials such as do, by a,’, that* 
l+e 36,’ EK — etan 2v+e 36,6, E — e tan } 
6 = ( —+— ; 4 op’ + (- a = ew —* o + 
E 2S;, Etang 2E 28,7 Etang 
(as, Se el | 
2E 2S Etang 7 
2v>+e 36,6, E — e tan l+e  36,?E — etan 
yn (- Met Bae astne) 4 (Lee, ome), | 
2E 2S, Etang E 2S;7 Etang 
2+e  36,¢,E —etang <a 
- one + 6st 
2E 2S; Etang { 
(13 { 9 ie « . 2s = . 
/ 2Qv+e 36,6, E — etang 2v+e 36,6, E — etang | 
= (- > ae ae ) oe’ + (- + See ee a,’ + | 
2E 2S;2 Etang 2E 2S Etang 
(- +e 36, E- < ten #) ; 
aes ~~ i A de . 
E 25> Etan ¢ , 
‘ , 2+ 2» + 3e F ? 2+ 2v + 3e ‘ y 2+ 2+ de , | 
t=. = E Tzy, Yur = E ines Ta = E ” Taz } 
aims : . 
Inversely expressing the differentials of the stresses in terms of the strains, expressions of the following form were 
obtained: 
15 oO,’ - Cree’ + Cire,’ + Cis€," | 
oy = Crier’ + Carey’ + Case,’ 
j a,’ = Csr’ + Case,’ + Case,’ \ (20) 
: Tay" = Curry | , 
| Tus" = Cissus’ 
Tex’ = Cure’ | 
in which the constants C are rather extensive and will not be written down here. However, 
| (16 * Consequently, value e = E/E, = (7/ee)/(o/€p) = €p/€e, in which ep and ¢, are the plastic and elastic strains, respectively, for the 





equivalent pure compressive stress ¢ = Vo, — ozoy + o,* at which the plate buckles. 
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Cu = Cos = Coe = E/(2 + Qo + 3e) 


Thus, apparently the deformation is governed by seven different constants. This anistropic behavior is be- 
tween the elastic behavior of a crystal of the rhombic and one of the regular system. In the regular system the 
number of constants is diminished from nine to three with respect to the rhombic system, because of the equival- 
ence of the three axes. In the present case, however, this equivalence exists only for the shearing stresses, so that 
only Cu = Css = Cos, by which 9 — 2 = 7 constants remain. 

In order to study three-dimensional cases of plastic stability, relations (20) may be inserted in the differential 
equations of the general theory of elastic stability of Biezeno and Hencky.*°—*? 

With plane states of stress, as occur in plates and shells, in Eq. (19) o,’ is zero, while ¢,’ does not interest us. 
Expressing ¢,, ¢, and S,* in terms of the principal stresses ¢, = p; and o, = 2, while o, is zero, it is seen that 

1 


< Diss - 4. 
a2 = 5 (2p: — p2), G, = 3 (2p2 — pi), Si? = 3 (P1” — Pipe + p2”) 


Expressing again the excess stresses in terms of the excess strains, it follows now from Eqs. (19) that 


on E(Ae,’ + Be,’), o,’ = E(Be,’ + De,’) 


inwhich A = ¢:/¢14, B= g2/¢i, D = ¢3/¢4 
and g = (1 — 28)?E + (4n? + 3e) tan ¢ 
go = (2 — B) (1 — 2B)E + (Arn? + 3e8) tan ¢ 
¢g; = (2 — B)?E + (4n? + 3e6?) tan ¢ > (2) 
g = [(5 — 4v) (1 + B*) + 2(5r — 4)6 + B8en?JE + 
}4(1 — v?)n? + 3e[n? + (1 — 2v)e8l|{tan ¢ 
Bp = p/p: and rv =—2-—-B+1 


while tay’ = [E/(2j+ 2v + 3e) lyzy’ = EF yxy’ 


In reference 13 the author derived these same equations in another way—-namely, by subtracting the initial 
stresses and strains before buckling from the total stresses and strains after buckling. 
In Eqs. (21) values A, B, and D may also be expressed in terms of the tangent modulus E, = do/de instead of in 


terms of 
tan ¢ = da/de, = da/(de — de.) = 1/[(1/E,) — (1/E)] = EE,/(E — E,) 


Introducing this expression for tan gin Eqs. (21) it follows that 


A= Vi/W, B= v2/V, D= v3 Ws 
in which yy (1 — 28)?E + 3(1 + ez, 
¥2 = (2 — B) (1 — 2B)E — [(2 — 8B) (1 — 28) — 4vn? — 3eB]E, 
v3 (2 — 8)*E + 3(1 + e)BE, 
We = [(5 — 4y) (1 + 6?) + 2(5y — 4)B + Ben? JE — 
(1 — 2») [(1 — 2») n?-3(1 +e) 8] £, | 


ll 


21a) 


while e = (E/E,) — 1, in which £, is the secant modulus. 
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DERIVATION OF THE DIFFERENTIAL EQUATION FOR tially elastic behavior has only academic value and 
PLASTIC BUCKLING OF PLATES does not influence at all the real buckling stresses with 

the existing eccentricities. It is, therefore, logical to 
' 


buckling, the stress distribution in the buckling plate take for the critical stress of columns o, = r°E,/X*, in 
has to be examined. As was shown by the author,!? which E, = do/de is the tangent modulus. This is in 
when buckling occurs at the yield stress, the entire agreement with general practice in aeronautical engi- 
plate deforms plastically. Furthermore, Chwalla neering. Shanley'* has shown that even a centrically 
showed that with rectangular bars a strain reversal in loaded column cannot remain straight above this 
the plastic domain according to the Engesser-Karman critical stress, while according to tests the maximum 
theory will only occur if the eccentricity is less than 1/50 stress does practically not surpass o;. Consequently, 

of the thickness. In practical cases one has totakeinto there is still more reason to assume that for plates and | 

account larger eccentricities, so that from the beginning shells the buckling stress is equal to that which would | 
the author contended that it was logical to base the al- occur if the entire plate deformed plastically. With 

lowable buckling stresses for bars as well as for plates the infinitesimal deflections assumed for buckling this | 

on a completely plastic behavior.'***°> The par- leads to a linear stress distribution, so that the bending | 

; 

: 

' 


In order to derive the differential equation for plate 
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PLASTIC STABILITY 





and twisting moments can be expressed in terms of the 
excess stresses. If, for example, a,’ is the stress at 
the boundary, the bending moment is M,’ = (h?/6)o;,’. 
After using Eqs. (21) and expressing the strains in 
terms of the deflections w, so that e,’ at the boundary is, 
for example, — (1/2)hO*w, Ox’, this leads to the bending 
and twisting moments 





i -» sabaoe +e =) 
alle _ - (. Ox oy? | 
O*w =) 
Mt se = RT TB D 22) 
’ ( Ox? + ra 9 4 ( j 
O*w | 
M,,'’ = —M,,! = 2EIF — 
, = dxdy } 


in which J is the moment of inertia 4°/12 of the plate. 


After insertion of these values in the well-known con- 
dition of equilibrium with reference to a translation in 
the z-direction, 





eM.’ OM,’ CM,’ } 0*w 44 *w 
—2—7% + -=h 1p, — 
Ox? Oxo" oy? Pt 2 - Oy? 


the following governing differential equation for the 
plastic buckling of plates is obtained :"* 


EI { O*u >/R PT Ow 4 D O*w 
a | A — + AB + 2F) —— — 
Ox! Oxy" oy! 
Ow O*w 
hp, +hoe— =0 (23) 
Ox* oy? 
In a similar way the differential equations for the 
buckling of shells may be derived.”**° 
BUCKLING STRESSES OF PLATES AND SHELLS WITH 


VaRIOoUS LOADING AND BOUNDARY CONDITIONS 


With these differential equations the buckling stresses 


for various loading and boundary conditions have been 


determined. In order to show that indeed the 
plastic buckling stress cx, is higher as compared to that 
in the elastic domain, c,, the more the ratio of the devi- 
as superposed by buckling, differs 
original deviator components. In 


ator components, 
that 
Table 1 the ratio cy, ¢; 


from of the 
is given for twelve cases for mild 
steel buckling just below the yield stress.*1;*?*4 

With Case | for example, an infinitely long plate 
subjected to compression in axial direction, of which 
one unloaded side is simply supported while the other 
one is free, with buckling pure shearing stresses are 
superimposed on the original pure compressive stresses, 
yielding a high ratio og/cg. With Case 2, besides the 
Shearing stresses of the twisting moments, bending 
stresses co,’ and c,’ are superposed on the original com- 
pressive stresses co,, sO that og/og¢ is somewhat less. 
For Cases 3 and 4, again, the twisting moments are less 
pronounced. With Cases 5 and 6, bending stresses 
o;' and og,’ are superposed on the original shearing 
stresses, but with Case 6 the shearing stresses due to 
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the twisting moments are more pronounced. There- 
fore, for Case 6, the ratio r,/7, is smaller than for Case 
5, because in both of these cases the resistance to shear- 
ing stresses is the smaller one. With Case 10, for ex- 
ample, with buckling biaxial bending stresses are super- 
posed on the original ones, yielding a small ratio og/ag. 
With buckling of a bar, Case 12, exactly the same state 
of stress is superposed on the original one, so that the 
ratio o/c, is as small as possible.* 

The theory proved to be in splendid accordance with 
the beautiful and systematic tests on plates of avional, 
an aluminum alloy, as published by Kollbrunner in 
Ziirich 8 years later.* The number of half waves in 
which the plates buckled proves that the plates behaved 
anisotropically to the same extent as predicted by the 
author's theory. 

In Tables 2, 3, and 4 the data for the Cases 2, 3, and 
4, as given in Table 1, with lengths a and widths 6, 
have been compiled. According to the theory the 
number of half waves is pg, while a, — a, represents 
the limiting lengths between which this number of half 
waves will occur. The number of half waves as found 
by Kollbrunner from his tests is given by p,, followed 

* When the author wanted to apply his theory to the buckling 
of the webs of built-up sections of compression members, he came 
to the conclusion that the existing theories of Bleich™® and 
Chwalla” on this matter were incorrect in the elastic domain, 
since they did not take into account in the right way the re- 
straining action of vertical and horizontal plates upon each other. 
Therefore he developed an exact theory of the buckling of those 
sections by writing down the general solution of the differential 
equation for the deflection surfaces of vertical, as well as hori- 
zontal, plates after which the conditions for the connection of all 
plates yield the buckling condition.2**-** This theory alse 
takes account of the fact that, in general, with buckling of the 
webs the intersection lines of the middle surfaces of the plates de 
not remain straight. For soft steel design, formulas were de- 
rived for the ratios b/h of the webs as a function of the slenderness 
ratios //r of the members, which were chosen in such a way that 
the buckling stress of the webs is equal to that of the members 


29 


as a whole.**? 























TABLE 2 
a j 3.20 4850 0450 8.16 
a, —a, | 2496-3400 | 4.316— 5.216 | 6.100 — 6.98b | 7876—8755 
Pa 3 5 7 9 
Pr 3(3) 33 7 (3) 9@) 
TABLE 3 
| a I8D 5.056 ) 9456 
| a, —a, | 330b—4.056 | 544b—6.126 | 7.506—8.206 | 888b-- 9.566 
| Pe 5 8 1 13 
| py 5 (1), 6 (3) 8 (3) 10 (1), 11(2) 13(1), 14(2) 
TABLE 4 
a 4.556 6.80 b 9106 11.35 
| a, — ay | 4.116 — 4.586 | 6436—683b,, 8.736—9.196 |11.04> —11.506 
Pp Q } 14 } 19 24 
Pr 9 (3) 13(2), 14) 18 (2), 19(1) | 24 (3) J 
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Comparison between theory and tests by Kollbruiner 
and Stiissi 


Fic. 11. 


by the number of tests in brackets. The conformity 
between theory and tests is better than could be de- 


manded. 


The stress-strain curve of the avional used is given 
by the full curve in Fig. 1]. In the same figure the 
buckling stress og, is plotted as a function of oy £, as 
is customary in this country. The tangent modulus 
curve gives the buckling stress for bars (Case 12, Fig. 
13), while the stress strain curve would apply for plates 
if the ratio og/o, were equal to the ratio /,/ /, in which 
k, is the secant modulus, as proposed by Gerard.** 


For Case 2 (Table 1) the plastic buckling stress is!*°4 


op = (29°EI/hb?) (VW AD + B + 2F) 


while in the elastic domain, with v = 0.3, the coeffi- 
cient of the last term is equal to 2.2, so that the ratio 
op/o, givenin Table 5 is obtained. 


The exact buckling stress for Case 4 follows from a 


transcendental equation'’**?° but may be written 


sufficiently accurately, as given in Table 5.*! 

In Fig. 11 the plate compressive curves (buckling 
stresses) according to these formulas are given for 
Cases 2 and 4. The curve for Case 3, not shown, falls 
between the curves for Cases 2 and 4. The average 
values of Kollbrunner’s tests are given by the open 
circles, denoted by the case in question. Tables 2, 3, 
and 4 refer to the tests indicated by the highest circle 
2 and circles 3 and 4, respectively. The discrepancies 
in the bucking stresses are apparently due to unavoid- 
able eccentricities. 
theoretical and experimental half waves may be ex- 
plained by the fact that these eccentricities affect the 


The exact conformity between the 


rigidity of the plates in both directions to practically 
the same extent. 
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New tests, executed in Ziirich with an improved ap- 
paratus, are given in Fig. 11 by the three full circles, 
and these give still better agreement.*! 

If the material had no memory at all for preceding 
plastic deformations (plastic flow), value e would have 
to be equated to zero in Eq. (21).'4 
tire plate to be plastic, this would lead, for Case 2, to 
According to the 


Assuming the en- 


the curve denoted by 2, e = 0, pl. 
tests this curve would apparently overestimate the 
buckling stress considerably. Still higher buckling 
stresses would result if at the same time elastic regions 
according to the Engesser-Karman theory were as- 
sumed (curve 2, e = 0, el. + pl.). 

A still better agreement with the author's theory is 
shown by the results of tests with Z- and H- sections 
and channels in extruded aluminum alloys, as published 
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PLASTIC 


by the N.A.C.A.“!~-** Using the method given by 
Schuette and McDonald,® the plate compressive curve 
that would follow from these tests for avional (with a 
stress strain curve according to the full curves in 
Figs. 11, 12, and 13) is given in Figs. 12 and 13 by the 
full curve connecting the eight circles. If the width 
to thickness ratio of the flanges is more than about one- 
third of that of the web, the flanges will be elastically 
built in by the web, so that the circles should be situated 
between the theoretical curves for Cases 1 and 13, 
Table 5, with one unloaded side free and the other 
simply supported and clamped, respectively. 

For Case 1 and sufficiently long buckling length, the 
buckling stress is °*** 


og = EFh?/b? 


while in the elastic domain F = 1/2(1 + v) = 1/2.6, so 
that the ratio og/og is as given in Table 5. The cor- 
responding plate compressive curve in Fig. 12 is de- 
noted by 1. If Poisson’s ratio for elastic deformations 
were equated to 0.5, according to Eq. (21) value F 
would be equal to 1/3(1 + e), while in the elastic 
domain, with e = 0, it would be '/3, which would yield 


a ratio 
Cp/CR = ¥ (1 + e= E,/E 


since in the author’s notation the secant modulus &, is 
equal to #/(1 + e). In this case the plate compressive 
curve for Case 1 would coincide with the stress-strain 
diagram. 

For Case 13 (Table 5) the buckling stress may be 
calculated approximately according to the energy 
method in order to find it in explicit form. As only 
the ratio between plastic and elastic buckling stress 
is needed, the accuracy of this ratio will still be suffi- 
cient. It follows that the ratio will be that given in 
Table 5. It may be stated that here, in contradiction 
to other cases of plate buckling, a higher value of B, 
representing the lateral contraction, involves a lower 
This is a consequence of the fact that in 
this case the plate buckles according to an anticlastic 


value of op. 


surface, so that the higher the ratio of lateral contrac 
tion, the more the curvatures in longitudinal and lateral 
The plate com 
pressive curve for this case is denoted in Fig. 12 by Case 


directions will promote each other. 


) 
su. 


The experimental N.A.C.A. curve is between the 
theoretical curves for Cases 1 and 13, indicating agree- 
ment between theory and tests. 

Assuming plastic flow (e = 0), the theoretical curve 
for Case 1 would be a straight line (line 1, e = 0 in 
Fig. 12), predicting apparently much too high buckling 
stresses. 

In case the width to thickness ratio of the flanges is 
less than one-third of that of the web, the latter will be 
less buckling proof and will behave as elastically built 


in by the flanges. The theoretical curve would be 
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TABLE 5 





BOUNDARY CONDITIONS|] A= 
AND LOADING 2/F, 


|—ET «| 


CASE RATIO o,/o, OR T/T, 





2 J CE [ ft) 0.455(,/A0+8+2F) 
4 qtC<*éE ST ry 0.598/aD+0.312(8+2F) 
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situated between the curves for Cases 2 and 4 in Fig. 11, 
so that here also good agreement between theory and 
tests exists. 

It follows, however, from Table 1 that for other cases 
the ratios ¢,/o, Of Tg/Te May be much lower. The 
general formulas for these ratios are given in Table 5 
for the Cases 6, 9, and 12 of Table 1,'*:**' for a cylindri- 
cal shell under uniform axial compression®**> (Case 
14), and for a uniformly compressed spherical shell 
(Case 15), both for buckling symmetrically with respect 
to the axis or a diameter, respectively. The correspond- 
ing curves for Cases 6, 9, 12, and 15 are given in Fig. 13, 
the lowest curve, that for a bar (Case 12), being deter- 
mined by the tangent modulus.* 

It may be observed that, in order to compute the 
ratios for Cases 14 and 15 as given in Table 5, the 
buckling stress was assumed to be a continuous func- 
tion of the wave length, while in the latter case smaller 
terms were neglected, both as usual in the elastic do- 
main. Therefore, these two ratios do not apply ex- 
actly for the extreme case of buckling at the yield 
stress, where tan ¢ is zero, for which both ratios appear 
This may be under- 
stood, since, for example, for a cylinder the correspond- 
ing most unfavorable buckling length (half wave length) 
is infinite, so that the strips parallel to the axis would 


to become zero just as for a bar. 


* For Cases 6 and 15 in Fig. 13, values of og and og refer to the 


y son 
equivalent pure compressive stresses ¢ = V o,?—ozoy+oy?+3rgy3 
under which the plate buckles. 
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have no proper buckling strength. If, furthermore, 
the shell buckles symmetrically with respect to the 
axis, the necessary strain in the circumferential direc- 
tion may be effected by a shortening in the axial direc- 
tion, which requires no excess stresses at the yield 
stress, so that with this mode of buckling no excess inner 
stresses are involved, as with buckling of a bar. In 
reality, the half wave length //p cannot be more than 
the length of the shell, denoted by /, so that the buck- 
ling force per unit breadth is***® 


(EIAr*/I2) + [D — (B2/A)](Eh/a?)2/m? 


N; = hog = 
At the yield stress D — (B*/A) = 0, so that 
N, = hog = EIA /I? 


in which J = h*/12, h is the thickness, and 2a is the 
diameter of the shell. At the yield stress the buckling 
stress is the same as that for a wide plate (Case 9 in 
Tables 1 and 5) with the same buckling length /. 


Case 9 also applies for a tube under external pressure, 
so that the theory also solves, for example, the plastic 
collapse of deep well casing under external pressure, 
which was dealt with before by using the same ratio 


¢,/o, as applies for bars.** 


Later Papers on the Same Subject 


In his first papers the author, in 1938, found that 
(assuming the plate to remain straight up to the buck- 
ling stress, as in the theory of Engesser-Karman for 
bars), even for a wide plate that buckles cylindrically 
(Case 9), the boundary surface between plastic and 
elastic region is curved.'*** This results im a varying 
flexural rigidity so that the plate will not buckle ac 
cording to a sine wave. As this makes an exact solu- 
tion practically impossible and as the author was of the 
opinion that the elastic regions do not influence the 
actual buckling stresses in practical cases, he did not 
go further in this matter and reckoned with entirely 
plastic plates.** Ilyushin (1944)** refers to the pres- 
ent author's efforts on the determination of this bound- 
ary surface and intended to give an exact solution for 
it. In deriving the differential equation of the plate, 
however, he assumes that the resulting normal and 
shearing forces in buckling do not vary, yielding a 
boundary surface parallel to the middle plane. It fol- 
lows, however, from the present author's earlier papers, 
that both of these conclusions are incorrect.!**° 
Ilyushin’s differential stress-strain relations for the 
plastic region are the same as the author's Eqs. (21) or 
(21a) if Poisson’s ratio v for the elastic region is equated 
to 0.5. Using this assumption gives, in general, too 
high buckling stresses. 

Handelman and Prager*® also conclude that the 
above-mentioned boundary surface is parallel to the 
middle plane of the plate. Furthermore, they con- 
sider only buckling under pure compression and plastic 
flow, by which their stress strain relations follow from 
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Eqs. (21) or (21a) by equating 8 = p2/p, and e to zero, 
This assumption is dealt with here (curve 2, e = 0, el. + 
pl. in Fig. 11 and line 1, e = 0 in Fig. 12) and proves 
to give too high buckling stresses. 


The stress strain relations of Stowell (1948)*" are the 
same as those of Ilyushin, and, where he also considers 
buckling under pure compression, they follow di- 
rectly from Eqs. (21) or (21a) by equating 3 = pe, p; to 
zero and vy to 0.5. He considers the entire plate plastic 
in his buckling analysis and thereby obtains results in 
agreement with tests. 
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Interference Between Wing and Body at 
Supersonic Speeds—Note on Wind-Tunnel 
Results and Addendum to Calculations’ 


CARLO FERRARI* 
Cornell Aeronautical Laboratory, Inc. 


ABSTRACT 


A correction is made to the numerical example worked out ina 
previous paper for the pressure distribution on a rectangular wing 
operating in the upwash from a pointed inclined body of revolu- 
tion in a supersonic stream. Wind-tunnel results confirm the 
amended calculated pressures and bring to light the drastic falling 


off of lift in the wing-body juncture 


INTRODUCTION 


— NOTE IS PROMPTED through realization that 
part of the numerical example used to illus- 
trate the theoretical results of a previous paper on 
wing-body interference by the author! contains an 
error in one of the boundary conditions. In order to 
assure a safe basis for proper application of the boundary 
conditions, it is necessary to review this facet of the cal- 
culations in detail. The boundary condition, 
which is applied incorrectly in reference 1, is depicted 
This oversight resulted from 


some 


in Fig. 10 of that paper. 
an attempt to bring together several steps in an itera- 
tion process, but it appears that a more carefully worked- 
out sequence should be followed. The basis of this 
iteration process lies in the following considerations. 


REVIEW OF THEORETICAL STEPS 


General Observations 

The solution of the wing-body interference problem, 
by the method of small perturbations that is carried out 
in the theoretical part of reference 1, consists essentially 
of finding a potential (or the velocity field derived from 
such a potential) for the wing-body system as a result of 
adding to the potential for the undisturbed impinging 

-” 

stream, I’,, the sum of the potentials related to the 
appropriate distributions of supersonic doublets which 
exist on the surface of the wing and on the body. These 
distributions are uniquely determined by the boundary 
conditions, which will be examined more critically be 
low. 

In conformity with the usual procedure followed in 
the linear theory, the distribution of supersonic doublets 


Received March 25, 1949. 

* This work, both experimental and theoretical, was carried out 
under Navy BuOrd sponsorship. 

t Member of C.A.L. consulting staff. 
di Aeronautica, Politecnico di Torino 


Professor, Laboratorio 


is obtained along the body axis and not at the surface, 
but this is entirely legitimate. Likewise, the distribu 
tion for the wing actually occurs on the x,y-plane in the 
region bounded by the wing surface, including the part 
extending inside the body, when the wing is considered 
as flattened down upon the x,y-plane. Here again the 
procedure conforms with the usual one used in linearized 
theory. 

Within the framework of the linearized theory, an 
exact solution of the wing plus body problem may thus 
be found by a process that depends upon the following 


precise formulations: 


Stepwise Detailed Procedure 


Let >, be the distribution of supersonic doublets in 
the region of the whole wing, including the extension 
within the body, and let 2. be the distribution along the 
body axis. Now represent by Vy1, Va2, V,. the com- 
ponents of velocity normal to the wing (flattened), 


which are due to the distributions 4), »., and free 
— 
stream, V’., respectively. In addition, let Vy’, Vn’, 


l’,,.’ denote the velocity components that are normal to 


the body surface, where these components also arise 


from the distributions ~,, D2, and V., respectively. 
Consequently, the supersonic doublet distributions 
will be determined through satisfaction of the equations 


V nt + Vi. + Vio =(Q and Var’ + Vago’ 4-Vio- 0 


(on the wing) (on the body (1 


When these expressions are written out in explicit 
form, they constitute a system of integral equations of 
the Volterra type, and they may be solved by the 
method of iteration, as follows: 

Step 1.—Calculate 2; by means of the boundary con- 
dition Vai + Vaew = 
(this is the zero approximation, equivalent to neglecting 


0, where V,2 has been set = 0 
all interference effect). 

Step 2.—Calculate 2 by means of V,2" + V, 
where V,;’ has been set = 0 (this is the zero approxima- 


- (), 


tion for the body distribution). 

Step 3.—Recalculate 2, using the first relation of Eq. 
(1), wherein the value of V,2 is that obtained from the 
S. calculated in Step 2 (this is now the first-stage ap- 


proximation for the wing). 
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Fic. 1. Upwash distribution out along wing span. 
Step 4.—Recalculate ~, using the second relation ot 


Eq. (1), wherein the value of V,,' is that obtained from 
the X; calculated in Step 1 (this is now the first-stage 
approximation for the body). 

Step 5.—This process is continued by employing, in 
the second-stage approximation, the distributions >, 


? 


and XY», which were obtained in Steps 3 and 4. 

Jecause the integral equations are of the Volterra 
type, it is assured that this iteration process will con- 
verge to give the same solution as would be obtained if 
both sets of boundary conditions were satisfied simul- 


taneously. 


Recasting of Step 3 as an Incremental Distribution 


In applying this procedure to the problem at hand, it 
is essential to note that Step 1 is the two-dimensional 
solution for a flat plate inclined at an angle 8:—e., 
there is a uniform down-wash of magnitude V.62 in 
the region confined between the two Mach planes 
emanating from the leading and trailing edges of the 
wing, Step 2 is the Tsien or Ferrari (see references 
given in original paper) derivation for the flow around 
an isolated body at an angle of attack, §;, to the free- 


stream flow. For the purposes of Step 3, it is found 
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more convenient to solve for the increment in the dis- 
tribution of doublets on the wing (flattened) rather 
than for the first-stage approximation to , directly. 
Consequently, the upwash increments on the wing due 
to the flow around the inclined body are just as calcu- 
lated previously, but in the region of the wing, which is 
covered by the body, the increment upwash distribution 
will, of necessity, be equal in magnitude but of opposite 
sign to the basic distribution for the wing, determined 
as Step 1—1.e., the total solution for this Step 3 requires 
that there be no flow through the body, where the 
boundary condition is again being applied in the x,y- 


plane. 


Numerical Calculations of Upwash Increments 


The correct boundary condition for finding this first 
stage approximation to the increment of lift on the wing 
due to interference is thus the one depicted here by 
means of a solid line (two pieces) in Fig. 1. This span- 
wise distribution of upwash is taken to apply inde- 
pendently of the fore-and-aft location on the wing with 
good approximation, because it is the average of two 
distributions, one calculated at a section near the lead- 
ing edge and the other near the trailing edge. That 
these two calculated spanwise distributions at opposite 
edges of the wing are almost identical, to the degree of 
accuracy expected from this linearized theory, is borne 
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out by the curves plotted in Fig. 2. The values of the 
upwash at the three spanwise stations, » = 0.25, 0.5, 
and 0.75, which were obtained previously (see Fig. 10 of 
reference 1) through assumption that the body diameter 
was infinitesimally small, are also reproduced on Fig. 2 
in order to re-emphasize that the upwash varies but a 
slight amount from front to rear locations in the region 
in which the wing is situated. In addition, these check 
points also agree so well with the more detailed step- 
wise calculations represented by the two smooth curves 
that they give good assurance that the former procedure 
of assuming the body diameter to be infinitesimally 
small [see Eq. (44) of reference 1] was justified for all 
practical purposes. 

For orientation it should be remembered that the 
configuration selected for computation is one having the 
general features as represented in Fig. 3, and the specific 
case of Mach Number = 2.0 and angle of attack = 
+8° is the one chosen for numerical illustration. Of 
course, under these conditions, the new correct bound- 
ary condition applicable to the incremental upwash in 
the region of the wing that is covered by the body 
necessitates a neutralizing value of upwash equal to 
— 0.1396 (—S° ), so that these additional normal velocity 
components actually constitute down-wash increments 


here. 
Trigonometric Representation of Upwash Distribution 


By means of a curve-fitting procedure that is entirely 
analogous to that already illustrated in reference 1, it 
will be found that a trigonometric series of cosine terms 
of odd multiples of the argument will represent the now 
amended boundary condition in a manner shown by 
means of the dashed line in Fig. 1. This trigonometric 
representation is required in order to carry out subse- 
quent steps in the theory. 


The formula of Clairaut-de la Vallée Poussin used in 
reference 1 gives a trigonometric series representation 
that passes exactly through a discrete number of selected 
points. The same series may also be obtained by using a 
tabular scheme such as devised by Runge and illustrated 
by J. P. Den Hartog in reference 2. The dashed curve of 
Fig. 1 is obtained by use of Runge’s 4S-Point Tabular 
Schedule, wherein the calculated upwash values for » = 
0 to 1.0 are extended antisymetrically from » = 1.0 to 
2.0, and the resulting distribution is then extended again 


to give an even function in the interval 7 = —2 to 7 
+2. This process thus defines an upwash distribution 
whose trigonometric representation is composed solely 
of cosine terms with odd multiples of the argument and 
the period of which is 4. The subsequent determina- 
tions of pressure distributions out along the wing span 
will also be expressed in terms of such periodic functions, 
but the detailed behavior of these distributions is only 
of interest in the quarter interval from 7 = 0 to 1, which 
corresponds to the actual wing semispan—that is, the 


pressure distributions are of the form 
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Fic. 3. Configuration used in the numerical application and 


upon which the wind-tunnel tests were performed. 


C, = 24> H,” cos [(1/2)rmn] 
m 
where the superscript (7) denotes that this particular 
spanwise distribution applies over the small chordwise 
strip, wherein J/,,"". = dh,»(&) dé is assumed constant, 
For purposes of the numerical application made here, 
the chord has been divided into only three sections of 
equal length, and the resulting distributions have been 
designated simply as applying to locations at one-sixth 
chord aft of the leading edge, at the wing axis (50 per 
cent chord points), and at one-sixth chord forward of the 
trailing edge. The notation follows exactly that em- 
ployed previously, and in particular » = 0 corresponds 


to the root chord (plane of symmetry), while » = 1.0 
represents the wing tip. 


Fourier Representation and Occurrence of the Gibbs 

Phenomenon 

Close scrutiny of Fig. | discloses that there are rather 
large discrepancies between the original upwash dis- 
tribution and the trigonometric representation in the 
region of the discontinuity at the body juncture, even 
though the 4S-Point Schedule allowed a determination 
of twelve terms (m7 = 23). It may be more satisfactory 
to use a Fourier representation and to take a great num- 
ber of harmonics in order to match closely the calculated 
upwash distribution, but such refinements are not 
needed to show the acceptable correspondence between 
theory and experiment, which is presented in Fig. 5. 

It may be noted that the behavior of such a Fourier 
representation in the neighborhood of the discontinuity 
is characterized by the Gibbs phenomenon, which pre- 
dicts that large discrepancies will be maintained between 
the given curve and the representing Fourier develop- 
ment regardless of the number of terms employed 
This difficulty can be obviated by an appeal to the 
smoothing effect of the boundary layer on the mathe- 
matical sharpness of the discontinuity—i.e., it 1s 
obvious physically that the discontinuity at the wing- 
body juncture is not actually as precise as depicted by 
the solid line in Fig. 1. Consequently, the discontinuity 
can be replaced legitimately by a continuous smooth 
curve that will deviate so slightly from the one having 4 
discrete jump that for all practical purposes the smooth 
curve can be employed without fear of any essential 
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WING-BODY 


discrepancy arising; in fact, a closer approximation to 
physical truth is likely to follow. Of course, such a con- 
tinuous curve can then be represented at every point to 
any desired degree of accuracy by a Fourier series con- 


sisting of a sufficient number of terms. 


Computational Results for the Pressures on the Wing 


Once the trigonometric representation of the upwash 
has been obtained to the desired degree of accuracy, the 
subsequent calculations required to determine the incre- 
ments in pressure distributions on a wing placed in 
this field of upwash may be carried out through 
following a procedure already illustrated in reference 1. 
The results of such a calculation in the present case turn 
out to be the three distributions plotted in Fig. 4, which 
apply to the three fore-and-aft sections into which the 
wing was arbitrarily divided. It is significant that a 
wider and wider region near the wing-body juncture 
shows a falling off of the interference-induced pressures 
as the trailing edge is approached, in distinction to what 
occurs at the leading edge; of course, the foremost tip 
of the leading edge supports the two-dimensional distri- 
bution, but no two-dimensional theory would account 
for the interesting phenomenon occurring farther aft 
in the wing-body juncture. (See Fig. 4.) 


EXPERIMENTAL RESULTS 
Description of Model and Stream Conditions 


Wind-tunnel tests have been made on a configuration 
possessing the same geometric characteristics as used in 
computing the above-mentioned numerical results. 
The wing was a 10 per cent thick double wedge with a 
flat upper surface that was aligned with the centerline 
of the body. The orifices were located on the flat sur- 
face, and thus 
negative ranges of angle of attack gave the equivalent 
The 
experiments were carried out at the Naval Ordnance 
Aerophysics Laboratory at The 
test section width is 19 in., while the height is 27.5 in.; 


complete run of both positive and 
of the flat-plate condition analyzed numerically. 
Daingerfield, Tex. 


thus, the 3.0-in. diameter (maximum along cylindrical 
section) model is accommodated with sufficient clear- 
ince to ensure the absence of shock reflection between 
the model and the tunnel walls at all angles of attack 
tested. The flow was close to 1J = 2.0, and calibration 
surveys reveal that this working section is remarkably 
free of anomalies. The Reynolds Number of the free 
stream, based on the 2.1-in. wing chord, was approxi- 


mately 1,270,000 (wing span from tip to tip = 12.0 in. 
and overall body length = 26.25 in.). The moisture 


content was controlled to ensure a dew point in the 
to —4°F. 


range +19 
Comparison with Theory 


The experimentally determined pressure coefficients 
along a spanwise row of ordinates located on the wing 
ixis (50 per cent chord point) are plotted in Fig. 5 for 
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+ §°, 4°, and 0° angles of attack, and the theoretical 
first approximation is presented for confirmation.* 


Build-Up of Pressure Coefficient from Basic Components 


The theoretical curve is obtained by use of the incre- 
ment pressure coefficient for the wing axis given in Fig. 
, to which is added the basic coefficients represented in 
Fig. 6. The basic distribution for any angle of attack 
is obtained by computing the two-dimensional pres- 
sures for the corresponding angle of attack (which 
change sign as the angle of attack changes sign) and 
adding the constant ‘“‘contour’’ effect shown in dotted 
line on Fig. 6. The ‘‘contour’’ effect is merely the ra- 
dial distribution of pressure at the particular longitudi- 
nal station along the missile axis at which the wing axis 
happens to fall; this radial distribution of pressure for 
zero angle of attack results from the longitudinal com- 
ponent of the impinging stream and accounts for the 
distortion to this stream which occurs as the flow makes 
accommodation for the bulk of the body. The distor 
tion is appreciable in the present example because the 
nose, with its rapidly changing contour, is so near the 
wing. The theoretical variations in the pressures en- 
countered at the wing tip are also included in the basic 
distributions. This use of theoretical basic pressures 
will not introduce any error in the final comparison of 
theoretical and experimental pressure coefficients, 
which include interference effects, because exact corre 
spondence was found to exist (within the degree of accu- 
racy of the experimental results) between wing-alone 
tests (not presented here) and the theoretical basic 
values just mentioned, which pertain to the wing before 
interference increments are added. Examination of 
Fig. 5 also discloses that the magnitude of the 
tour’’ effect from the body, for zero degree angle of at- 
tack, is also adequately confirmed by the experimental 
data. 

The correspondence between experimental points and 
the theoretical result obtained by summing the angle of 
attack and “contour” effects of the basic distribution 


with the increment due to interference occurring along 


“con- 


the wing axis, as shown in Fig. 4, is extremely close, even 
for this first approximation step. The adequacy of the 
small perturbation theory to deal with this phase of the 
wing-body interference problem is manifest, particu- 


*The author wishes to express his indebtedness to R. H 
Cramer, of the Cornell Aeronautical Laboratory, for planning and 
performing the experimentation, as well as the analysis which pro 


vides this comparison of theory with test data 
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distribution components for build-up of span 
wise pressure variations. 


Fic. 6. ‘“Basic’’ 


larly because the experimental data of Fig. 5 confirm 
falling off of pressure in the wing-body juncture. 


Influence of Wing upon the Body 


The calculations previously given for determining the 


influence of the wing upon the body need no revision, 


but it should be recognized that the numerical data 
adduced are for the first approximation step only 
More detailed comparisons of wing-body interference 


for the configuration tested are expected to be made 
available for publication at a later date in a separate 
paper to be devoted entirely to presentation of such ex- 


perimental results. 
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Temperature and Velocity Profiles in the 
Compressible Laminar Boundary Layer 
with Arbitrary Distribution of 


Surface Temperature 


DEAN R. CHAPMAN* ann MORRIS W. RUBESIN? 
Ames Aeronautical Laboratory, N.A.C.A. 


SUMMARY 


An analysis is presented which enables the temperature pro 
files, velocity profiles, heat transfer, and skin friction to be calcu- 
lated for laminar flow over a two-dimensional or axially sym- 
metric surface without pressure gradient but with an arbitrary ana- 
lytic distribution of surface temperature. The general theory is 
ipplicable to a gas of any Prandtl Number, although the numeri- 
cal results given herein have been computed for air (Pr = 0.72). 
The predictions of the theory for the special case of constant sur- 
face temperature are compared with the calculations of Crocco. 
On the basis of this comparison, it is inferred that the present 
theory enables heat-transfer and skin-friction calculations accur- 
ite to within about 5 per cent to be made for flight conditions 
up to Mach Numbers near 5 and to within about 1 or 2 per cent 
for supersonic wind-tunnel conditions up to considerably higher 
Mach Numbers. . 

A particular effort has been made to present the results, which 
re simple considering their generality, in a form that can be used 
readily in practical applications. From the mathematical point 
of view, the theory is applicable to an arbitrary analytic distribu- 
tion of surface temperature, but in any given practical case it is 
necessary that the surface-temperature distribution be approxi 
mated by a polynomial. The only unknowns in the final equa 
tions developed are the coefficients of this polynomial, so that the 
work involved in applying the theory in any given case depends 
entirely on the work involved in approximating a given surface 
temperature distribution by a polynomial 

An example is worked out in detail which illustrates some of the 
principal effects of variable surface temperature. It is shown 
that both positively infinite and negatively infinite heat-transfer 
coefficients can occur. The anomaly of infinite and negative 
heat-transfer coefficients is discussed and attributed to the cus 
tomary definition of the heat-transfer goefficient, which is shown 
to be fundamentally inappropriate for flows with variable surface 
temperature. In the particular example considered, a conven 
tional method for calculating the net heat transferred yields com- 
pletely incorrect results. A brief qualitative discussion of the 
possible effects of the heat transfer on flow separation is given. 
In order to facilitate the use of the results, all of the principal 
equations developed are collected and summarized in the section 


‘ntitled ‘‘Practical Use of Results.”’ 


INTRODUCTION 


R' CENT THEORETICAL INVESTIGATIONS of the stabil- 
ity of laminar boundary-layer flow have shown 
that it is possible for extensive regions of laminar flow 

Presented at the Fluid Mechanics Session, Seventeenth Annual 
Meeting, I A.S., New York, January 24-27, 1949. 

“ High Speed Research Division. 

t Full Seale and Flight Research Division 


to be maintained at supersonic velocities even though 
the Reynolds Number is extremely large. The calcu- 
lations of Lees’ * indicate, for example, that in the ab- 
sence of adverse pressure gradients the cooling due to 
radiation alone at a Mach Number of 3 and at 50,000 
ft. altitude can be sufficient to stabilize the laminar 
boundary layer regardless of Reynolds Number.  Al- 
though these predictions of the stability theory have not 
as yet been quantitatively verified by experiments, 
there is little doubt as to their qualitative correctness. 
Consequently, the characteristics of the laminar bound- 
ary layer have recently assumed a much greater degree 
of importance to the engineer than was generally the 
case in the past. The maintenance of extensive laminar 
flow is, of course, highly desirable in order to minimize 
the friction drag as well as the cooling requirements re- 
sulting from aerodynamic heating. 

The practical applications wherein the laminar bound 
ary layer appears to be of most importance are those 
flows that simultaneously involve both heat transfer 
This complicates a theoreti 
cal analysis in First, the well-known 
boundary-layer equations for an incompressible fluid 
are inapplicable, and the more complicated equations 
for a compressible fluid involved. 
Second, the heat-transfer and stability characteristics 
of the boundary layer are intimately connected with 
that is, with both the velocity and tempera 
This means that 


and high Mach Numbers. 


two ways. 


are necessarily 


each other 
ture profiles in the boundary layer. 
temperature and velocity profiles cannot be considered 
separately and that imprudent assumptions about either 
of these profiles cannot justifiably be made. Thus, 
velocity and temperature profiles have to be solved 
simultaneously (at least indirectly) in each particular 
problem wherein both heat transfer and high Mach 
Numbers are involved. 
Several studies of the characteristics of laminar 
boundary layers in the presence of heat transfer have 
been made, but all are subject to rather severe limita- 
tions. For the particular case of constant surface tem- 
perature, a complete set of calculations for both velocity 
and temperature profiles was made by Crocco* * and 
to some extent also by Hantzsche and Wendt.*®° In 
many practical applications, however, the surface tem- 
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perature is not constant, and appreciable errors are 
introduced by basing heat-transfer calculations on 
values obtained under the assumption of constant sur- 
face temperature. In an attempt to evaluate these 
errors, several investigations have been made of certain 
types of surface-temperature distribution.’~* These 
analyses, however, are too restricted to be of use quan- 
titatively in the problems that at present appear to be 
of considerable importance in aerodynamics. One 
severe limitation of these analyses is that they are appli- 
cable only to low-speed gas flows. Moreover, in each 
of these investigations certain assumptions have been 
made which are admittedly only rough approximations. 
As a result the existing analyses are not entirely satis- 
factory even for low-speed flows. 

As far as high-speed flows are concerned, probably the 
least accurate assumption that has been made in the 
previous investigations is that the velocity profile does 
not change with different rates of heat transfer through 
the boundary layer. This, of course, cannot be the 
case, since cooling and heating have markedly different 
effects on the velocity profiles, particularly at high Mach 
Numbers. Analyses that assume an invariant velocity 
profile are, therefore, restricted to low rates of head 
transfer as well as to low velocities. In references 7, 8, 
and 9, different types of invariant velocity profiles and 
different surface-temperature distributions were used, 
and, hence, it is not surprising that some of the quanti- 
tative results of these investigations are not in complete 
agreement. Nevertheless, these investigations are 
qualitatively consistent with each other and with the 
present investigation in indicating a pronounced effect 
of a variable surface temperature on the heat-transfer 
characteristics. 

The purpose of the present investigation is twofold: 
first, to avoid the restriction to low-speed flows and low 
rates of heat transfer which is inherent in existing theo- 
retical analyses; and, second, to obtain a solution ap- 
plicable to an arbitrary surface-temperature distribu- 
tion. It is emphasized that the present investigation 
makes no assumptions about the velocity profile. 
Likewise, no restrictive assumptions regarding the dis- 
tribution of surface temperature are made aside from 
those that are necessary in order for the usual boundary- 
layer theory to be applicable. The solution given 
herein, therefore, enables the velocity profiles, tempera- 
ture profiles, and, hence, heat-transfer characteristics 
to be calculated for a compressible laminar boundary- 
layer flow without pressure gradient but with an arbi- 
trary analytic distribution of surface temperature. 

Although the analysis is formally worked out for two- 
dimensional flow, the end results can also be applied 
to bodies of revolution by using the general relationship 
between axially symmetric and two-dimensional bound- 
ary-layer flows given by Mangler in reference 10. It is 
to be noted that the results of Mangler are valid 
for boundary-layer flows with variable surface tempera- 


ture. 
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The absence of assumptions (in the present investiga- 
tion) about the velocity profile and distribution of sur- 
face temperature would at first appear to hopelessly 
complicate the solution. It turns out, however, that 
the analysis yields surprisingly simple end results, 
which can be readily used. Since the mathematical 
development leading to these results is somewhat 
lengthy, it will be advantageous to outline briefly in 
the remainder of this introduction the more salient 
features of the analysis by means of which the simple 
yet general results are obtained. 


The first salient feature consists of the use of the 
stream function y, instead of the normal space coordi- 
nate y, as one of the independent variables, thereby 
transforming the boundary-layer equations to a much 
simpler form. This transformation was first used for 
incompressible flow by von Mises! and for compressible 
flow by von Karman and Tsien.'* It has since been 
used in obtaining solutions to other compressible flow 


problems.” ** 


A second and less essential feature of the analysis 
consists of the use of a parameter C in the assumed 
linear relationship between viscosity ratio and tempera- 
ture ratio (u/u. = CT,/T.). This parameter was 
employed in the analysis of laminar mixing given in 
reference 13, wherein its use was shown to yield accu- 
rate laminar-mixing velocity profiles at Mach Numbers 
up to 5. The merit of using such a parameter stems 
from the fact that it does not complicate the analysis, 
yet for any free-stream temperature it allows the ap- 
proximating equation for the viscosity to be more ac- 
curate in the region of boundary-layer flow near the 
surface rather than in the less important region near 
the free stream. As is discussed later on, the additional 
flexibility provided by the parameter C allows calcula- 
tions of velocity and temperature distributions, accu- 
rate to within a few per cent, to be made at high Mach 
Numbers over an unusually wide range of free-stream 
temperatures. This should prove useful in comparing 
the laminar boundary-layer characteristics obtained in 
high-Mach-Number supersonic wind tunnels with those 


to be expected in flight. 


SYMBOLS AND NOTATION 


coefficients defined by Eq. (56) 


c = factor of proportionality in the equation 
ght. = C/T a) 

$8 = the factor V 7,/T.,.(T.. + 5S)/(Tw + S) (a 
function of x) 

Cag = section skin-friction drag coefficient of a plat 
of length L[ fo” + dx/(1/2)p.U~°L] 

( = specific heat at constant pressure 

h = local heat-transfer coefficient [¢/(Ty. — Te ] 

h = local heat-transfer coefficient for the special 
case of constant surface temperature 

k = thermal conductivity 

[ = characteristic length of body 

M = Mach Number 


Number (hx/k o) 


local Nusselt 








restiga- 
of sur- 
elessly 
rT, that 
results, 
matical 
newhat 
efly in 
Salient 


simple 


of the 
coordi- 
hereby 
much 
sed for 
essible 
> been 
le flow 


lalysis 
sumed 
npera- 
Tr was 
yen in 
accu- 
mbers 
stems 
alysis, 
1€ ap- 
re ac- 
ir the 
| near 
tional 
Icula- 
accu- 
Mach 
tream 
yaring 
ied in 
those 


uation 


1 plat 


T.)] 


special 








TEMPERATURE AND 


local Nusselt Number for the special case of 


Nu = 
constant surface temperature (ix/k) 
b = pressure 
Pr = Prandtl Number of the gas (cpu/k), 0.72 for air 
q = local rate of heat transfer from the surface, 
j per unit area 
Q = total rate of heat transfer from a surface of 
length L, per unit width 
Re; = Reynolds Number based on the length 
x(U.x/v..) 
Re = Reynolds Number based on the length 
iAU_b je) 
S = Sutherland’s constant 
r = absolute temperature ‘ 
re, = average wall temperature (fo” T.dx/L) 
iz = velocity outside the boundary layer 
u, v = velocity components in x and y directions, 
respectively 
x, ¥ = coordinates parallel and perpendicular to sur- 
face, respectively 
= mass density 
= coefficient of viscosity 
7 = ratio of specific heats 
= coefficient of kinematic viscosity (u/p) 
= thickness of boundary-layer velocity profile 
bt = thickness of boundary-layer temperature pro- 
file 
w = exponent in the equation u/u. = (T/T,)® 
T = shear stress [u(0u/dy)] 
y = stream function 
~ = dimensionless independent variable occurring 
in the Blasius function f(y) [this variable is 
related to y by Eq. (42)] 
f(n = the solution to Eq. (17) (the Blasius equa- 
tion) 
V,(n), ¥n(n) = characteristic functions of the problem, de- 
fined as the solution to Eqs. (33) and (44), 
respectively, for a given value of n 
Vin = function defined by Eq. (28) 
r(n), #(n functions defined by Eqs. (25) and (48), re- 
spectively 
~ = order of magnitude 
Subscripts 


ee = conditions outside of the boundary layer 


= conditions at the wall or surface of the body 
= equilibrium conditions of an insulated sur- 
face 
Superscripts 
= dimensionless quantities as defined in Eq. (11 
= ordinary differentiation 


Basic EQUATIONS AND ASSUMPTIONS 


Basic Equations 


In the classical boundary-layer theory, temperature 
gradients, as well as velocity gradients, in the direction 
of flow are neglected in comparison to those perpendicu- 
lar to the direction of flow. (See Fig. 1.) It is the 
purpose of the present investigation, however, to ac- 
count for large temperature gradients in the direction of 
Consequently, the first step in the analysis is to 
examine the basic differential equations of laminar 
motion of a viscous compressible fluid in order to see 
what restriction on the magnitude of the temperature 
gradient along the wall is necessary in order to be con- 


flow. 
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Figure | —Schematic dramng of the flow. 
sistent with the approximations of the usual bound- 


ary-layer theory. Since these considerations are some- 
what lengthy, the details are presented in Appendix A. 


It will suffice here to state that the analysis 
given therein shows that the usual boundary-laver equa- 
tions, * 
Ou Ou O Ou\ momentum 
pu T+ - u ; 
ox ov ov ov equation =| 
: " . | 
fe) 4 continuity 
- (pit) - (pv) = i 
ae ay equation 
: (1) 
O(c,T’) O(c,T) 
pu —~ — + pu * 
Ox Oy 
o («S) -— (S*) energy 
Oy \ Ov Ov equation 


are applicable to flows with large temperature gradients 
in the x direction provided 


(OT /O0x), ~ (T, — T.)/6 or less (2) 


i“ 


At supersonic Mach Numbers the surface temperature 
T,, is normally of the same order of magnitude as the 
equilibrium temperature T, = T.(1 + 0.169M.?), 
and the boundary-layer thickness is approximately 
6L/V/Re. 


comes 


Under these circumstances condition (2) be- 


oT Ta) 


~ 0.03.1 .2+/ Re or less 2a) 
O(x/L) 


This means, for example, that at a Mach Number of 
1 and a Reynolds Number of 10° the maximum surface 
temperature gradient consistent with the usual bound- 
ary-layer approximation is equivalent to a linear gradi- 
ent that would produce a thirty-fold increase in T/T, 
over the length of the surface. Thus, for high-speed 
flows, condition (2) does not impose any restrictions of 
practical importance on the present theory for reason- 
ably smooth and continuous surface-temperature dis- 
tributions. For low Mach and Reynolds 
Numbers, however, condition (2) may in some cases 
restrict the applicability of the theory. 


subsonic 





* Throughout this report all equations are written in absolute 
units; hence, the acceleration of gravity and Joule’s constant do 


not appear in the analysis. 








550 JOURNAL OF THE 





5.0 

















20 











Figure 2.—Variotion of viscosity with temperature. 


Assumptions 


If no additional assumptions are made other than 
those implied in condition (2), the problem under con- 
sideration then involves the solution of the three simul- 
taneous partial differential equations [Eqs. (1)] for the 
three dependent variables u,v, 7. It is apparent that 
some additional assumptions are needed in order to 
make Eqs. (1) amenable to solution. The important 
assumptions, including those already made, which are 
used in obtaining a solution to the basic equations of 


conservation of momentum, mass, and energy are: 


i) 6/L < < 1 (usual boundary-layer approximation), 
ti) (O72 /Ox),, ~ (T, — T.)/6 or less. 

(iii) Pr = c,u/k = constant (0.72 for air). 
constant. 


iv) ¢, = 
CT/T.., where C can depend on 7,, 


V) B/ be = 
7 .,and M,.. 


Assumptions (i) and (ii) are discussed in Appendix 
A. For normal aerodynamic applications these two 
assumptions do not restrict the usefulness of results 
based upon them. Likewise, assumptions (111) and 
(iv) do not unduly restrict the applicability of the 
theory, since c, and Pr for air are constant to within a 
few per cent over a wide range of temperature. These 
first four assumptions are usually made. Assumption 
(v), however, is not made so often and requires some 
discussion. 

It is known that an accurate equation for the varia- 
tion of viscosity with temperature is given by Suther- 


land’s formula 
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u/pe = (T/T.) “(Te + S)/(T + S)) (3) 
where S is Sutherland’s constant and is approximately 
216°F. absolute for air. Curves representing this 
formula are shown by solid lines in Fig. 2 for two free. 
stream temperatures that are considered to be represen- 
tative of the two extremes encountered in aerodynamic 
applications. One free-stream temperature is 18S8°F_, 
which is typical of the high temperatures encountered 
in the atmosphere near 160,000-ft. elevation. The 
other free-stream temperature is —3S8S8°F., which corre- 
sponds approximately to the free-stream temperature 
of a supersonic wind tunnel operating in the Mach 
Number range near and beyond 5. Also shown in Fig 
2 are the approximating formulas wy. = (T T.,)%% 
and wu. = (7/T.). It is apparent from a glance at 
Fig. 2 that no single value of w in the functional rela- 
tionship w/w. = (7/T..)° can adequately represent 
the viscosity over the wide range of free-stream tem- 
peratures encountered on passing from supersonic wind- 
tunnel to flight conditions. Of course, by using a dif- 
ferent value of w for each free-stream temperature and 
Mach Number, a sufficiently accurate description of 
the actual boundary-layer characteristics could be ob- 
tained. Since the structure of the boundary-layer 
equations is such that this latter procedure would re 
quire a prohibitive amount of work for flows where the 
surface temperature varies, it is necessary to use some 
other approximate equation. 

It has been shown in references 2, 3, and 5 that the 
assumption w = | greatly simplifies the boundary- 
layer equations yet provides a sufficiently accurate solu 
tion in many cases. This fact has also been used to 
advantage in the analyses given in references 4 and 
14. The linear approximation used in these references 
is w/po = 7 /T..,* which always gives the correct value 
for the viscosity in a region near the outside of the 
boundary layer but unfortunately gives grossly inac 
curate values in the much more important region near 
the surface. (See Fig. 2.) Asa result, the assumption 
uu. = 71/7. does not predict satisfactorily all of the 
important boundary-layer characteristics at high Mach 
Numbers and high free-stream temperatures. This 
shortcoming of the usual linear equation for viscosity 
can be largely overcome by the simple artifice of adding 
to the linear approximation a factor C' that can depend 
explicitly on the free-stream temperature as well as the 
average surface temperature. Thus, regardless of 
Mach Number, free-stream temperature, or average 
surface temperature, the approximation 


u/te = C(T/T.Q) 4 
can be made to yield nearly correct values for the vis 
cosity in the important region near the surface. If, 

* Hantzsche and Wendt use p/p = T/T y in reference 5 and 
Near the 


ipply this approximation at all steps in their analysis 
T/T, so that division o 


free stream this becomes p./uy = 
these two equations yields u/u, =7/T,. 
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TEMPERATURE AND 


jor example, the approximate Eq. (4) is matched 
with Sutherland’s Eq. (3) at a temperature equal to 


the average surface temperature 7,, then C is given 
bv 


C=VT./TATe+S/(To+ S16) 


This constant is not to be confused with the variable 
Ce = (ue/uo)(T2/T), which is introduced later. 

It is shown in the discussion of results that, even at 
Mach Numbers near 5, the use of Eq. (4) [assumption 
v)] in simplifying the differential equations does not 
limit the accuracy of the final results by more than 
about 5 or 6 per cent over the entire range of free-stream 
temperatures encountered both in extreme supersonic 
wind-tunnel and flight conditions. This is not surpris- 
ing in view of the fact that the very approximate as- 
sumption w/w. = T/T. (Fig. 2) predicts the skin- 
friction and heat-transfer characteristics to within 
about 15 per cent at 17. = 5 and to within a smaller 
percentage at lower Mach Numbers. 


SOLUTION OF EQUATIONS 


Method of Solution 


Using the assumptions that c, and Pr are constant, 
the basic Eqs. (1) can be written as 


Ou | Ou re) ( Ou ) 
pul + pv = Lh 
Ox Ov oy Ov 


momentum 


y) 
pu) + (pv) = 0 continuity 
Ox j 
(6) 
o7 ol 1 oO ( o7 ) 
pul + pv =~ u T 
Ox Oy Pr oy Oy 


(o“) 
iv energy 
oy j 


This set of partial differential equations can be 
greatly simplified by transforming them to a new coor 
dinate system. As is indicated in the introduction, the 
transformation that is used employs the stream func 
tion Y, instead of y, as one of the independent variables. 
By using y as one of the independent variables, the con- 
tinuity equation is automatically satisfied, and the 
velocity components are given by 

Pao OY Da OY " 
“= , v ‘ 
p Oy p Ox 
The formulas for making the transformation from the 
independent variables (x, y) to an arbitrary and com 


pletely new set (s, ¥) are 
Oo Os O Oy O 
Ov Ov Os Oy Oy 


re) Os O Oy O 


Ox ry Ox Os Ox Oy 
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Using Eqs. (7) and letting s = x givest 


0/Ovy = (pu/p.)(0/dp) 
(0/Ox), = (0/0x), — (pv/p.)(0/Op) 
Hence, 
pu(O0/Ox), + pv(0/ Oy), = pu(d/Ox), 
and by using the relation 
(9) 


pu = Cpoau 


which is a consequence of the assumption that yu is pro- 
portional to the temperature, then the basic Eqs. (6) 
Ox Po OW 


( 2 a 

Cu momentum 
OY 

oT ; Me re) 


Lp (cy?) 4 
ox Prp.. Ow . oy 


baCu fou\? 
energy 
Cro \OW 


The above system of equations can be further sim- 


become 


Ou Me 


(10) 


plified by introducing the dimensionless variables 


a =e. GP me BFE 
T* = T/T... p* = p/pa (11) 
u* = w/w, ¥* =y/V ».U,LC 


Using these dimensionless variables, Eqs. (10) become 


ou* o , ou* P 
— oy* u oy* momentum (12a) 
or™ 1 O ( 27) 
= — ———1 w 
ox* Pr oy* oy* 
: - ou* 
(y — 1)M.*u* oy" energy (12b) 


for which the boundary conditions are 


u* = 0 at y* = O(y = 0) 
g° = lat ?* = oe) (y = 0 ) 
T* = T,,* aty* = 0 
I* = laty* = « 
The dimensionless surface temperature 7),* is, of 


course, a function only of x. 

It is important to note that the momentum Eq. 
(12a) and the boundary conditions on u* are completely 
independent of the energy Eq. (12b) and the boundary 
conditions on 7*. This is, in fact, the fundamental 
point of the entire theoretical analysis, since it is only 
through the independence of the momentum equation 
from the energy equation in the coordinate system em 
ployed that a simple solution to the problem is possible. 
It also may be noted that, if the viscosity were not 
taken as proportional to the first power of the tempera- 

t The variables held constant in a differentiation process ar« 


explicitly indicated in those cases where ambiguity could result 


if the subscript notation were not used 
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- = | ] [/ | or 
i el a oe aS TT . T Sf’ +f" =0 (17) 
/e-—_—_—_ |_| A 2 a SS a ‘ ; is Pon ° e . ar ee 
| {Al yg | | which is the familiar equation of Blasius satisfying the 
Mou ae tt} 4 . boundary conditions 
Pe ee ee ek 2p Biden 
| il } | | f(2) =2 
ae ete f'(0) = f() = 0 
a ie aan Tan a oT a .. . 2. . ef These boundary conditions are, of course, consistent 
—— t t—A—_+ tft — with those imposed on Eq. (12a). Since complete 
12}-——_+_ + _4__ + +1 1 1 | | | tables of the functions f, f’, and f” can be found in 
ere EIA CI many places (e.g., reference 16), the numerical solution 
oe } | to Eq. (12a) is known. Curves of these functions are 
a a ar ame i ae a shown in Fig. 3. 
ST aESERIERAEEEEEREEEEE Gao . a ae The continuity and momentum equations are solved 
8 | ee a hed by Eqs. (7) and (13), respectively. Consequently, the 
only equation now remaining to be solved is the energy 
y a a | | 1] Eq. (12b). Substituting Eq. (13) into Eq. (12b), 
6 + + + “ + 4 
o7* 1 oO (¥ oT*\ y—1M.?(f" 
aL L { Ox* 2Pr dy* ay iss. ™ 
ee A | as a) Because f’ and f” are functions of 7 alone, the independ- 
ent variables (x*, ¥*) are changed to (x*, 7). This is 
od a aa i x done by the transformation formulas 
1 or. 7 a ae 
NI (2) -(2) +22-(S) 
°O 5 10 is 20 25 30 Ox*] On" 7. Ox* On ox. 
° f oO 
Figure 3 Variation of Blasius functions f, f' and f" with are 19 
boundary-layer characteristic variable 9. ae 
ra) ) _ om Oo 1 a) 
Fe 2 OF On -Wx*f! On 


ture, the density p* (and hence 7*) would appear in 
Eq. (12a). This, of course, would excessively compli- 
cate matters even in the (x, Y) coordinate system, since 
the equation representing conservation of momentum 
would not be independent of the energy equation, and, 
hence, the energy equation would not be a linear differ- 
ential equation 

The fact that the density p* 
momentum equation means that the function that 
solves Eq. (12a must be the well-known Blasius func- 


does not appear in the 


tion for incompressible flow. Hence, 


u* = (1/2)f'(n) (13) 
. ke 
where 7 is defined by 
f(a) = V°/V' x" (14) 


solves Eq. (12a) can easily be 


Thus, differentiating 


The fact that Eq. (13 


verified by direct substitution. 


Eq. (14 
i * 
See = — J (15) 
Ox” ee xn* 
f'(On, OW*) = 1//x* (16) 


e . 
and‘substituting into Eq. (12a) gives 
4 


Substituting Eq. (19) into Eq. (18) gives the following 
final form for the energy equation 


oT T* oT* Pr 
: : + Prf ne 2Prf’x* o— x 
On? “ On : on” 4 
(y — 1)M,.? (f’ 20 


with boundary conditions on the function 7*(x*, 7 
which can be written as 


T*(*, 0) = T,* + He" 
T*(x*, o) = 1 


where 7.* is the dimensionless equilibrium temperature 
that would exist if the surface were insulated and ¢(x* 
is an arbitrary function that specifies the surface- 
temperature distribution. Eq. (20) is a linear, second 
order, inhomogeneous partial differential equation (of 
parabolic type) with inhomogeneous boundary condi- 
tions. The standard mathematical methods of sepa 
rating variables are thus applicable for its solution. 

It may be noted that the transformation from (x, ¥ 
to (x, ¥) is not the only transformation that produces 
the same end result [Eq. (20)]. The transformation 
introduced by Crocco* and also used by Hantzsche and 
Wendt! will lead precisely to Eq. (20). Likewise, the 
transformation used recently by Howarth" will also 
make the momentum equation independent of the 
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TEMPERATURE AND 


energy equation and thus eventually lead to Eq. (20) 
when applied to the present problem. 


Solution for Temperature Profile 


Particular Solution of Inhomogeneous Equation.— 
Since Eq. (20) is linear, a particular solution of the form 


T*(x*, n) = X(x*) + N(n) (21) 


can easily be found. Substituting Eq. (21) into Eq. 


(20) gives 


1 : oF, P , 
— | NV’ Pr fN’ + (y — 1)M.2(f’ ] = ¢*X’ 
zl + SF; 7,0 ) (J) 

(22) 
Hence, both sides of Eq. (22) must be equal to a con- 


stant A. This yields as a solution for x 


X = K log x* 
which will give a logarithmic singularity at.the origin 
unless A = 0. Setting K = 0 the function N(») is de- 
termined from 


N” + PrfN’ = —(Pr/4)(y — I)M.7(f")? (23) 


This ordinary differential equation is the same equation, 
apart from the constant factor (y — 1)M.’, which has 
already appeared in the investigation of Pohlhausen.'® 
The solution to Eq. (23) is obtained by the usual 
methods of solving first-order ordinary differential 
equations. It is 


N(n) = 1+ [(y — 1)/2]M.2r(n) — (y — DC gt) 
(24) 


where 


P ad g 
r(n) = f Uf” (e) J?" f [f"(0)]2-"rdedt (25) 


and 
g (n) = fir@rra (26) 
” 
The function 7(y) is plotted in Fig. 4 for Pr = 0.72 


(air). Since Eq. (24) is a particular solution to Eq. 
(20), the arbitrary constant C, can, for convenience, 
be taken as zero. It is easily verified that C; = 0 
corresponds to an insulated surface, since Eqs. (24) 
and (25) show that (ON/On),, is zero. Hence, (0) 
is the dimensionless equilibrium surface temperature 
Ree 

N(0)= T.* = a =1+4+7(0)2 . MMe? (27) 
This, of course, is the well-known expression for the 
surface equilibrium temperature in a compressible 
fluid. The “‘recovery factor” r(0) is independent of 
Mach Number and equal to 0.845 for Pr = 0.72. 
(See Fig. 4.) It may be noted that Eq. (25) immedi- 
ately gives r(0) = 1 for any gas wherein Pr = 1. 
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Figure 4.— Variation of r and F with boundary-layer 
characteristic variable 9, for B,#0.72 


Summarizing, a particular solution to the inhomo- 
geneous Eq. (20) is 
N(n) = 1+ [(y — 1)/2]M =’r(n) (28) 
which satisfies the boundary conditions 
T*(x*, 0) = T7,* [((ON/On)» = 0) (29) 
T*(x*,0) = 1 
Solution to Homogeneous Equation.—The homogene- 
ous part of Eq. (20) is 
o°7* o7* o7T* 
+ Pr — 2Prf'x* — =0 30) 
On? f On f ox* ( 


Solutions to this equation can be obtained by separation 
of variables in the form 


T* = X(x*)T(n) 
which, when substituted into Eq. (30), gives, 
(1/f’Y)(¥" + PrfY’) = 2Prx* X’/X (31) 
Hence, both sides of Eq. (31) must be equal to a con- 
stant, say 2Prn. Then 
X ,(x*) = x*" (32) 
and Y,,(») is determined from the equation 
Y,” + PrfY,' — 2Prf'nY, = 0 (33) 
The boundary conditions most convenient to impose 
on the functions J’,(m) are 


Y,(0) = 1) 


Y,(2) = of (34) 


Since the homogeneous equation is linear, solutions of 
the form a,x*"Y,(n), where a, is an arbitrary constant 
can be superimposed to give as a general solution 


> anx*"Y,(n) (35) 


n=0 








D4 JOURNAL OF THE 
which satisfies the following boundary conditions on 
fii (cal 


T*(x*, 0) = > a,x** 
0 


T*(a*, ©) = 9Q 


These functions }’,(y) that solve the homogeneous part 
of the energy equation have been calculated numerically 
for values of n = 0, 1, 2, 3, 4,5, and 10. The method 
of numerical computation used is described in Appendix 
B, and curves of the functions Y,(n) are shown in 
Fig. 5. 

Complete Solution. 
complete solution is obtained by adding the particular 
solution (24) of the inhomogeneous equation to the 
general solution (35) of the homogeneous equation, 


Since Eq. (20) is linear, the 


T* (x*, n) = N(n) + Yoanx*" V,,(n) (30 
0 


The complete solution satisfies boundary conditions 
equal to the sum of the boundary conditions imposed 


on Eqs. (23) and (35)—-namely, 


T .* = T*(x*, 0) = T.* + > anxr™ (37) 
T*(x*, ©) = ] 


Comparing this with the boundary condition for an 
arbitrafy surface-temperature distribution, which, for 
convenience, was written in the form 


Lg =o iis") (38) 


it is obvious that the coefficients a, are determined by 
the arbitrary function ¢(x*) through the relation 


t(x*) = }oa,x* 


0 


(38a) 


Thus the coefficients a, are simply the coefficients in the 
Taylor series for t(x*), so that if (x*) is known analvti- 
cally, 

a, = t"’(0)/n! (39) 
In practical applications, however, Eq. (38) would be 


replaced by a finite sum 
proximated by a polynomial 


that is, ¢(x*) would be ap- 
the coefficients a, of 
which are the coefficients appearing in the final Eq. 
(36) for the temperature distribution. 

Although the temperature distribution is given in 
terms of yn, the corresponding value of y can be found 


ia simple manner. Thus, Eq. (7) can be written as 


¥ - fo. pu)dy (40) 


where the integration is carried out with x held con- 
stant. Substituting Eqs. (13) and (14) into Eq. (40) 


gives 
(y/2)V U./».xC = - 1*dn (41) 
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Figure 5. —Variation of Yn and Yq with boundary-layer 


characteristic variable 7, for F =0.72. 


This can be integrated in closed form simply by sub- 
stituting Eqs. (36) and (28) into the integrand of Eq. 
(41) and then integrating term by term 
F a U. : ote I <2 : kn 

| =nt+- M.2r(n) + Do anx*"V,(n) (42 
Z voxC 2 0 


The fixed functions ry and Y, are the integrals of r and 


Y’,, respectively, 


> 


” 
r(n) = / r(n)dn +3 
*” 
[vane tH 


Curves of these functions are given in Figs. 4 and 5 so 
that the value of y corresponding to a given value of 9 
can be found readily from these curves and Eq. (42). 


Il 


Y,(n) 


Heat Transfer 


The convective heat g transferred from a unit area of 
surface per unit time at a given point is 


gq = —k(OT/2y) w (45) 


Although this equation does not include a term ac- 
counting for the heat radiated from the surface, it is 
an easy matter to do so by the conventional method. 
The equations that follow, therefore, represent the heat 
transferred through the boundary layer by convection 
and conduction. 

Differentiating Eq. (36), substituting into Eq. (45), 
and noting that (Or/0n), = 0, 


re) = — , 
q= ~ eT .(3") > aug ¥ (6) (46) 


By using Eqs. (7), (13), and (16), On/Oy can be evalu- 
ated to give 
On On dy \ ym 4 


coat pets : (47 
¥.eC 2T* 


dy wWdy Vy.U.LC 
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hence, 
On lu 
) YE 
OV) voxC 2T,,* 
The ratio k,/7,,* can be written as 
ky w @ m 7 
«ft 2 feb. = Ceke (49) 


i ios x (es Rar 


where, according to Eq. (3), 


CoH pu®/To* = VT o/Tol(To + S)/(Tw + S)] (49a) 


The variable C, is a function of x since 7,, depends 
upon x, and is not to be confused with the constant C 
defined by Eq. (5). For the special case of constant 
surface temperature C, is equal to C. The rate of heat 


transfer is 


Re i. ; l ; i rr \ 
ne Ce y = L a Ya n (O) (50) 
2 78. 6 


Usually g is expressed as the product of a heat-transfer 
coefficient and the difference between the surface tem- 
perature and the equilibrium temperature. The usual 
expression, however, involves a constant surface tem- 
perature, so that a slightly different expression will have 
to be used for the case of a variable surface tempera- 
ture. Actually, there are many possible definitions 
that will reduce to the customary one in the special 
case of constant surface temperature. The most 
logical equation to choose for defining the heat-transfer 
coefficient h/ is 

q = h[T(x, 0) — T,.] = h(T, — T,) (51) 
where 7(x, 0) is the local surface temperature 7),. 
Eq. (51), of course, is identical to the customary equa- 
tion if the surface temperature 7, is constant. Com- 
paring Eqs. (50) and (51) and noting that the leading- 
edge surface temperature is 7*(0, 0) = 7,* + apo, it 
follows that 


j [T*(0, 0) an ig oe 
, a 4 
(Te sil T,*)C 
VU. /(v0xC) S agxt Y 


T2= F9 a 


h= 


Rules 
9 x 


» (O) (52) 
where / is the heat-transfer coefficient for the case of 
constant surface temperature and is given by 


h — (k./2)V U.C/v.x Yo'(0) 


Thus, Eq. (52) simply gives the additional corrective 
terms necessary to account for the effects of variable 
surface temperature on the local heat-transfer coef- 
ficient. The use of the temperature 7, in the defining 
Eq. (51) for the heat-transfer coefficient does not enable 
the final expression (52) to be completely independent 
of Mach Number, as is the case for flows with constant 
surface temperature. Nevertheless, Eq. (52) is sub- 
stantially independent of the Mach Number, since the 


(53) 
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only effect of Mach Number would be its indirect ef- 
fect on T,, which indirectly would have a small effect 
on C and Cy. 


The relation between the local Nusselt Number and 
Reynolds Number is 


hx - [7*(0, 0) — T.*)C. 


Nu = Nu com or ge = 
k.. (72° ~ TO 
i V Re, 2 : A 
2VSC(Te* — T.*) x, ee oe 
where Nu = hx/k.e = —(VW/2)~VRe,Yo'(0) is the 


local Nusselt Number for constant surface temperature 
and Re, is the Reynolds Number based on free-stream 
conditions and the length x. The total rate of heat 
transfer Q per unit width of a two-dimensional surface 
is obtained by integrating Eq. (50) from x = 0 tox = L, 


a Rol Vol > anbn¥,'(0) 
2 Cre n=O 


where the coefficients 


“ 
/ Cot is *dx . 
0 


can be evaluated for any given distribution of surface 
temperature (given C,). 

The quantities Y,,’(0) that appear in all the fore- 
going equations for heat transfer are given in Table | 
for x = 0, I, 2, 3, 4, 3, and 10. 


(55) 


Q= 


(56) 


b,, = 


TABLE | 
Numerical Values for ¥n’(0) 

nN Y,’(0) 

0 —0.5915 

1 0.9775 

2 1.1949 

3 - 1.3680 

4 ~1.4886 

5 —1.5975 
2.0121 


Solution for Velocity Profile 

The dimensionless velocity is related to the independ- 

ent variable by the Blasius function 
u* = (1/2)f’(n) 


so that the velocity is immediately known for any 
given value of 7. Moreover, the value of y correspond- 
ing to a given value of 7 is obtained from the relation 


(13) 


y i y—1 7 

; =nt M.*r(n) + 

2 VvixC 2 
2, a,x*" Y,,(n) (42) 
0 

which was previously derived. Since the functions 


r(n) and Y,,(n) are fixed functions for all problems (they 
are given in Figs. 4 and 5) and since the coefficients 
a, are determined once the surface-temperature dis- 
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tribution is specified, it is therefore a relatively simple 
matter to calculate velocity, as well as temperature, 
profiles for an arbitrary distribution of surface tem- 
perature. 


Skin Friction 
By differentiating Eq. (13), the shear stress can be 
evaluated at any point in the boundary layer 
Ou I On .. 
uw = -f"(n) U. 
oy of oy 
Using Eq. (47), the shear at a given value of n becomes 


(57) 


r= f"(n)V Va/vaxC(uU./4T*) 
which reduces to 


(57a) 


Te = (f°(0)/4)UCuV 9p. U aa /Cx 


for the shear stress at the surface. Integration from 
x = 0 tox = L gives the section skin-friction drag 


coefficient 
P oa . 
Cay = [f"°(0)/V ReC](1/2) | Cydx*/Vx* (58) 
ia] 


which reduces to 


Cay = [f" (0) V Re] VC = (1.328/V Re)VC_ (58a) 


for flows with constant surface 


later part of this paper, this formula is compared with 


temperature. In a 
the results of other analyses. (See discussion of re- 
sults.) 


Relation Between Skin Friction and Heat Transfer 


Although the equation relating heat transfer and 
skin friction for the general variable surface-tempera- 
ture flows is too complicated to be of much practical 
use, it is derived here not only for the sake of complete- 
ness but also for the purpose of illustrating that the 
familiar heat-transfer and skin-friction analogy is not 
applicable to flows where the surface temperature 
varies. This conclusion deduced directly 
from the dissimilar boundary conditions imposed on 
that is, 
the 


could be 


the energy and the momentum equations 
Uu 0 at T,. varies—-however, 
following gives a quantitative indication at the relation- 


surface while 
ships. 

The ratio between local rate of heat transfer and 
local skin friction (shear) is obtained by dividing Eqs. 
(50) and (57a 

2kel . 2 : 
_- - - a,x* "[—Y,'(0)] 
7 f (O)L cdo n=0 
This can be put in a more illustrative form by evaluat- 
ing Yo’(0). Thus, it is easily verified by direct sub- 
stitution into Eq. (33) that 


Yo(n) = frre f Uf" n) |Prdn 


hence, 


Y,’(0) = 


—[f"(0) ]”" / I “Lf"(n) |Prdn 


The numerical calculations of Pohlhausen show, how- 
ever, that for 0.6 < Pr < 1.0 an approximation valid 
to within less than 1 per cent is 


J [f"(n) }?"dn = 21f"(0) 1?" ~ 1/Pr' 
so that to within 1 per cent, 
¥o'(0) = —[f’(0)/2]Pr'” (59 


Using this equation, the ratio g/r becomes 


q fel « = ‘en Y,’(0) , 
= ~ x — 7 (60 
7 Upp be Oat" Fg) 


This equation shows that there is no simple relation 
between heat transfer and skin friction valid for flows 
with variable surface temperature. For the special 
case of constant surface temperature, however, Eq. 
(60) reduces to 


Tl? ry ; r 2, . 
g/t = Chl —_ io/Uaukt (61) 
which is the same relation as that given by Hantzsche 
and Wendt.° 

The ratio between total rate of heat transfer and the 
integrated skin friction also reduces to a simple expres- 


sion for the particular case of constant surface tempera- 


ture. Thus, dividing Eq. (55) by Eq. (58a) and re- 
membering that for constant surface temperature 
a, = T,* — T,*, 
) asV Re 
Ha = OS _¥)] 62 
ko(T» — T.)V Re f’(0) 


Substituting Eq. (59) into Eq. (62) gives, finally, 
Q/Ro(T — Te)V Re = (Ca Re/2)Pr 63 
This expression also is applicable only to flows with 
constant surface temperature. 
PRACTICAL USE OF RESULTS 


Statement of Complete Solution to Problem 


The foregoing solution to the laminar boundary- 
layer equations has been obtained for an arbitrary sur- 
face-temperature distribution, which for convenience 
is expressed in the dimensionless form 


i le/le = i,* + is” 38 
where 


M. 


| ig : 1 + r(0) 2 — 
0.72 


1 + 0.169 \/ =? (for air with Pr = 


ll 


is the dimensionless equilibrium temperature and ¢(x*) 
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TEMPERATURE 
js an arbitrary analytic function of x* =x/L. Since 
any regular function such as t(x*) can be expressed in a 
power series (or approximated by a polynomial for 
practical applications), the coefficients a, in the Taylor 
series 

t(x*) (38a) 
may be considered as known for any given surface 
temperature distribution. It is important to note that 
the final results of the analysis, which are summarized 
below, involve as unknowns only these coefficients a,. 
The work involved in applying the theory, therefore, 
depends only on the work involved in approximating a 
given surface temperature distribution T,,*(x) by a 
If ¢(x*) is given in analytic form, then 


a, = t'"(0)/n! (39) 


polynomial. 


Many of the equations summarized below are ex- 
pressed in terms of the characteristic boundary-layer 
variable n, which is related to the physical coordinates 
through the relation 


y | : — | 
5 | a fee py M. r(n) + - c* Y,,(n) 
« VnXC 2 n=0 


(42) 


This relation allows those equations that are expressed 
in terms of » to be converted to physical coordinates, 
since r(n) and Y,,(n) are known functions (Figs. 4 and 
5). The principal equations developed in the course 


of the analysis are: 


ee. ee r 5 1 r(m) + 
I*(e*, 9 Tr =1+ T.°r(n) 4 
> a,x*"¥,(n) (36 
‘i u Ls : 
(xs, @) 7a = 5 f(a) (15 





oT k + 7 a w 
q = —k,, ( -) = a 
Oy / w Cr, 


S apx*"Y,'(0) (50) 


See ane tt z [Tw*(0, 0) — Te IC, 
Ty ro T, p ge" -_ i; ™C 
> aax*"¥,,"(0) 
RaoCee Ua. as . a , (0 (52) 
2 Cy.x T,* — T,* 
h = —(2., 2)VU.C yx Yo'(0) (53) 
, hx — "ty ” Tv." (0, 0) a 7. 
Vu Nu —= ; = 
k, rd Py a 2 tg 
~ [7 a,x*"Y,,'(0) 
é, +> © ow 


2 = r.* — T,"*) 
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Vu = hx/ke = —[Yo'(0)/2]VRe,C 


l 
cay = [f"(0)/V ReC](1/2) f (Cydx*/V/x*) (58) 
0 
Cay = [f"(0) V/RelV/C ) (for flows with 


>» constant surface 


(58a) 
= (1.328/¥y ‘Re) ~ iy \ temperature) 


It is to be remembered that these equations are for two- 
dimensional flow but may be applied to axially sym- 
metric flows without pressure gradient by making the 
modifications indicated in reference 10. 

Numerical values of Y,’(0) applicable for all flows 
with variable surface temperature are given in Table | 
Curves of the fixed functions f’(n), f’(n), r(n), r(n), 
Y,(n), and Y,,() appearing in the above equations are 
shown in Figs. 3, 4, and 5. The expression 7,,*(0, 0 
is, of course, the leading-edge surface temperature. 
The variable C, appearing in these equations is defined 
by the relation 


Coe . Mw*/ a (49) 


Using Sutherland's formula [Eq. (3)], C, is given by 


Cy = VT u/To[(Te + S)/(To + S)] 


The constant 


C= VT,,/T.[(T. + S)/(Tx + S)] (5) 
matches the viscosity as given by the approximate Eq. 
(4) with Sutherland’s formula at the particular tem- 
perature 7,. For general engineering use, Eq. (5) 
should be sufficiently accurate to account for the tem- 
perature level of a particular boundary layer. 

Once the surface temperature has been approximated 
by a polynomial, the procedure for calculating a given 
involves only the substitution of the coef- 
appropriate 


quantity 
ficients (a,) of the polynomial into the 
equations and the use of the appropriate curves in 
Figs. 3, 4, and 5. 
In practical applications, the local Nusselt Number 
\u for constant surface temperature, 
’(0) , 

Nu = VRe,C = 0.296 VRe,C 
could, if desired, be replaced by some other expression 
that explicitly includes the effect of Mach Number. 
(The above expression only implicitly includes the ef- 
fect of Mach Number through the constant C.)  Al- 
though such a procedure would possibly improve the 
accuracy slightly, a rigorous justification of it has not 
been attempted. The most complete set of calcula- 
tions available to date for Nu has been made by 
Crocco.* 4 

The foregoing equations are applicable if the dis- 
tribution of the rate of heat transfer is prescribed rather 
than the distribution of surface temperature. For 


such cases, Eq. (50) determines the coefficients a,, 
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Figure 6. —Skin friction drag of an insuloted plate using 
oitferent variations of viscosity with temperature. 


which, once determined, can be substituted into the 
remaining equations to calculate the temperature dis- 
tribution, velocity distribution, and other quantities of 


interest. 
DISCUSSION 


Reduction of Present Theory to the Solutions of Blasius 
and Pohlhausen for Incompressible Flow over an Insu- 
lated or Heated Plate 
For flow over an insulated plate, the surface tempera- 

ture is equal to the equilibrium temperature 7, and all 
coefficients a, are zero [Eqs. (38) and (38a)]. More- 
over, for low-speed subsonic flows, the absolute tempera- 
ture within the boundary layer never departs appre- 
ciably from the free-stream temperature, so that C -> |] 
as M.,—» 0. Therefore, the velocity distribution is 


u = (1/2)f'(m) 


where the relation between and y [Eq. (42) | reduces to 
jens 
(y/2)V U./vox = 7 


This is the well-known solution of Blasius for incompres- 
sible flow. 

Writing the free-stream Mach Number as U,,/- 
Vy RT.., Eq. (36) for the temperature distribution be- 
comes 

T/T. = 


1+ [(y — 1)/2)(U.?/yRT.)r(n 


or 
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(T — T.)/To = 7(y)(Ua?/2yT @) 
which is a solution of Pohlhausen.!® The 
theory, therefore, contains as special cases the familiar 
solutions of Blasius and Pohlhausen for low-speed flow 


present 


over an insulated plate. 

Another low-speed solution for heat transfer from a 
temperature plate to an incompressible 
boundary layer is also a particular solution of the pres- 
Thus, as V,, — 0, C — 1, and Eq. (33 


constant 


ent analysis. 
becomes 


Nu = 0.296~/ Re 


which is identical with the expression obtained by Pohl- 
hausen for Pr = 0.72. 


Comparison with Other Results for Compressible Flow 

Surface of Constant Temperature 

Before discussing some of the effects of a variable 
surface temperature, it will be well to compare the re- 
sults of the theory given herein with those of other 
investigations of compressible laminar boundary-laver 
flow Such comparisons must, of course, be restricted 
to cases of constant surface temperature, but neverthe- 
less will serve to illustrate the degree to which the 
present theory [based on the approximation y/u. = 
C(T/T..)| agrees with other more precise and much 
more complicated calculations of boundary-layer flow 
at high Mach Numbers. 

Skin Friction. 
skin friction coefficient of an insulated surface as de- 
termined from Eq. (58a) and also as determined from 
In this figure the various methods 


Fig. 6 shows a comparison of the 


other calculations. 
of calculation are compared for two free-stream tem- 
peratures: 188° and —3SS8°F., which are considered 
to be representative of extreme conditions in the 
stratosphere and in high Mach Number supersonic 
wind tunnels, respectively. For each free-stream tem- 
perature three different methods have been chosen for 
comparison with the present theory. One method is 
based on the assumption y,u. = (7T/T..)°%*, often 
considered to be adequate for general aerodynamic ap- 
plications. A second method is based on the simple 
assumption w/u. = 7/7. used by Crocco and Lees. 
The third method chosen for comparison is labeled 
‘“‘Sutherland”’ on the figures, and is taken from the cal- 
culations of Crocco.‘ This latter method of calcula- 
tion is considered to be the most reliable method avail- 
able since it uses Sutherland’s Eq. (3) for the variation 
of viscosity with temperature. 

It is apparent from Fig. 6 that the assumption u* = 
7* is a first approximation. The assumption »* = 
7 *°76 is seen to be good at flight temperatures (for 
which it was proposed) but unsatisfactory at the low 
temperatures of tunnels. On the 
other hand, the assumption u* = CT* yields good re- 
sults for both free-stream temperatures at Mach Num- 


supersonic wind 


bers up to 5. From these curves in Fig. 6, it can be 
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seen that for an insulated surface the skin friction is 
predicted by the present theory to within about 1 per 
cent for free-stream temperatures encountered in high 
Mach Number supersonic wind tunnels and to within 
about 5 per cent for extreme free-flight conditions up 
to Mach Numbers near 5. For noninsulated surfaces 
heated or cooled within reasonable practical limitations, 
the order of accuracy of the present theory is about the 
same as that illustrated by Fig. 6. 

Heat Transfer.—For the special case of constant sur- 
face temperature, the relationship between heat trans- 
fer and skin friction is 


QO/Re(Tw — TV Re = (Cap Re/2)Pr** (63) 


as was previously derived. The calculations of Crocco 
show that, even if Sutherland’s equation is used, Eq. 
63) is still valid to within about 0.5 per cent for any 
free-stream temperature. Therefore, the order of ac- 
curacy with which the present theory predicts the skin 
friction, as illustrated by Fig. 6, is also the order of ac- 
curacy with which it predicts the rate of heat transfer. 


General Effects of Variable Surface Temperature 


Heat-Transfer Coefficient.—Probably the most impor- 
tant single effect of variable surface temperature appears 
in the heat-transfer coefficient and, hence, Nusselt 
Number. Thus, the conventional definition of heat- 


transfer coefficient by means of the Newtonian cooling 


law 
q = k(T,, — T, 
vields 
_{T.*0, 0) — T.*1Ce he Ci 
sii ite ~ 7,° 2 * 


>, a,x* *Y,’(0) 


| » 

| l “ l 
Vow T.- 
Unfortunately, at each 
w* is equal to the 


(52) 


T,* 


as was previously derived. 
point where the surface temperature 
equilibrium temperature 7,*, this expression will be- 
come both positively infinite and negatively infinite 
in all but a few exceptional cases. (The physical ex- 
planation for this is given later.) The exceptional 
flows where an infinity does not appear are those with 
constant surface temperature, wherein Eq. (52) re- 
duces to h = h, and those wherein the surface-tempera- 
ture distribution possesses the extraordinary property 


that the limit 


lim 


Tw* — Te* 


/ 

> a,x*" Y,'(0)/ (T.* — T,*) 

3 | 

is finite at every point where 7,, = 7,. Of course, this 
latter condition generally is not possible. Since the 
equations for all boundary-layer characteristics other 
than h(T, u, g, and Q) do not become infinite under any 
physically possible conditions, it is clear that the cus- 
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tomary definition of heat-transfer coefficient is funda- 
mentally inappropriate for flows with variable surface 
temperature. For variable surface temperatures it ap- 
pears more convenient to deal directly with the dimen- 
sionless group gx/k_7., rather than hx/k,. 

Total Heat Transferred.—As far as the total rate of 
heat transfer is concerned, an illustrative comparison 
can be made between the present theory that accounts 
for the effects of a variable surface temperature and a 
customary calculation that assumes that the heat- 
transfer coefficient is that determined for a constant 
surface temperature (4). In order to simplify this 
comparison, C and C, temporarily will be taken as 
unity, since their actual values are unimportant for the 
purpose at hand. The customary method gives g = 


h(T, — T,.), and 


* i - Tk. Yo'(0) 
Q= h(T, — T.)dx = —- ; 
J0 2 

Pa 5 aa 
\ = > A,X dx* = 
Ps 0 eo 0 
7 hate) U oL = Ge 
- \ yo 9 n+ (1/2) 


The present theory gives, from Eq. (55), when C = 
Cw = i. 


ee a 
oe <> “_._ ¥",()) 
7 v. > n + (1/2) 


Thus, the ratio of the total heat transferred as given 
by the two methods is approximately 


QV > Qn ” / : ; Gea 
ea 5 ——*_ y,'@)/ ¥'@ Pn 
Q BZn+ (1/2) ‘ / o'( UD 


(64) 


This ratio obviously is not unity except in special cases. 
In fact, for different surface-temperature distributions, 
this ratio can be either positive, negative, zero, or in- 
finite. The particular example discussed later on 
yields a negative value for this ratio. 

Flow Separation.—-The calculations for the velocity 
profile allow some useful qualitative information to be 
obtained with regard to the effect of a heated surface 
on flow separation. Eq. (57), when evaluated at the 


surface, yields 


(S") (*) = f"(0) U. ; U, 

Oy/o  \wle * Nv xO 4Ty* 

In theoretical calculations the separation point of a 
flow is taken as that point where (0u/Oy),, = 0. Hence, 
the above equation shows that, for a heated plate in 
the absence of a pressure gradient, the slope of the 
velocity profile at the surface will vary approximately 
as the reciprocal of the absolute temperature. Conse- 
quently, it would theoretically take an infinite surface 
temperature under these conditions to cause flow separa- 


tion. 
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Dimensionless distance along plate, x*= " 


Figure 7 —Surface femperoture distribution for 
example problem. 


If, on the other hand, the upstream part of a surface 
is heated and the downstream part is subjected to an 
increasing pressure (adverse pressure gradient), then 
a smaller increase in pressure will cause separation 
than would be necessary if no heat were added to the 
upstream part. Likewise, if an adverse pressure gradi- 
ent simultameously exists over a heated region of a 
surface, flow separation will occur more readily than if 
that same surface were not heated. In general, 
it can be concluded that the addition of heat can 
contribute to but cannot by itself cause flow separa- 
tion. 


DISCUSSION OF AN EXAMPLE 


The particular example to be considered in detail is 
designed to illustrate conveniently some of the more im- 
portant results of the theory. Thus, for simplicity, 
two-dimensional flow is considered with C,, and C set 
equal to unity. This means that the small indirect 
effect of Mach Number on C,, and C and, consequently, 
on the heat transfer is neglected. The other boundary- 
layer characteristics (u and 7°) are calculated for two 
Mach Numbers, J, = 0.5 and M, = 3.0, which are 
considered to be representative of subsonic and super- 
sonic flows, respectively. 

The distribution of surface temperature selected for 
the example is given by the quadratic 


t(x*) = T,* — T,* = T,*(0.25 — 0.83x* + 0.33x*?) 
so that 


ao = 0.257;,* 


a, = 


AERONAUTICAL 


SCIENCES—SEPTEMBER, 1949 
a2 = 0.337,* 
lL, = 0, na>2 


This distribution is plotted in Fig. 7 with the ratio of 
surface temperature to equilibrium temperature as the 
ordinate. This distribution is considered to exist at 
both M, = 0.5 and M, = 3.0, and, therefore, the 
actual level of surface temperature becomes a junction 
of M., through the explicit dependence of T, on Mach 
Number. The purpose of this figure is to indicate the 
relative change of surface temperature with respect to 
the overall effective temperature difference (7, — T,) 
between the surface and the free stream. It will be 
noted that a rather severe case has been analyzed, since 
the overall temperature variation in the plate is one- 
holf the equilibrium temperature—or looking at it in 
another light, it is twice the maximum difference be- 
tween 7, and T,. 


Heat-Transfer Coefficient 


By substituting the coefficients dy, a, and dy into 
Eqs. (52) and (53) and using the values of Y,,’(0) 
given in Table 1, the heat-transfer coefficient is obtained 
as a quadratic function of x*. The resulting effect of 
the surface-temperature distribution (Fig. 7) on the 
local heat-transfer coefficient is illustrated in Fig. 8a, 
where the Nusselt Number divided by the square root 
of local Reynolds Number is plotted (solid line) as a 
function of the dimensionless distance along the plate. 
It may be noted from Eq. (52) that the effect of Mach 
Number is eliminated when C = 1 and the ratio 7,,/T7; is 
independent of M/.. Fig. Sa, therefore, applies for 
both M. = 0.5 and \/.. = 3.0 in the example under 
consideration. The dotted line in Fig. 8a represents 
the local heat-transfer coefficient determined by the 
conventional method wherein the effect of the variable 
surface temperature is not included. The dotted line 


is represented by 
Nu __ hx _ _ _ Yo'(0) 
VRe, kV Re: 2 


Several interesting effects may be observed from this 
figure. In the region of decreasing temperature poten- 
tial (0 < x* < 0.35), the local heat-transfer coefficient 
is less than that for a constant surface temperature 
plate. The heat-transfer coefficient is reduced so 
markedly that it becomes zero (x* = 0.2) and then even 
becomes negative (0.22 < x* < 0.35), which means 
that heat flows into the plate although the overall 
temperature potential is in the direction that would 


= 0.296 


normally cause heat to flow from the plate to the free 
stream. Since the heat flows in a direction opposite to 
the local overall temperature potential, the heat-trans- 
fer coefficient / in Newton's law of cooling 


h q iF « — i 


is necessarily negative over part of the surface. 
anomalous behavior continues downstream to x 
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Figure 8—Effect of variable surface temperature on the 
heat transfer chorocteristics of the example problem. 


0.35, where h becomes first negatively infinite and then 
suddenly becomes positively infinite. As was stated 
earlier, this behavior is due to the inapprdpriateness 
of the conventional heat-transfer coefficient when ap- 
plied to flows with variable surface temperature. From 
the physical point of view, this inappropriateness exists 
because the definition of h at a given point involves only 
conditions at that particular point (7.e., the local value 
of T,, — T.); whereas the boundary-layer properties 
at a given point do not depend only on the local tem- 
perature potential but, instead, depend in a complicated 
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manner on the entire distribution of surface tempera- 
ture upstream of the point in question. (The distribu- 
tion of surface temperature downstream of a given 
point has only an extremely small effect on the bound- 
ary-layer flow upstream.) Thus, if the surface tem- 
perature at a fixed point were maintained at the 
equilibrium temperature 7, and if the distribution of 
surface temperature upstream of that point were varied 
widely, then the changes in the velocity and tempera- 
ture profiles that would take place upstream would 
propagate downstream, thereby causing the tempera- 
ture profile, and, hence, g at the point concerned also 
to vary widely. The local temperature potential 
(1, — T,), however, is always zero in such a process, so 
that / is necessarily both positively infinite and nega- 
tively infinite at isolated points where 7, = T,. 

In the region of increasing temperature potential 
(0.35 < x* < 1), the heat-transfer coefficient (Fig. 
Sa) is higher than that for a plate of constant surface 
temperature. This is attributed primarily to the 
infinity at x* = 0.35. Actually, all of the wide varia- 
tions illustrated by Fig. Sa are a consequence of the 
fundamental inappropriateness of / when applied to 
variable surface temperature flows. 


Heat-Transfer Rate 

The local rate of heat transfer (g) for the example 
under consideration is obtained by substituting the 
coefficients do, a, and a2 into Eq. (50) and using the 
values of },,’(0) given in Table 1. The result is plot- 
ted as a solid line in Fig. Sb where the dimensionless 
group @x/k. T.V Re, is used as the ordinate. Also 
shown in Fig. 8b is the relation 


@ 


mre: wt GE St x 
Re T,% Re, 2 0 T.* 

ne (T,,* — T.*) 
which is obtained by neglecting the effects of variable 
surface temperature on the heat-transfer coefficient. 
The difference in the two methods of calculation is less 
pronounced in terms of the local heat-transfer rates 
than in the local heat-transfer coefficients. This is due 
to the fact that in the example considered the overall 
temperature differences are a minimum where the coef- 
ficients differ by a maximum. 

By substituting the coefficients ao, a, and dz into 
Eq. (64), it is found that the ratio between the total 
heat transferred as given by the two methods is Q/Q = 
—1.9. This means that not only is there a difference 
in the absolute value of the total heat transferred from 
the plate but that the actual direction of net heat trans- 
fer is different. Thus, for the example under considera- 
tion, the conventional method gives a completely incor- 
rect picture of the total rate of heat transfer by the 
plate. 
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Figure 9.- Temperature distribution in boundary layer of example 
problem 


Temperature and Velocity Distributions 


The temperature distribution as a function of 7 is 
obtained by substituting the coefficients do, a, and a» 
into Eq. (36), using the curves for r(n) and Y,(n) given 
in Figs. 4 and 5. The value of y corresponding to each 
value of n is calculated from Eq. (42), using the curves 
for r(n) and Y,,(n) given in Figs. 4and 5. The resulting 
temperature profiles are shown in Fig. 9. Although 
the thermal boundary layer is considerably thicker at 
the higher Mach Number, there is more heat trans- 
ferred than in the low Mach Number case. Because 
the local rate of heat transfer can be expressed as 


— &) " ae ay 
q = 7a oy a — w+e oy - 


it is apparent that this phenomenon is caused by the 
larger thermal conductivity at the high temperatures of 
the supersonic flow and also by the higher equilibrium 
temperature. 

The boundary-layer velocity distributions for the ex- 
ample under consideration are indicated in Fig. 10 at 
three positions along the plate, x* = 0.05, 0.35, and 
1.0. These profiles are obtained in a simple manner 
[Eq. (13)] by first using the curve of f’(n) in Fig. 3 
and then using the same relation between yn and y as 
was used for the temperature profiles. Although the 
velocity profiles at MM. = 0.5 very nearly correspond 
to the Blasius profile for isothermal incompressible 
flow, it is seen that the thicker profiles at M, = 3.0 
depart considerably from each other and also from the 
Blasius profile. The smaller velocity gradients and, 
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hence, larger boundary-layer thickness at MM. = 3.0 
are attributable to the lower densities (higher tempera- 
tures). The decrease in density causes the streamlines 
to spread more for the supersonic flow than for the sub- 
sonic flow. Fig. 10 also shows that the cool portion at 
the rear of the plate tends to make the velocity distri- 
butions more convex, possibly indicating increased 
stability. The high temperatures at the leading edge 
produce an inflection point in the velocity profile at x* = 
0.05 (barely noticeable in the figures) and probably indi- 
cate decreased stability. 

It is stressed that the above statements on stability 
are of a qualitative nature and are included primarily 
to make the reader aware of the application of this solu- 
tion in the consideration of boundary-layer stability. 
It is known from reference 2 that the boundary-layer- 
density distribution is also important in stability con- 
siderations. The density distribution (see temperature 
distributions) is more markedly affected by the surface- 
temperature variations than the velocity distribution 
The maintenance of stable boundary layers by adjust- 
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Figure /0.- Velocity distribution in boundary 
layer of example problem. 
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TEMPERATURE AND 
ing surface-temperature distributions, therefore, may 
have practical possibilities. This solution provides a 
basis for such studies but they are omitted here be- 
cause they are beyond the intended scope of the present 


paper. 
APPENDIX A 


Maximum Surface-Temperature Gradient Consistent with 

Boundary-Layer Theory 

Since the effects of relatively large temperature gradi- 
ents along the direction parallel to the wall are to be 
included in the analysis, it is necessary to re-examine 
the boundary-layer momentum and energy relation- 
ships in order to see what restrictions or what additional 
terms, if any, are needed in the usual boundary-layer 
equations. The Navier-Stokes equations expressing 
conservation of momentum for the steady flow of a 
viscous compressible fluid are (using the summation 
convention for conciseness, wherein a repeated index in 
, term designates summation): 


Ou; Op ‘ °| (oe 4 ) 
Mi ex, | Oxy NOx, | Ox; 
2 ° ( ot) 
uM (Al) 
3 Ox:\ Ox; 


lhe left side of this equation represents the mass times 
acceleration of a unit volume of fluid, and the right side 
represents the unbalanced pressure and viscous force 


per unit volume. Thus, for the case of two-dimensional 


flowi=1(m=%,wm=y, m% U, Ue v), the above 
equation reduces to 

Ou | Ou Op 4 O ( ; =} 4 
ou + pu = -— a) 

Ox oy Ox Ox Ox 


°t ou ‘ *) | _2 ot (> 4 | (A2) 
oy us y ow 3 oxL Vax * dy a 


rhroughout the following analysis the pressure will 
be taken as constant along the x direction, which is 
parallel to the surface (Fig. 1). By tracing along a 
streamline such as shown in Fig. 1, and repeating the 
method given by Prandtl for estimating the order of 
magnitude of the various terms in Eq. (A2), it is seen 
that 
Ou/Ox ~ U/L, Ou/oy ~ U ,/6 
Moreover, the continuity equation 


(0/Ox)(pu) + (0/Oy)(pv) = O (A3) 


allows the remaining terms in Eq. (A2) to be estimated 
is follows: 


Ov/Oy ~ Ou/Ox ~ U/L 


6 
v~ fo ox)dy ~ (U,/L)é 
0 


Ov/Ox ~ v/L ~ (U,/L)(6/L) 
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Terms in Eq. (A2) such as 0u/Ox can be neglected 
compared to Ou/Oy, since the foregoing estimations 
show that this amounts to neglecting terms of order 
6/L compared to unity. This, of course, is valid for 
Reynolds Numbers of the order encountered in aerody- 
namics since 


6/L = 6/V Re 


Neglecting the terms in Eq. (1), which are of the order 
5/L compared to unity, gives the familiar momentum 
equations for laminar boundary-layer flow 

Ou Ou o ( ot) 


s Pa =: 
ao TP dy ody r ov 


i= 2; odp/oy = 0 


| 


(A4) 


It is now necessary to repeat the above analysis using 
the energy equation, since it is this latter equation that 
explicitly involves temperature gradients in the x direc- 
tion. The differential equation giving the balance of 
energy for steady flow of a viscous compressible fluid is 


qT) _ | , 2 ( = ‘ 
ti On, Ox, * Da\ Oxy 


Ou; / OU; Ou; 2/0u;\? . 
M + a & (A5) 
Ox; \Ox,; Ox; 3\Ox; 
It is to be noted that this is a single equation, since> 
energy is a scalar quantity, whereas the momentum 
Eqs. (1) actually represents three separate equations, 
The left 
side of Eq. (A5) represents the rate of increase in en- 


one for each component of the momentum. 


thalpy per unit volume of fluid, and the three terms 
on the right side represent the rate of work done by 
pressure forces, the rate of heat transfer by conduction, 
and the rate that heat is dissipated by viscosity, re- 
spectively. Eq. (A5) written out in full for two space 
dimensions becomes 
O(c,T) O(c,l°) Op Op 

id Ox “a Oy = Ox "7 Oy 


3) (: —) 4 a) (« ~) 4 of a(S) 4 
Ox\ Ox oy\ Oy Ox 
Ou [Ou * Ov/Ov . Ou 
= “i a + =) | r 
nf 2(5°) . (= 4 =) | (A6) 
oy 3\0x Ov 
With constant pressure and with those velocity terms 


neglected which have already been shown to be of 
order 5/L compared to unity, Eq. (A6) reduces to 


O(c,T) O(c,T) ) ( -) 3) ( ) 
v - = k k 
- Ox “a oy Ox\ Ox i Oy\ Oy * 


(S*) (A7 
rm a Av) 


This differs from the usual energy equation for bound- 
ary-layer flow only by the presence of the first term on 
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the right-hand side. Considering order of magnitudes 


again, 


(0/Ox)[k(OT/Ox)] ~ k(OT/Ox) L 
: a T.-T 
a eee Pe. 
oy\ Oy (6,)? 


0: 
For gases such as air, where the Prandtl Number is 
nearly unity, the thickness (4,) of the thermal boundary 
layer is of the same order of magnitude as the thickness 
5 of the velocity boundary layer. Hence, 


(0/Ox) [k(OT/ Ox) ] Al (OT /Ox) » | tain 

(0/dy)[R(OT/Oy)}] LL(Ty, — T,) 6 poe 

For applications wherein the temperature gradient 

OT /Ox),, along the surface is of the same order or less 

than the gradient (7,, — T.,) ‘6 normal to the surface— 
that is, for applications wherein 

(07 /0x)_e ~ (Ty — T. (A9) 


)/6 or less 
Eq. (A8) becomes 


(0/Ox)[k(OT/Ox)] 6 


(d/dy)[R(OT/dy)] — L 


Terms of order 6/Z compared to unity, however, have 
already been consistently neglected. Hence, within 
the order of accuracy of the momentum Eqs. (A4), the 
first term on the right-hand side of Eq. (A7) may be 
neglected, provided, of course, that Eq. (A) is satis- 
fied. 

The results of the foregoing considerations on order 
of magnitudes may be summarized by stating that the 
usual boundary-layer equations 


Ou Ou a) Ou 
pu — + pv i = > (A st momentum 
oo . oO we 
a (pu) + dy (pv) = 0 continuity | 
(1) 
pu tees pv —_ i mG =) + 
Ou\? 
(s) energy 


are applicable to flows with large temperature gradients 
in the x direction, provided 
(2) 


(OT /0x)» ~ (Ty — T.)/6 or less 


APPENDIX B 


Numerical Solution of the Functions Y ,(n) 


The temperature distribution functions }',(n), are 
determined by integrating Eq. (33), 


VY,” + PrfY,’ — 2Pmf'Y, 


Il 


(33) 
with the boundary conditions 


Y,(0) = 1 


OF THE AERONAUTICAL 


SCIENCES—SEPTEMBER, 1949 


Y,(o) =0 


Because the Blasius functions f(n) and f’(y) cannot be 
expressed in closed form, Eq. (33) is intractable by the 
ordinary analytical methods, and, hence, numerical 
integration is indicated. In order to start the usual 
methods of stepwise integration of a second-order dif. 
ferential equation, it is necessary to know the values of 
the dependent function and its first two derivatives at 
one of the boundaries. When the derivatives are not 
known at a boundary, it is necessary to estimate them 
and then check the validity by the degree to which the 
other boundary conditions are satisfied. Such a process 
usually reduces to one of trial and error, necessitating 
many tedious computations. Fortunately, this process 
can be avoided in the solution of Eq. (33) and similar 
boundary-layer equations where the asymptotic solu- 
tion can be found. It will be noticed from Fig. 3 that 
the functions f(n) and f’(n) become linear and constant, 
respectively, for large values of 7». Expressions for 
these functions accurate to four significant figures in 


the range 7 > 4.1 are:!® 


f'(n) = 


9 
f(n) = 2(n — 0.86038) 


If these equations are introduced into Eq. (33) there 


results the asymptotic equation for Y,, 


Y,” + 2Pr(yn — 0.86038) Y,’ — 4nPrY, = 0 (BI 


letting x = V Pr(y — 0.86038), Eq. (B1) becomes* 


Y,” + 2xY,’ — 4nY, = 0 (B2 
The general solution of the linear differential Eq. (B2 
can be determined, if one particular solution (U’) is 
known, by means of the expression 


. 


. ‘e—* dx c 
yY = uo(c f 7 + C2) B3 


where C, and C; are arbitrary constants. A particular 
solution U of Eq. (B2) is obtained readily when it is 
recognized that Eq. (B2) can be transformed to Her- 
mite’s differential equation. Letting ix = £ and m = 
2n, Eq. (B2) is transformed to the Hermite equation 


V,” — 2EY,,’ + 2mY,, = 0 B4 
for which the following Hermite polynomial is a partic- 
ular solution 

= m(m — 1) 


ym (9t\"—2 


1! 


m(m — 1)(m — 2)(m — 3) 


H,(§) = 2 


Str 


(9e"—" — ... (BS 
1 s/ 


<> 


Since m(= 2n) is always an even number in the present 
application, the polynomial is always real and cari be 





* The use of x as a dummy independent variable throughout 
the remainder of this appendix is not to be confused with th 


space coordinate x. 
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written 


2n(2n oll 1 ) 


= H.,(ix) = (—1)"(2x)” — 


1! 
2n(2n — 1)(2n — 2)(2n — 3) 
~1)"-1(2,)2"-? ls i. 9) 
| 2! 
(—1)*-*(2¢)"—-* — (B6) 
Substituting Eq. (B6) into Eq. (B3), there results 
*e-*" dx 
Y, = He,(tx) af te 2 c| (B7) 
‘ | ‘Jo H2,,?(1x) ; 
Imposing the boundary condition Y, = 0 asy— ©, 
it is found that 
C C f. — BS) 
9 = _ ( ‘ 
‘ ‘J. Ho,2(ix) , 
therefore, 
4 CH (i f th (B9) 
n = —OjlTen\tx, comes ( 
nS, 


is the asymptotic solution to Eq. (33). The first two 

derivatives needed in order to start the numerical inte- 

eration are obtained from Eqs. (B9) and (B2), 

my , aa n° @-* dx " el % 

‘eel | a) f Ho,2(ix) 7 aa | > 
He,'(ix) .. me 
HiAés) ° * Haaléz) 


Y,” = —2xY,' + 4nJ/, 


(B10) 


(B11) 


Eqs. (B9), (B10), and (B11) were evaluated numeri- 
cally for values of x corresponding to 7 = 4.1 and y = 
12. Once the functions Y,, Y,’, and Y,” were ob- 
tained at these values of , the stepwise numerical inte- 
gration was started toward » = 0. The numerical 
integration method of Runge-Kutta was used in the 
calculations. The value of Y,(0) obtained at the end 
of a numerical integration process was never unity 
the necessary boundary condition), since the correct 
C,; was not known beforehand. This does not cause 
complications because Eq. (33) is homogeneous and the 
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function obtained in the stepwise process can therefore 
be multiplied by that constant factor that makes 
Y,(0) = 1. 
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Readers’ Forum 





RIEF REPORTS of investigations in the aeronautical sciences 
and discussions of papers published in the JOURNAL will be 
presented in this special department. The publication will be com- 
pleted 6 to 8 weeks after receipt of the material. No proofs will be 
sent to the authors. The Editorial Committee does not hold itself 
responsible for the opinions expressed by the correspondents. Con 
tributions should not exceed 800 words in length 


Addition of Heat to a Gas Flowing in a Pipe 
at Sonic and Subsonic Speeds 


Maurice Roy 
Professor, Ecole Polytechnique, Paris, France 
April, 1949 


- THEIR PAPER ‘‘Addition of Heat to a Gas Flowing in a Pipe 
at Subsonic Speed”’ (JouURNAL OF THE AERONAUTICAL SCI- 
ENCES, Vol. 16, No. 2, p. 84, February, 1949), J. V. Foa and G 
Rudinger contest, together with the opinions of many other au 
thors, the one I expressed* when I stated that a steady-flow proc- 
‘ss at sonic speed is inconsistent with an external exchange of 
heat or with a change in the chemical state of the flow 

While developing an interesting conception, the authors, as a 
matter of fact, investigate a case quite different from the case of 
the steady flow that alone was considered by me. They con 
template, in fact, the transient conditions brought about by 
supplying heat locally to the initial flow, and then they discuss 
the final steady flow to which the disturbances thus produced 
would lead. But this conception of the final flow could not es- 
cape the inferences (and, in particular, the one above mentioned) 
that are directly drawn from the study of the final steady flow. 

How is it then that there is a disagreement? 

Concerning the physical hypotheses about the flow, I have 
issumed that the fluid was governed by the laws of perfect gases 
Cp and y constant) and that its viscosity and conductivity were 
negligible, assumptions that Foa and Rudinger seem also to admit 

In that case, the divergence between our results could be due 
to the method used in determining the permanent final flow 
Perhaps, too, this divergence is at least partly to be explained 
by the fact that the authors localize the supply of heat in a cross 
section of the pipe—a conception difficult to realize. 

As for me, I have assumed this supply to be permanent and 
distributed in the mass of the flow with an intensity per unit of 
mass and per unit of time varying along the pipe. It must be 
noted that a one-dimensional flow and an external supply of heat 
do not easily fit in with the supposed absence of internal con- 
ductivity. Moreover, as I have emphasized, this lack of con- 


ductivity excludes the so-called ‘‘deflagration”’ type of flame 
propagation. 

I think today that my theory, published in 1946 when I did 
not know of other concordant studies, is far too schematic and 
should be completed, taking into account the internal conduc 
tivity, the turbulence, and also, of course, the viscosity of the 
fluid. 

In addition, in the case of chemical reactions, more accurate 
physical data ought to be known concerning the speeds of the 
reactions. 


* Reference 13 of the authors and previously a note ‘‘Théorie unidimen 
sionnelle des tuyéres-foyers’’ (Comptes rendus de l’'Académie des Sciences, 
Paris, Vol. 222, p. 835, April 8, 1946 
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For the sake of simplification, when considering the case of 
sufficiently large pipes, the assumptions of the flow being steady 
and one-dimensional could still be retained. , 

However, even so simplified, the problem raised remains q 
difficult one and calls for new studies. 


Heat Addition to a Flowing Gas 


J. V. Foa and G. Rudinger 
Cornell Aeronautical Laboratory, Inc., Buffalo, N.Y. 
June 21, 1949 

UR PAPER ON GAS FLOWS with heat addition! has been the 
O subject of several comments?~‘ in the Readers’ Forum 
Some of these comments seem to imply that the disagreement be. 
tween the conclusions of our paper and those of previous steady. 
state analyses is only apparent and essentially reducible to a 
matter of semantics. This implication is clearly expressed in 
Hall’s comment,? where it is pointed out that we have used the 
phrase “initial condition” to describe the state in time preceding 
the addition of heat anywhere in the system, whereas other inves- 
tigators have found it convenient to use it as referring to the 
state in space preceding heat addition. 

One is, of course, at liberty to deal with a steady state where 
heat is being added at a constant rate and to refer to the upstream 
conditions as “‘initial’’ conditions. This is, as Hall points out, 
a widely established use of the term and indeed a valid one if it is 
adhered to consistently from the formulation of the propositions 
to that of the final inferences. When, however, the results of 
such a steady-state analysis are interpreted to mean that the heat 
transfer coefficient becomes zero, or that combustion stops wher 
the gas reaches sonic velocity, or that the upstream condition 
must be modified when the amount of heat added exceeds a cer 
tain “limiting value,’ then it is clear that in the final inference 
the word ‘“‘initial’’ has come to refer to the state preceding the 
addition of heat in time at a fixed position in space and that the 
rules of valid inference have thus been violated. 

The 
should, as we stated in our Introduction, ! 


equations derived in previous steady-state analyses 
“only be interprete 
to relate the conditions upstream and downstream of the heating 
region.’ These relations [which are restated in our paper 
Eqs. (8), (9), and (10)] are well-established and have been plotted 
in various forms (such as the one presented by Barry and Ful 
ton*). The three steady-state relations, together with thre 
boundary conditions, are sufficient to determine all flow and stat 
parameters in a steady régime of heat addition when the heating 
rate is known. Obviously, any change of this rate will result 1 
a change of flow and/or state parameters, upstream and/or 
downstream of the heating region, and therefore in the establis! 
ment of a new steady state 

Of the three cases mentioned by Barry and Fulton,* the firs 


We < 1.0, and 


were discussed 1! 


two—(a) constant m, 7\, and also p. when 
constant 7, 7;, and also po when JM, < 1.0 
our paper as illustrations of a treatment that is completely gener 
within the framework of the assumptions that are listed there 
The third case (constant JJ, and 7}) is not sufficiently define 
and could not be applied to a physical system without specifica 
tion of the missing boundary condition. Thus—e.g., in the 
systems considered in detail in our analysis (constant-area ducts 
discharging directly into the surrounding atmosphere)—the third 
boundary condition is p2 = const. for Mz < 1.0, or M2 = 1.0, ane 
it can then be verified that when heat is added at M > 1/¥V1 
the downstream static temperature does indeed decrease, but 
the amount of heat that must be removed from the air supply 


system in order to maintain the prescribed upstream boundary 
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conditions is larger than that of the heat that is added to the 
fow. This mechanism is one of no interest and was for this rea- 
son specifically excluded in our study where it was concluded 
page 90 that “except in the irrelevant case where the addition 
of heat is more than offset by a simultaneous heat removal, the 
addition of heat will always result in an increase of the static 
temperature of the gas above its initial value, regardless of the 
initial Mach Number.” 


More generally, it is true that, if the upstream conditions 
could be forced to remain constant when the rate of heat addition 
is changed, then the downstream static temperature would de- 
crease when the heating rate is increased at M > 1/Vy. It 
must be noted, however, that in order to achieve such a situation, 
the control of the upstream conditions could not be effected di- 
rectly but only through a special adjustment of the downstream 
boundary conditions, so that the latter would become functions 
of the amount of heat added. Since the boundary conditions 
are essential characteristics of any physical system, it is clear 
that any dependence of the boundary conditions on the heating 
rate would of necessity imply a continuous modification of the 
physical system itself as the heating rate is changed——a situation 
that is distinct from that treated in the problem under dis- 
boundary 


cussion. A with automatically controlled 


conditions could of course be made to behave in any desired 


system 


manner. 


Roy’s statement? that a ‘‘steady flow process at sonic speed is 
inconsistent with an external exchange of heat or with a change of 
the chemical state of the flow’”’ is correct and in full agreement 
with our own conclusions. On the other hand, his recommenda- 
tion’ that the Mach Number be kept as low as possible in the 
region of ignition because practically no combustion could take 
place in a mixture that is ignited at M > 0.6 reflects an erroneous 
interpretation of his results. Similarly, Barry and Fulton’s re- 
cent discussion,® where they properly recognize the importance of 
the boundary conditions, can hardly be reconciled with Barry’s 
previous statement® that ‘‘the amount of fuel which may be 
burned is limited by choking.’’ In this connection, it may be 
pertinent to recall the cautionary note that was contained in our 
paper! (page 85), where it was explained that ‘‘when the inlet 
boundary conditions are assumed to be fixed, special care must be 
exercised in drawing conclusions from the results. For instance, 
the imaginary solutions obtained in certain cases must be under- 
stood to imply the physical impossibility of the assumed condi- 
tions rather than a limitation to the amount of heat that can be 
added.”’ 


Hall? and Roy‘ have also attempted to explain the disagree- 
ment between their own and our conclusions on the basis of argu- 
ments (heat transmission within the gas, localization of the heat 
source, consideration of transient phenomena, effects of viscosity 
and turbulence, mechanism of the chemical reactions) which 
Actually, this 
disagreement does not arise from the manner in which the rela- 


are completely foreign to the problem as treated. 


tions governing flows with heat addition are derived and treated 
but only from their physical interpretation. 
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FORUM 


A Discussion of the Theories of Plasticity 
D. C. Drucker 


Brown University, Providence, R. |. 
March 1, 1949 


gens EXCELLENT INDEPENDENT PAPERS by Dr. Batdorf (July 
JouRNAL) and Professor Bijlaard (this issue) deal with the 
same subject and come to essentially the Same conclusion—a de- 
formation theory of plasticity predicts the correct plastic buck- 
ling load for plates; a particular incremental, or flow, theory 
does not. Their conclusions are substantiated in all cases in- 
vestigated experimentally. 
ported by Dr. Batdorf are just as conclusively against deformation 
and in favor of an incremental theory of plasticity. Dr. Batdorf 
therefore states that neither type of theory is proper. Professor 
Bijlaard adopts the point of view that the answer justifies the 
means. 

Before analyzing the main issue, it seems desirable to attempt 
a clarification of some of the basic ideas and terminology of plas- 
ticity theory included in the papers. 
wording used to indicate that a relation exists between finite 


However, crucial direct tests as re- 


Deformation theory is the 


changes in stress and strain, usually between stress and strain 
Incremental, or flow, theory is used to emphasize 
The 
stress-strain relation is in terms of the differentials of strain 


themselves. 
the necessity for considering the history or path of loading 


Total strain or a finite strain increment is strongly dependent on 
the manner in which the final state of stress is reached. In all 
the forms written down to date, the two types of theory coincide 
only when all components of stress are increased in ratio, oy = 
Koi;. 

Deformation theory can be shown to be mathematically incon- 
sistent in general and therefore, of necessity, often physically 
impossible. As theories either do or do not involve the history of 
loading, they must be of one type or the other and therefore must 
be incremental to be valid. 

Some confusion seems to have arisen in this respect because of 
both authors’ identification of incremental theory with its simplest 
possible form which they employ. In tensor notation, to save 
space, 


dey = dey’ + dei” = deij® + FSydf 


where de, the total strain increment, is composed of an elastic 
increment, de;,°, and a plastic dejj’. 
further plastic deformation, f, is Jz = 1/2SjSi or equally well the 
octahedral shearing stress. F = F(f), and Sj is the stress devia 
tion, Siy = aij — 1/soxx5y, the stress with the average normal 
stress taken away. Such a theory is the easiest to handle of all 
incremental work-hardening types, but it is of the class termed 
isotropic stress hardening. As it is isotropic, it implicitly assumes 
that if the yield point is raised by tension in one direction it is 
raised equally for tension in any other direction. Much worse, it 
assumes a rise in the compressive yield point of the same amount 
as for the tensile and so is by no means in good agreement with 


The criterion of loading or 


experiment. 
Far more general incremental theories have been formulated, 
the most general at present may be written in the form 


of ( of 
dey = dey’ + G—(— 
— 7 + {= 


doy 

where the criterion of loading, f, may be as anisotropic as desired 
and may depend upon stress, strain, and their history, as may G. 
Such dependence obviously makes for tremendous mathematical 
complication, but simple experiments such as alternating plastic 
tension and compression show beyond doubt that it does exist. 
It is likely that simple forms will be adequate for most problems. 

In particular, the choice of 


f= 


where SS is the stress deviation, ¢” is the plastic strain, and m is a 


(Sz — mez”)? + (Sy — mey”)? + (Sz; — me”)? +... 
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constant, gives both a Bauschinger effect and a lower yield point 
for tension in the y-direction than in the x-direction, following 
tension in the x-direction. The major objection of Dr. Batdorf 
to incremental theory is thus seen to be an objection to the par- 
ticular isotropic stress theory f = J2. Nevertheless, the theory 
proposed by Dr. Batdorf apparently does differ basically from 
the G(Of/deiy) [Of /Oexde,;) form, which assumes a linear relation 
between differentials of stress and of strain. Such an assumption 
is by no means a necessary one, but the great difficulties now 
encountered indicate that the possibility of ever solving stress 
analysis problems in plasticity with a nonlinear theory is certainly 
remote. 

The assumption of linearity can be given direct physical mean- 
ing. In the case of shear imposed on an initial tension or com- 
pression which is held fixed, line A on Dr. Batdorf’s Fig. 3, incre- 
mental theory predicts, and experiments confirm, that the first 
increment of shear strain is purely elastic. Unlike Dr. Batdorf’s 
theory, linearity further assumes that this will be true when the 
normal stress is increased, line B of Dr. Batdorf’s Fig. 3. This is 
an extremely important point, and tests will soon be under way 
at Brown University (Contract N7onr-358, Office of Naval Re- 
search and Bureau of Ships) to check or disprove linearity by 
direct measurement. 

Now to return to the basic question. 
that deformation theory (which is invalid in general) gives correct 
plastic buckling loads while incremental, or flow, theory (which 
agrees with many physical facts ) does not? In the first place, 
the comparison that has been made is not entirely fair. The 
deformation theory results are based on the Shanley idea of 
buckling proceeding as the load increases and so should corre- 
spond to the start of lateral displacement and not quite the meas 
ured buckling load. Secondly, if the Shanley effect is taken into 
account, answers obtained with incremental theory for initially 
perfectly flat rectangular plates supported on all sides become 
reasonable. Bending of plates with a slight initial bow would 
lead to similar results. 

On the other hand, the Shanley effect alone is not sufficient to 
account for the present discrepancy in the start of twisting of a 
long cruciform section or rectangular plate hinged on one edge 
parallel to the compression direction and free on the other. This 
is the case of shear added to compression, and, as described pre- 
viously, according to linear theory the initial modulus in shear is 
the elastic modulus whether the compression increases or remains 


Why do the authors find 


constant. 

Two possibilities exist. 
in a sense corresponds to Professor Bijlaard’s suggestion that 
deformation theory may be admissible for the small strain incre 
ments considered in buckling, as well as to Dr. Batdorf’s theory. 
If not, another factor must be taken into account in the unex- 
plained case. For example, an extremely small initial deflection 
or twist will greatly lower the effective shear modulus under in- 
creasing axial load and may well provide the explanation. 

Both possibilities should be explored, but even if linear incre- 
mental theory is found to be an oversimplification and if deforma- 
tion theory should prove adequate for buckling, the conclusion 
cannot, and should not, be extended to general plasticity prob- 
Direct experimentation only can establish the proper 


Linear theory tay be incorrect, which 


lems. 
incremental theory. 


On Nonsteady Motion of Delta Wings 


John W. Miles 
Department of Engineering, University of California at Los Angeles 


May 27, 1949 


W: CONSIDER HERE the nonsteady motion of a pointed wing 
in a supersonic stream directed along the positive x axis. 
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Let the port and starboard edges be yp(x) and ys(x), respectively, 


where 
yp(O = y;(0) = 0) 
ys'(x) > tan 8 . 
¥p'(x) < — tane | 
6 = sin~! (¢/U (2 


which is to say that the wing is a ‘‘simple’”’ plan form and the 
problem is “purely supersonic.’’ 6 is the Mach angle, U is the 
free-stream velocity, and cis the velocity of sound. For ap 
isosceles, delta wing y,(x) = x tan u, Vy(x) = —x tan wu, uw >6 
We are interested primarily in two problems, corresponding to 
the response of this wing to downwash distributions exhibiting 
harmonic and step function time dependences, respectively. 

If the prescribed downwash velocity normal to the plane of the 
O+ 


wing (z = is w(x, y) exp(iwt), positive down, the resulting 


(linearized) velocity potential at the wing is given by! 





d(x, vy, 0+) = dt n) w’ (§, n) tané 
0 
1 9 9 — 
g(x, y) = — (x* tan* 6 — y* exp(— Aa x 
us 
cos [\ sin 0(x? tan? 6 — y?)/?]-l(x tan @ — |y 
\ = (wb/U) sec? 0 4 


where b is a characteristic length. Unfortunately, the integra- 


tions demanded by Eq. (3) have successfully resisted reductior 
to known (i.e., tabulated) functions, even for the simplest w(x, y 
However, it is possible to reduce spanwise, weighted integrals of 
this potential to integrals arising in the two-dimensional prob- 


lem, Thus, we introduce 


¥(x) = o(x, y, OF) f(y)dy o) 


yp(x) 
where f(y) is an appropriate weighting function. Introducing ¢ 
from Eq. (3) and making certain transformations, we obtain th 


result 
x “ye lé 
v(x = tanég dé dn GC en E, n)w £, ‘ if 
0 S yplé 
: ' 1 ; ; 
G(x, vy) = exp(— 2dAx)1(x cos (Ax sin @ sin ¢) 
r Jo 
veos¢ + yid¢g (7 
The (dimensionless) weighted force on spanwise strip is given by 
™ 1 ys (x 
pl = 
y(x) = ttan @| - p(x, y, 0-) — 
9 7 
= 4 - Vz 
p(x, v, OF) flaid 8 
m) w\. 
= i y(¥) cot é Q 
Ox [ 
The most important cases in practice are fo(y) = land fi(y) =) 


corresponding to the calculation of the force and rolling moment 
respectively, on a spanwise strip. Making these substitutions 


in Eq. (7), we find that Eq. (6) may be rewritten 


. 


x 
tan @ / G — f)w,(E)dé, n = 0, 1 4 
e 0 


Xx sin @ l 


IL - = 


¥xaA\X) = 


Go(x) = exp(—idAx 


W(X) = w(x, y) y"dy 
vplx 

Since Gy(x) is essentially the two-dimensional Green’s function, 
the calculation of the lift and moments on our wing is reduced to 
an equivalent two-dimensional problem. Similarly, the problem 
of finding the wing loads due to the deflection of a delta control 
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surface with inboard end y = a may be reduced to an equivalent 
two-dimensional problem, it being necessary only to assume 
w(x, y) = w(x)1l(y — a) in Eq. (11). In essence, Eq. (9) states 
that “‘strip’’ theory is exact for m = 0, 1. 

In the case of step function loading, we assume the downwash 
in the form w(x, y)1[t¢ — fo(x)]. The result is expressed as in 
Eq. (6), but G(x, y) is replaced by its Bromwich integral, viz., 


x ya(€) 
y(x, t) = tan 6 dt dnA[x — £&, 7, t — to(t)] w(é, n) 
0 yp (§) 
(12) 
1 di (iwt)G 
> exp(2 x, y) 
A(x, %, j=. — lim SS: (13a) 
27 0 — e+ tw 
1 Tr 
= — def(x cos ¢ + y){1[ Ut — x sec? @ X 
2r 0 


1 +sin @ sin g)] + 1 [Ut — x sec? 0(1 — sin @sin g)]}  (13b) 





For the cases fo(y) = landf,(y) = y, we find 
hee ; ; 
J = tan @ Ao _——____— wn(é)di, n = 0,1 14) 
0 Ct — Uto(t 
w) = 1[(1 — sin @) — w] + 
I . . , ‘ R w — cos?@ 
l[w — (1 — sin @)J1[(1 + sin@) — w] cos~' | ————— 
- w sin @ 
(15) 
For a uniform gust entry (at ¢ = 0), Uto(x) = x, while for a sud- 


den change in angle of attack (at ¢ = 0), &(x) =0. Either sym- 

metric or asymmetric gusts may be considered, and we can also 

find the response of the wing to sudden symmetric or asymmetric 
deflections of the control surfaces. 

As an interesting special case, we remark that if we take \ = 

n = Oin Eqs. (9)-—(11), corresponding to steady symmetric flow, 
the lift coefficient for the wing is 

Cr = 4 tan 6[w(x, y) Jave. (16) 


a result previously obtained by Professor Lagerstrom of GAL- 
CIT.® 


REFERENCES 

! Garrick, I. E., and Rubinow, S. I., Theoretical Study of Air Forces on an 
Oscillating or Steady Thin Wing in a Supersonic Main Stream, N.A.C.A. 
T.N. No. 1383, 1947. 

? Miles, J. W., The Aerodynamic Forces on an Oscillating Airfoil at Super- 

nic Speeds, Journal of the Aeronautical Sciences, Vol. 14, No. 10, pp. 351- 
358, October, 1947. 

* Verbal communication. 


On Transient Two-Dimensional Flows at 
Supersonic Speeds 


Charles E. Watkins 

Research Scientist, Langley Aeronautical Laboratory, N.A.C.A.. 
Langley Air Force Base, Va. 

May 31, 1949 


N THE May, 1949, issue of the JOURNAL OF AERONAUTICAL 

SCIENCES, Biot! treats the chordwise lift distribution on a 
two-dimensional supersonic wing penetrating a sharp-edged gust 
by two methods: one based on a surface distribution of sources 
in the fluid at rest and the other by considerations of conical 
flow. In the present note it is shown, with the aid of the Laplace 
transformation, that this lift distribution, as well as that due toa 
two-dimensional wing undergoing a sudden change in angle of 
ittack, follows directly from the moving source solution given 
by Garrick and Rubinow.? This approach is of particular inter- 
est in that the different phases of the time histories of the lift 
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Shape functions in terms of boundary conditions 
for sudden change in angle of attack. 


distributions automatically present themselves as mathematical 
limitations in the inversions of the Laplace transformation. 
Furthermore, more general transient problems follow directly. 
Before proceeding with the analysis it may be pointed out that 
Heaslet and Lomax have used the source distribution method 
employed by Biot and independently obtained the lift distribu- 


tion for both gust problems. 


ANALYSIS 


Consider the velocity potential due to a distribution of moving 
sources which satisfies the boundary conditions for a thin oscil- 
lating wing in supersonic flow referred to a uniformly moving 
coordinate system as given in? 

1 


p(x,s,t) = — i 
28r Jo JO 


&, O)w(t)d 8 dé (1) 


where W(é, 6) and w(t) = [W(t — 71) + W(t — 72)] are, respec- 
tively, spatial and time functions of airfoil shape to be chosen 
so as to satisfy given boundary conditions. The variables x, 2z, 
£, 0 refer to rectangular coordinates, ¢ is time, 7; and r2 
lag functions associated with the physical mechanism of a moving 


ire time- 


supersonic source and defined by 


7} o 2 S(u- sin@), mm = ‘— 5 (M+ sine 
CB? cs? 

M is Mach Number, c is velocity of sound, and 8 = ¥ ‘M2? — 1. 
Eq. (1) is by no means limited to an oscillating airfoil. In other 
words, the time function need not be sinusoidal but may be chosen 
in such forms that Eq. (1) will represent the velocity potential 
for a two-dimensional wing in other types of time-varying super- 
sonic flow. In particular, the velocity potential for a two-di- 
mensional wing undergoing a sudden change in angle of attack 
or upon penetrating a sharp edged gust may be readily obtained 
from this equation: 

(1) Sudden Change in Angle of Attack—Consider first a wing 
undergoing a sudden change in angle of attack. The shape func- 
tion in terms of the boundary conditions for this case may be ex- 
pressed as follows (see Fig. 1) 

W(é, 9 — Va over the range occupied by the wing} (2) 
= ( outside the range occupied by the wing \ a 

where V is the velocity in free stream direction and a is the change 
in angle of attack. 

The corresponding time function w(t) is the unit step function; 


hence, 
[ z( — 71) t+ wt — rr ] = 0, when t<n 
= 1, whenr, < t < re 3) 
= 2, when t>n 


can be performed 
The Laplace trans- 


The integration with respect to @ in Eq. (1) 
with the aid of the Laplace transformation. 
form of Eq. (3) is 

Sti 4 e572) (4) 


w(s) = (1/s)(e7 
where s denotes the variable or operator in the Laplace trans- 
formation. Making use of Eqs. (2), (3), and (4) and the defini- 
tions of 7; and 72, the Laplace transformation of Eq. (1 


ered independent of z) is 


consid- 
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Ma ee 
Va é cB2 x- & i iC 
gis) = — cosh { s - sin @ }d@dt 5) 
Br Jo 0 s cB? 
Making use of the known Bessel function relation 
Bg 
(1/2) cosh (A sin 0)d@ = Io(d ti 
0 
the integration with respect to @ yields 
S. =i= , > : 
- va e ct x—-cé 
¥\(s) = TI s = dé 7 
B 0 § cp? 


For the integration with respect to &, it is convenient to revert to 
the inverse of the transform of Eq. (7), which may be found in 
tables of Laplace transforms.‘ Reverting to the direct function 


and performing the integration yields 
¢9(x, t) = 


va Mx = 
x cos a 


or x 


CB*t x — Met 
+ ¢Btl r — cos™! 
ct 


where the inverse trigonometric terms in Eq. (8) must be inter 
preted to be: O when the arguments are > 1, and 7 when the 
arguments are < — 1. 

The expression for the ratio of superstream pressure to dy- 


namic pressure is 


} l O¢ L 0¢ 
V of Ox 


— (9) 
q V 


from which Eq. (8) gives 


Ap 4a Mx — Cpt 3 x — Met 
-— = — tees? ——— 4 e — cos! smc 
x VW ct 


q Br 
(10) 


where the interpretations of (or limitations on) the inverse 


trigonometric terms in Eq. (8) apply. From the considerations 
of the arguments, three distinct time phases discussed physically 
in a subsequent paragraph, can be established: 

Phase 1.—If 


(Mx — cB*t)/x > 1 
then 
(* — Met)/ct > | 


ind the pressure coefficient reduces to 


Ap/q = — 4ta/M 11) 
Phase 2.—lf 
—1 < (Mx — cp*t)/x <1 
then 
—1 < (x — Met)/ct < 1 
and the pressure coefficient is given by Eq. (10) directly. 
Phase 3.—If 
Ux — CS*t)/x S$ — I 
then 
x — MCt)/at < 1 
and the pressure coefficient reduces to 
Ap/g = —4a/B 12 
Physically, Phase 1 is bounded by ¢ = O and ¢ = x/(V + ¢), the 


latter value denoting the time required for a point x on the wing 
to meet an oncoming disturbance originating at the leading edge. 
Phase 2 is bounded by ¢ = x/(V + c) andt = x/(V — c), the latter 
value denoting the time required for a point x on the wing to 
overtake the outgoing disturbance originating at the leading 
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ae W=O m W=-Var * W=0 oz 
V ow wm oe 
sth vv y / vv 
(a) SPATIAL FUNCTION 
5 
w(t) | 
| ————_ 
| | 
91 ‘ 
qT, tz t 
(b) TIME FUNCTION 
Fic. 2. Spatial shape function for wing-penetrating 
sharp-edged gust. 
TABLE 1 
Present note Biot Heaslet & Lomax 
x v + Met x + Mt 
t t t/¢ 
edge. Phase 3 refers to values of ¢ > x/(V — c) and corresponds 
to the steady Ackeret results. 
2) Wing Penetrating a Sharp Vertical Gust.—For a wing 


passing into a sharp vertical gust of vertical velocity wo, the spa- 

tial shape function in terms of the boundary conditions may be 

written (see Fig. 2) 

W(é, 0) = wo, for the portion of the wing within the gust } 
= (), for the portion of the wing outside the gust \ 


> 


The time function w(¢) is used as given in Eq. (8). 
of the velocity potential follows directly from the preceding dis- 


Derivation 


The corresponding pressure coefficient is 


Mx — CB*t 
—— 14 
x 


rhe pressure coefficients in the three time phases previously dis 


cussion 


q rB I 


cussed are: 


Phase 1. 
Phase 2. 
Phase 3.- 


state condition for constant vertical motion 


The pressure coefficient is zero. 
The pressure coefficient is given by Eq. (14) directly. 
The pressure coefficient reduces to that of the steady 


Ap/q = 4wo/BV 15 

For convenience to the reader the relations between the co 

ordinates used herein, those used by Biot,! and those used by 
Heaslet and Lomax are given in Table 1. 

By direct substitution it is easily shown that all the results are 


in agreement. 


REFERENCES 


Wing Striking a Sharp-Edged Gu 
12, No. 5, p. 296, May, 1949. 

of A ir Force 

N.A.C.A 


1 Biot, M. A., Loads 
Journal of the Aeronautical Sciences, Vol 

? Garrick, I. E., and Rubinow, S. I., Theoretical Study 
Oscillating or Steady Thin Wing in Supersonic Main Stream 
No. 1383, 1947 

3 Heaslet, M. A., and Lomax, Harvard, Two-dimensional Unsteady Lift 
Problems in Supersonic Flight, N.A.C.A. T.N. No. 1621, 1948 

4 Campbell, G. A., and Foster, R. M., Fourier Integrals for Practical Appl 
Van Nostrand Co., 1948 


m a Supersont 


on on 
TN 


tions; D Inc., 








NS) 
whil 
in tl 
thro 
that 
entr 


The 


and 1 


in th 
are ¢ 
hy 

by tl 
ber « 


ating 





TTLAX 





esponds 


a wing 
the spa- 
may be 


rivation 
ing dis- 


14) 


isly dis- 


irectly. 
steady 

15) 

the co- 


ised by 


ults are 


ed Gusi 
49. 

es onan 
A. TN 


ady Lif 


1 Applt- 





On Unsteady One-Dimensional Flows 
with Heat Addition 


William Swartz 

Department of Mathematics, Massachusetts Institute of Technology, 
Cambridge, Mass. 

May, 1949 


N INTERESTING PAPER of Kahane and Lees! treated, by the 
A method of characteristics, the problem of addition of heat 
toan initially isentropic flow through atube. They have utilized 
a form of the characteristic equations containing more than two 
dependent variables, necessitating approximation other than 
those involved in taking finite differences. This has 
pointed out by Lin,? who also derives a simple form of the char- 
.cteristic equations permitting the numerical integration to pro- 
At his suggestion, the pres- 


been 


ceed in a straightforward manner. 
ent writer has derived these same equations in a manner anal- 


ogous to that of Meyer.* For a perfect gas, these equations are 


ds/dt = q* along dx/dt = wu (1) 
d du ° dx 
p = at dt along =uFa (2) 
yp a 7 dt 


where s is the entropy divided by the specific heat c, and g* is 
the rate of increase of entropy—i.e., the heat generated per unit 
mass divided by c,T. Other symbols are those used by Kahane 
ind Lees 
In view of the additional approximations used by Kahane and 
Lees, it was deemed interesting to solve the same problem with 
the simplified method. The relationship between the notations 
Q* and K of reference 1 and g* is 
q* = (1/¢e,T) (dQ*/dt) = K/c, (3) 
It is convenient to work with all quantities in dimensionless form 
The value of K employed by Kahane and Lees gives g* (dimen 
sionless) = 0.386y. Dimensionless forms of the variables a, p, 
and p are obtained by dividing by their respective initial values; 
the velocity is to be divided by ao, the initial value of a. The 
distance and time parameters are taken as x/L and aot/L, respec 
tively, where L is the length of the section of the tube which is 


subjected to heating. In terms of dimensionless quantities, the 


Egs. (1) and (2) remain unchanged, but the equation of state be 
comes 
p = pre (4 
and the v« locity of sound is 
f =) 
a= V p/p (oO 


Since the calculations are exceedingly laborious, it seems worth 
while to point out that the solution may be explicitly integrated 
in the region bounded by the x-axis and the Mach lines passing 
through (0, 0) and (1, 0) (Fig. 2 of reference 1). It can be shown 
that p and u are constants throughout the region, and that the 
entropy and pressure are given by 


s = q*t (6 
ind 
p = et! (7) 
The Mach lines are 
x = 0.8 = (2/q*)e1/20") + const. (8) 


ind the path lines are the straight lines 
x = 0.8t + const. (9) 


The calculations are also materially shortened by noting that, 
in the region to the right of x/ = 1, the pressure and velocity 
are constant along each of the family of Mach lines characterized 
by the plus sign. Likewise, they are constants along each mem- 
ber of the minus family of Mach lines in the area to the left of 
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x/L = 0. This is easily shown, remembering that g* = 0 in 
these regions. 

The characteristic and path line diagrams obtained agree well, 
in general, with those of Kahane and Lees. Several Mach lines 
of the minus family at the left of the disturbed area were ex- 
tended until they nearly intersect with the straight Mach line 
through (0, 0). An estimate of the coordinates of the beginning 
of the upstream-moving shock wave was thus found to be in the 
neighborhood of apt/L = 1.2 (presumably more accurate than 
Fig. 3 of reference 1, which shows no shock at aof/L = 1.5 
Qualitatively, of course, there is general agreement. 

Of interest is the behavior of temperature along the path lines 
Using the equation of state and Eq. (5), the absolute tempera- 
ture in dimensionless form is merely 


T = a* 


The temperature was calculated along three path lines originating 
atx/L = 0,0.5,and 1. Although the last path line is in a region 
of no heat addition, a slow steady increase in temperature is 
found. This is associated with the compression indicated by the 
density calculations. Along the other path lines, temperature 
increases at a more rapid rate within the area of heat addition, 
with a slight drop in temperature upon crossing the boundary 
x/L =1. This drop of temperature after an addition of heat to «a 
fluid element is not contradictory to the Second Law of Thermo 
dynamics. It has been mentioned before —e.g., in reference 4 
Foa and Rudinger® emphasized the change of temperature asso 
ciated with a flow field and raised objections against statements 


not made by previous authors 
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The Problem of Stringer Roll 


Paul E. Sandorff 

Group Engineer, Research, Lockheed Aircraft Corporation, Bur-. 
bank, Calif. 

May 23, 1949 


wi REGARD TO THE PROBLEM OF stringer roll treated by 
G. C. Best in his paper in the May, 1949, issue of the 
JouRNAL (Vol. 16, No. 5, p. 289), the writer wishes to direct 
attention to the article ‘Bending Rigidity and Column Strength 
of Thin Sections,’’ A..S.M.E. Transactions, November, 1947. 

The reader will find there a ready-to-use method for treating 
rolling tendencies of unsymmetrical flanges, together with sub- 
stantiating test data and charts that facilitate solution. The 
mathematics begin from the same basis as that used by Mr. Best, 
but they are developed to show that flanges act with an efficiency 
determined only by geometrical considerations. Revised or ef 
fective section properties may therefore be computed which indi- 
cate directly the behavior of the structure under column or bend- 
ing load. 

It is believed that the stress analyst will find this approach 
better suited to his purposes, though Mr. Best's treatment does 
make available certain information otherwise not obtainable 
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Note on Life Expectancy Analysis of Structures 


George Gerard 
Guggenheim School of Aeronautics, 
York, N. Y. 


lune 17, 1949 


New York University, New 


HE CUMULATIVE DAMAGE CONCEPT! has received rather wide 
paisa because it provides a simple method of analysis in 
which the effects of loading at various stress levels can be ac- 
counted for. Experimental evidence as to the validity of the 
cumulative damage concept, however, is somewhat contradictory. 
Fatigue test data obtained by loading at only two stress levels 
indicate that metals may be appreciably damaged by overstress- 
ing and strengthened by understressing. On the other hand, 
experiments by Miner! at a number of stress levels offer data in 
support of the cumulative damage concept. 

It is to be noted that relatively little attention has been given 
to the frequency distribution of loads imposed on an airplane in 
considering the validity of the cumulative damage concept. It 
was found that, by considering both the loading spectrum and the 
fatigue properties of the material, the distribution of damage in- 
crements was of the nature of a probability curve. Hence, it is 
proposed that the validity of the cumulative damage concept 
may be established for life-expectancy analyses by using a similar 
loading spectrum in laboratory tests. 


LIFE-EXPECTANCY ANALYSIS 


For particular airplane characteristics, the loading spectrum 
usually specified in the form of load relative frequency as a func- 
tion of load factor can be assumed to be a load relative frequency- 
stress level relationship for the structure. The fatigue property 
spectrum may also be represented as a single curve when drawn 
for a mean stress corresponding to the pertinent load factor. 
Methods for constructing such curves are presented in reference 
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Fic. 1. Life expectancy analysis using cumulative damage con 


cept. 


The analysis of life-expectancy proceeds from the loading an 
fatigue property spectra. The criterion for fatigue life according 
to the cumulative damage concept! is 

=n/N = 1.0 


In the analysis, the loading spectrum is given in terms of relative 
frequency and contains the fatigue life, F, as a factor. Hence, 
from Fig. 2, the per cent of fatigue life used at various stress levels 
contains the F term which when summated has the following 
form 

SD] 


1/F)S (n/N) =6 2 


In general, 6 ~ 1.0, as required by the cumulative damage con- 
cept and, hence, by adjusting the value of F, the condition im- 
posed by Eq. (1) can be satisfied. Thus, the fatigue life, F, in 
number of repeated loads encountered is obtained as 







































































Fic. 2. Distribution of damage increments for analysis of reference 2 
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DAMAGE INCREMENT DISTRIBUTION 


From fatigue analyses that utilize the cumulative damage con- 
cept, the increments of damage, n/N, which are attributed to 
each stress level, can be plotted to obtain a damage curve such 
as shown in Fig. 1. .Damage curves plotted from the gust analy- 
ses conducted on an hypothetical airplane in reference 2 are shown 
in Figs. 2. Analyses on several other airplanes indicated curves 
of a similar nature and, consequently, those shown in Figs. 2 may 
possibly be representative of this type of analysis. 

It is significant to note that the curves all peak sharply, indi- 
cating that the fatigue life is strongly influenced by a particular 
stress increment referred to herein as the ‘‘critical damage incre- 
ment.” In fact, analysis of the data used to construct Figs. 2 
indicated that, in all cases analyzed except one, over 70 per cent 
of the fatigue life could be attributed to the two increments at 
the peak of the distribution curve. This trend remains to be 
established for many more cases, but, if it is correct, it indicates 
that almost all damage can be attributed to a rather narrow band 
of stress, the critical damage increment. This trend approaches 
the type of fatigue analysis in which only one stress level is as- 
sumed responsible for fatigue failure. 

Of greater significance than the critical damage increment is 
the fact that the distribution of damage increments closely ap- 
proaches a class of curves known as probability curves. A dis- 
tribution of this nature tends to minimize the effects of damaging 
and strengthening at high and low stress levels, respectively, since 
the various stresses are disposed about the critical damage incre- 
ment and of a low magnitude relative to the critical increment. 

It appears that, because of the probability distribution of the 
damage increments, laboratory tests should be run using a loading 
spectrum of this nature. Tests run in this manner may establish 
the validity of the cumulative damage concept for the particular 
type of life-expectancy analysis discussed herein, although it is 
probable that its validity cannot be established in general. 
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End-Fixity of Columns 


George Winter 

Professor and Head, Department of Structural! Engineering, 
Cornell University, Ithaca, N.Y. 

May 13, 1949 


I THEIR PAPERS ON end-fixity of columns in the February, 
1949, JoURNAL OF THE AERONAUTICAL SCIENCES (Vol. 16, 
No. 12, pp. 116, 117-119), N. M. Newmark and G. C. Best have 
given attractive approximate formulas for determining the ef- 
fective length coefficient of columns with known rotational end 
restraints 

It may be of interest, in this connection, to note that simple 
charts for the determination of this same quantity were published 
in the following: ‘‘Elastically Encastred Struts,’’ by N. J. Hoff 
J. Royal Aero. Soc., Vol. 40, p. 663, November, 1936), and 
Buckling of Trusses and Rigid Frames,’’ by George Winter, 
P. T. Hsu, B. Koo, and M. H. Loh, (Bulletin No. 36, Cornell 
University Eng. Exp. Station, April, 1948). These charts allow 
the determination of c directly on the basis of the rotational 
spring constants at both ends—i.e., Dr. Newmark’s & and ke 
or Mr. Best’s K,; and K.—without the necessity of first computing 


iv) 
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from them auxiliary quantities, such as in Mr. Best’s method 
Ways of accounting for the influence on these spring constants 
of axial forces acting in the restraining members are simplified 
by additional charts in the second of the above publications. 
The main chart of that paper for determining c accounts not only 
for positive but also for negative restraints. The latter feature 
is important because a complete stability analysis of a frame- 
work is frequently impossible without accounting for negative 
restraints of which Mr. Best notes that ‘they are not considered 
to actually occur.” 

Indeed, in the simplest truss, consisting of two compression 
members rigidly joined at their intersection and for simplicity’s 
sake hinged at their other ends to the third, tension member nega- 
tive restraint does occur for the compression bar with the smaller 
L/r at the end joined to the one with larger slenderness if axial 
forces are equal in both of them. In other words, the more 
slender bar, tending to buckle, induces an end moment into the 
more rigid bar rather than to provide a restraining moment. 

The importance of this factor is brought out in Mr. Best’s 
papers, in which is shown the application of these charts to an ex- 
act solution, by successive approximation, of the stability of a 
rigid framework composed of any number of bars. 


Errata—The Rolling and Yawing Moments 
Due to Yaw of Unswept Wings 


R. F. Anderson 
Northrop Aircraft, Inc., Hawthorne, Calif. 
July 14, 1949 


—— ERRORS appeared in my Readers’ Forum item on pages 
445-447 of the July issue of the JouURNAL OF THE AERONAU- 
TICAL SCIENCES: 

(1) Page 447, first paragraph after Eq. (5). The last sen- 
tence of this paragraph should read: ‘‘Values of m were ob- 
tained from reference 4 using mio = 27.’ 

(2) Page 447, third paragraph after Eq. (5), fourth and fifth 
The corrected sentence should read: ‘‘For example, the 
o0C,? = —0.125/A could be used.” 
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Stability of Sandwich Plates 


P. P. Bijlaard 

Professor, Institute of Technology, Delft, Holland, on leave; 
Associate Professor, Cornell University, Ithaca, N.Y. 

June 15, 1949 


[Ts AN INTERESTING PAPER, Seide! calculated the compressive, 
buckling stress of sandwich plates with clamped unloaded 
edges, assuming the faces to act as membranes. His computa- 
tions may be used as a check on a method the present writer de- 
veloped some time ago.” * In order to do this the thickness of 
the faces has to be taken into account in Seide’s calculations too, 
since, apart from disregarding the flexural rigidity of the faces, 
its neglect underestimates the shearing rigidity of the core. With 
a thickness ratio t;/h. = 1/20 of faces and core, it already results 
in a 9 per cent smaller shearing rigidity. 

The thickness of the faces may be taken into account by insert- 
ing other values for the shearing stiffnesses Dg; and Dg, in Eqs. 
(6) of reference 4. As was derived in references 2 and 3, think- 
ing all other deformations to be zero and considering the faces 
to have no flexural rigidity, as in reference 4, but having their 
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actual thickness, with the notations of reference 1 the shearing 
force per unit width is 


QO, = G.[(he + ty)?/h-] (Ow/Ox ] 
so that 
OC. . (h + ty)? 
De, = = ( = JD (2) 
”" ww / OX ” h. ” 


This leads for example in the elasti¢ range to the differential equa- 


h. D ZWyocr’ Ow , 
Vw +(1—- — —V? - = () (3) 
(he + ty)? G. D. ox? 


tion 


instead of Eq. (A5) of reference 5. In these formulas, G, is the 
shear modulus of the core material and D and o,,’ are the flexural 
stiffness and buckling stress of the sandwich plate, both under the 
assumption of zero flexural rigidity of the faces. In the plastic 
range, the more exact differential equations are also obtained by 
replacing value h, by (h, + t;)?/he. 

The total critical thrust 2t;o., is equal to the sum of the thrust 
2tyoer’ from Eq. (3) and that of the single faces, 2t;o0, calculated 
for the same wave length, analogous to what was done in refer- 
ences 2 and 3, and has to be minimized. 

In references 2 and 3, the present writer proved that the critical 
load for a sandwich plate may be written as follows: 

Po = Pe + (PaW + P,-?)-? (4) 
in which Py is the joint critical load for the single faces; Pm is 
the critical load for the sandwich plate, if G. with respect to the 
2 and ry, is assumed to be infinite, and P, is ob- 
tained by assuming only G, to havea finite value. In calculating 
P,, and P,, the proper flexural rigidity of the faces is neglected. 
All loads are calculated for the same wave length, which has to be 


shearing stresses 7 


chosen such as to make P?,, a minimum. 

This method can be considered as exact if the deflection sur- 
faces applying for the computation of Ps, ?,,, and P, are identi- 
cal, as for example for cylindrical buckling of a simply supported 
sandwich plate under end thrust and for buckling of a rectangular 
sandwich plate with simply supported edges, subjected to com- 
pression in both directions.2 * If the deflection surfaces have 
not the same form, the method yields conservative results, be- 
cause in reality they have to coincide everywhere, giving an ex- 
cess restraint and so increasing P,,. The latter case obtains 
for the case considered by Seide in reference 1, so that, taking 
into account the changes proposed above here, Seide’s paper 
provides a possibility to check how much the present writer’s 
method is conservative here. Even for this case, where the de- 
flection surfaces for computing P, and /P, differ greatly,’ it 
proves to be only conservative by, at most, 4 per cent. Thus it 
may be advised to use it also in other cases, since the calculations 
involved are much less than those based on the differential equa- 
tions of the plate as a whole. 

With the notations of reference 1 the method gives an example 


for a rectangular plate under pure compression or pure shear 
(P.- or See) = nkor?D, h2 (5) 


in which 7 is the reduction coefficient for plasticity of the faces, 
ko is the elastic buckling stress coefficient for a plate with infinite 
shearing rigidity as far as stresses 7,;, and 7,,. are concerned, and 


bis the plate width. In Eq. (5), 


D. = (D;/De + ¢)D: (6) 
in which 
D, = Eyt;(he + ty)? 211 -— uy?) 
D; = Eyt?/6(1 _ uy?) (7 
p = a/(a + nrko) 7 
r = (#*D,/b*G.)/[he/(he + ty)*] 


while Ey and yy are Young’s modulus and Poisson’s ratio for the 
faces. For loading by pure compression and with simply sup- 
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ported or clamped edges 
a = 1 + (B?/m?* 8 


in which £ is the plate aspect ratio a/b and m is the number of 
half wave lengths in the direction of loading. For pure shear in 
a long plate with simply supported or clamped edges 


a = [1 + (L2/b*)]"? (9) 


in which L is the half wave length in the long direction. The ey. 
pression (D;/D; + ¢) in Eq. (6) has to be minimized by variation 
of the wave length. 
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On Damping in Pitch for Delta Wings 


John W. Miles 
Department of Engineering, University of Calitornia at Los Angeles 
May 24, 1949 


magic HAS RECENTLY CALCULATED the damping in pitch 
for thin triangular wings on the assumption of an equiva 
lent steady flow over a parabolically cambered wing.' Using his 
results, we have obtained a quasi-stationary solution including 
all first order terms (in frequency). 
following relation: If # is a linear operator for a given super 
sonic plan form which gives the potential (#) at the wing due to 
the downwash (w) there [viz. @(x, y) = &{w(x, y)}], om the as 
sumption of linearized steady flow, the quasi-stationary solution 


This solution is based on the 


is given by? 
P(x, y, t) = &) wer, y, t} — sec? 6[0/d Ut)| [xb{ w(x, ¥, bj - 
b{xw(x, y, )}] + O[(wb/U)? 


where U is the free-stream velocity, directed along the positive 2 
axis; 6 is the Mach angle; and w is the angular frequency 
(This result may require modification in those configurations 
where a Kutta condition must be invoked.) 

If we consider a pitching oscillation of a delta wing about the 
centroidal axis [x = (2/3)c, where c is the maximum chord), 
the downwash distribution is given by 


4g jot 2 
w(x, y, t) = Ua ‘ + aa x— 7 exp (twt) 2 
4 Sf 0 


where we choose the average chord (c/2) as our characteristic 


length. The lift and moment [about x = (2/3)c] coefficients 
due to this distribution are found to be 
: pU*\_|-1oL 
Cy, = —— IS =m oli 
= 2 Oa 
= mf, «<1 3b 


= 2rtany + O0(o3), o<1 (8 
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Si 
Oo. 3 tan y + Ofo' logo), c< 1 (4c) 
Cu, = () (9) 
oV?\ cl! oOoM m Q 
Cy. = - 1S — = —— (2 —tan’?@), o > 1 
. 2 2 Tart a a) 
re) > . 
2l (6a) 
m : ' : : , 
—. (2(2f1 — fo) — (8fo — 2f;) tan? @], o <1 (6b) 
if . . a 
8x : ; 5 
“a tan y + O(e* logo), o< 1 (6c) 
m = 4tané 4) 
: (r/2)o 
fe =feilo) = > (8) 
E’(o} 
: ‘ (3r/4)o(1 — o?) 
fi = file : -_ : (Q) 
. ; (1 — 267)E’(c) + o?K'(c) 
o = (tan y/tan @) (10) 


where 2y is the apex angle, m is the two-dimensional lift curve 
slope, and K’ and E’ 
integrals of the first and second kind, respectively. 


are the complementary complete elliptic 
The cases 
7>1,0<1,and¢ < 1 correspond to the cases where the leading 
<dges are outside, inside, and close to the axis of the Mach cone, 
respectively. 

If the axis of rotation is at x = ac/2, we have 


Crg(a) = Cr4(4/3) — [a — (4/3)] Cr, (11) 


(4/3) — 


[a — (4/3)1C,, (12) 


Cug(@) = Cuy(4/3) + fa — (4/3)1Cr, 
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Upon examination, we find that below a certain critical Mach 
Number Cug will be positive (i.e., ‘‘negative’’ damping) for axes 
located between two points lying fore and aft of mid-chord (a = 
1). At the Mach Number specified by 


tan@d@=1, M= V2, (13a) 


4f, — 3fo\'”? 
tan @ = 
3fo — 2h; 


these two points coincide at a = 1, and for greater Mach Num- 


(13b) 


bers, Cay remains negative. The maximum value of this critical 


Mach Number is V 2(¢ = 1). We also find that Cy is negative 
for all axes aft of a = 3/2, regardless of apex angle. 

In the limiting case of the extremely slender wing, (o < 1), we 
find 


“ ) = —Pri(D 
( Mg (a ae 


— a)? tan y + Oe log a (14) 
so that we have the rather interesting anomaly that the damping 
vanishes for rotation of an extremely slender delta about the base. 
We remark that this same situation holds for any slender body. 

The functions f; and fy are given by the curves mg[0C,/- 
O(cqg/V)] and (O0C;,/da), respectively, in Fig. 5 of Robinson’s 
paper,' where. represents o 
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